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SUMMARY 


Computer programs have been developed to implement the computational scheme 
arising from Van Holten's asymptotic method for calculating airloads on a helicopter 
rotor blade in forward flight, and a similar technique which is based on a discretized 
version of the method. The basic outlines of the two programs are presented, followed 
by separate descriptions of the input requirements and output format. Two examples 
illustrating job entry with appropriate input data and corresponding output are included. 
Appendices contain a sample table of lift coefficient data for the NACA 0012 airfoil 
and listings of the two programs. 


INTRODUCTION 

The computer programs described in this report were developed during the course 
of an evaluation of Van Holten's asymptotic method (ref. 1) for the calculation of 
airloads on a helicopter rotor blade in forward flight. The validity and computational 
feasibility of the approach were investigated (ref. 2), and numerical results for specific 
flight conditions were compared with corresponding experimental data and 
computations based on other analytical methods. Program ASYMP1 was written to 
implement the computational scheme of reference 1 (the relevant equations and 
expressions are given in ref. 2). As an extension of this investigation, the above 
computational scheme was made more efficient by discretizing the variation of the 
doublet strength distribution, g(rj 3 > Y^), utilizing both piecewise constant and piecewise 
quadratic representations for the spanwise variation. The details of the discretized 
scheme are presented in reference 3 and the corresponding computational method is 
applied in Program ASYMP2. 

The general organization is similar for both programs, since the two schemes 
differ only in the manner in which the velocity induced by the blade pressure field is 
calculated. The basic unknown to be determined is the doublet strength function, 
g(r. , Y. ). Its solution is effected by a collocation technique. In the original scheme 
(ret. 1/ this unknown function appears as a continuous modal representation for the 
spanwise variation and a finite Fourier series for the azimuth variation. In the 
discretized scheme the unknowns are the values of g at the midpoints of the spanwise 
segments, with a finite Fourier series for the azimuth variation of each. In either case, 
the solution reduces to the determination of the coefficients in the collocation 
representations, which is accomplished by setting up a system of simultaneous 
equations. 

The main programs and the various subprograms will be discussed in the next 
section, but the general sequence of program steps is as below. 

(1) Read and write input; define auxiliary parameters required for the 
computation. 

(2) Start loop for collocation points. 

(3) Test the tangential velocity U T ; if U T <( u T ) min > set U P the condition for 

zero lift and go to the end of the collocation loop. !i ’ 

(4) If airfoil data tables are used, determine the lift curve slope for the current 
collocation point. 

(5) Start loop for the number of blades and define the first azimuth interval. 



(6) Compute the induced velocity contributions for the current interval and add 
to the corresponding coefficient matrix elements in the system of equations. 

(7) Increment the azimuth interval; if the azimuth limit has been reached, go to 
the the next step; if not, go to Step 6. 

(8) End loop for number of blades. 

(9) End loop for collocation points. 

(10) Set up the extra equations for the blade motion parameters. 

(11) Solve system of equations and write solution. 

(12) Compute and write output. 

Some mention must be made of the way in which the airfoil data are used in the 
computational scheme. The basic equation to be set up at any collocation point is of 
the form 


w b (r bo’ 'W * w ind 


(1) 


where w^. is the velocity induced by the pressure fields of all the blades and w^ the 
normal velocity due to blade motion at a collocation point (r^ , ). This is rewritten 

as ° 1 ° 

w b<W V - '»in/ D * [*hd - (w in/ D ] 


= < w ind> 2D + 


(2) 


x2D 


where (w. .) is the induced velocity corresponding to a steady, two-dimensional flow. 
Accordingly, corrections based on airfoil data are only applied to this term and, after 
modification, equation (2) becomes 


2 n 

w, (r, , T. ) = (w. ,) . + AW 

b bo bo ind mod 


■ [ <w ind )2 mod - <w ind )2D ] 
For steady, locally two-dimensional flow, 


+ w 


ind 


(3) 


<w ind> = -S< r bo’' l 'bo )/pU T (W 

At the local incidence and Mach number corresponding to the collocation point, the 
airfoil data are interpolated in the form of a locally linear relation given by 


C H = a “ + C lo 

To account for such a relation, the expression in equation (4) is modified as 
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( 5 ) 
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Using equations (4) and (5) in equation (3), the modified form of the boundary condition 
is obtained as 


w b (r bo- V 


(2 it 

a 
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- U 
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+ w. 


ind 


( 6 ) 


Comparing equations (1) and (6), it can be seen that the modification due to the airfoil 
data is purely an additive term that is easily incorporated into the basic computational 
scheme. 

The next section presents an outline of the main program and the various 
subprograms for ASYMP1 and ASYMP2, listing the input and output of each subprogram 
along with a brief description of its function. The sections following this deal with 
descriptions of the input and output for the two programs, and samples of job entry with 
the corresponding output. 


PROGRAM OUTLINE 
Program ASYMP1 

Main Program . - The program steps listed in the previous section will now be 
discussed in more detail. 

(1) The details regarding input data and the format in which they are to be 
entered are given in the next section. The data include parameters describing rotor 
geometry, the flight condition and blade motion. Specification of the blade motion 
parameters (collective pitch, coning angle and the two frst harmonic flapping 
coefficients) is optional. If they are not specified, they will be calculated as part of the 
solution by generating additional equilibrium equations. The input also includes the five 
spanwise locations of the collocation points, the normal and reduced azimuth intervals 
to be used for the numerical integration, the minimum value of the local onset velocity 
below which the zero lift condition is to be used and an integer specifying whether 
airfoil data are to be used. The program then writes the input data as part of the 
output and defines auxiliary quantities such as the induced velocity from simple 
momentum theory and certain factors occurring in the induced velocity contribution 
from the trajectory segment immediately adjacent to the collocation point. 

(2) The collocation loop consists of an outer loop for the eleven equally spaced 
azimuth locations and an inner loop for the five spanwise locations. 

(3) For the current collocation point, the local onset velocity (UT) is compared 
with the specified minimum value (UTMIN). If UT> UTMIN, the next step is executed. 
If not, the zero lift condition is set up and the next collocation point is taken up. 

(4) If airfoil data are not used, this step is skipped. If they are used, the local 

values of a/2 tr and are interpolated from the tables. The local values of Mach 

number (MLOC) and incidence (ALOC) are first defined. Subroutine TABSCH is called 
to find the position of MLOC and ALOC in the arrays MCL and ACL that were read as 
part of the input. The required values of a/2ir (SLCR) and C& (CLQ)) are then 
determined by linear interpolation. J 

(5) As mentioned in reference 2, the near field has a square root singularity 
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corresponding to the leading edge, and this is dealt with by stopping the numerical 
integration just in front of the leading edge and accounting for the remainder 
analytically. Before starting the numerical integration, therefore, this contribution is 
entered into the coefficient matrix, A. The system of equations is of the form 

Ax = B 


where x is the array of unknown coefficients (A.., B.., eq. (9) of ref. 2 X If 
UT <UTMIN, the zero lift condition is set up at this point. ‘ 

In order to start the numerical integration, a loop for the number of blades is set 
up. For the current blade, Subroutine DMIN is called to determine the positions along 
the fluid particle trajectory at which it is directly over the blade, within a distance 
DMAX. Around these locations (stored in array PMIN), the reduced azimuth interval 
DP2 will be used. The first azimuth interval is defined with its ends at PB1 and PB2. 

(6) For the current azimuth interval, a loop is set up for the 5-point Gauss- 
Chebyschev integration (p. 12 of ref. 2). With 55 unknown coefficients, 4 blade motion 
parameters and one right-hand side, there are 60 integrations to be carried out over 
each interval of azimuth. The corresponding functional values are sequentially obtained 
from function FUN2. 

(7) PB1 and PB2 are incremented for the next azimuth interval. If the azimuth 
limit PLIM has been reached, the integration is stopped and the next blade is taken up. 
If not, the interval is tested to check if it includes one of the "close" locations stored in 
PMIN. If it does, the reduced azimuth interval is used. 

(8) At the end of the loop for the number of blades, the terms in the blade 
normal velocity, w b (r bo , ’i' ), are entered in the corresponding coefficient matrix 
elements. 

(9) At the end of the loop for the collocation points, certain spanwise integrals 
required for the total blade lift and the moment about the hub due to the lift are 
calculated (p. 42 of ref. 2). 

(10) To close the system of equations, four additional equations are set up (eqs. 
(E8), (E9), (E10), (Ell) of ref. 2). If the four blade motion parameters are specified in 
the input, these equations are replaced by equations of the form 


0 

o 


(0 ). + 

o input 


= (a ) . etc. 
o input’ . 


(11) The completed system of equations is solved by calling Subroutine GELIM. 
This is a library-supplied routine that uses the LU decomposition. The solution is 
overwritten on the vector B and the program prints out the values of the collocation 
coefficients and the blade motion parameters. 

(12) With the basic solution complete, the program oomputes and prints out 
various output quantities in tabular form. The output format is described separately in 
a later section. 


Subroutine DMIN. - 

Input: MU(y), LAM(X), DB(AT.), RB0(r bQ ), PBfl'F^), PLIM, DMAX 

Output: I, P 
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Comments: This subroutine locates those azimuth positions along the fluid particle 

trajectory at which the particle is "close" to the j n blade. The axis is fixed to the 
blade and rotates with it (fig. 1 of ref. 2) so that the x. coordinate along the trajectory 
periodically goes to zero whenever there is an intersecuon with the z^ axis. At some of 
these locations the particle may be too close to the blade (within a specified distance 
DM AX). These are the locations (I in number) that are returned in the array P. The 
routine locates the positions by scanning the trajectory with small azimuth increments 
checking for a change in sign of the x. coordinate. When such a position is located, the 
distance from the blade is compared with DMAX. 

Subroutine TRAJ . - 

Input: R0(R /Rj), AR(A), MU, LAM, RB0, PB0, DB, PB^) 

Output: R(r), SX(sin \), CX(cos x)> SHP(sinh T), CHP(cosh T), 

ST(sin 0), CT(oos 0), SHE(sinh n)> CHE(cosh n), 

SP(sin 4>), CP(cos <J>) - returned through common block TRAJ1 
ZS(z b /s) 

Comments: Given a point ^ on the trajectory relative to the j** 1 blade, the 

corresponding coordinates in various coordinate systems are calculated. The 
coordinates (r,x, z b ), (T,0,x) and (n,<|>) are respectively of the point in cylindrical, 
prolate spheriodal and plane elliptic coordinate systems. 

Function FUN1 . - 

Input: R0, AR, I 

N(n), PI(kt) - through common block MAIN1 

R, SX, CX, SHP, CHP, ST, CT, SHE, CHE, SP, CP - through 
common block TRAJ1 

Comments: The values of various functions required in FUN2 are calculated. 

Specifically, the six derivative expressions to (Appendix C of ref. 2) are returned 
for values of I from 1 to 6 in the input. 

Function FUN2 . - 

Input: R0, AR, TW( e), MU, LAM, RB0, PBQ), DP(A¥ b ), DB, 

X(\, P- 12, ref. 2), I 

N, PI - through common block MAINl 

Comments: As described in Step 6 of the main program, the coefficient matrix 

elements corresponding to the various unknowns and the right-hand side require 
numerical integration, for which the necessary functional values are returned by 
function FUN2. Given the azimuth position X, Subroutine TRAJ and Function FUN1 are 
used to set up the relevant expressions. The value, 1 = 1, corresponds to the function 
multiplying A Qo (eq. (9) of ref. 2). The value, I = 2, corresponds to the function 
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multiplying A (n = 1,2, 3, 4). The values, I = 3, 4, 5, 6, 7, respectively correspond to the 
functions multiplying 0 , a Q , b, and the right-hand side. The forms of these 
integrands are given in Appendix C of reference 2 as induced velocity coefficients. The 
nonintegral parts of these expressions are the result of analytical integration of the 
near field over a small interval AT. adjacent to the collocation point and are defined in 
the main program (Step 5). Over tnis interval, therefore, the near field contribution is 
skipped in Function FUN2. 

Function FUN3. - 


Input: R0, AR, X, I 

N, PI - through common block MAIN1 

Comments: The expressions for the total bladp l^t ^nd^the moment about the hub due 
to the lift involve certain spanwise integrals (I , I , I , I (n = 1,2, 3, 4) see p. 42 of ref. 2) 
and the corresponding integrands are set up inrl^Ni! fhe radial ^os^ioy, r^/Rj, is X. 
The notation, I = 1,2, 3, 4, corresponds respectively to the integrals I Q , I Q , I n , l£ 

Function PNM. - 


Input: N(n), M(m), X(x) 

Comments: This function generates the associated Legendre function P m (x) over the 
ranges 0 <n *4, 0 £m<3, |xl<l. Although the relevant recursive relations could be 
used, the function defines P™ explicitly in terms of x for all the above values of n and 
m (Appendix D of ref. 2). 

Function QNM . - 

Input: N, M, X 

Comments: The associated Legendre function Q^Oc) is calculated over the ranges 

l<n<4, l<m^.2, lx I >1. As listed in Appendix D of reference 2, the exact definitions 
are used for |x |<3 and asymptotic expansions are used for |x | >3 to avoid the 
accumulation of roundoff error. 

Subroutine TABSCH . - 

Input: X, N, XT 

Output: 11,12, INT 

Comments: This routine searches an array X, of dimension N, for the position of a 

value XT. If XT lies between X(I1) and X(I2), INT = 0. If XT is outside the range of X, 
the subroutine returns INT = -1 for XT<X(1) and INT = 1 for XT >X(N) (it is assumed 
that the elements in X are arranged in increasing order). 

Program ASYMP2 

Main Program . - Here again the previously listed program steps will be discussed 
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in more detail. 

(1) This step is similar to Step 1 of Program ASYMP1. There are some changes 
in the entry of input data (described in the next section). Instead of spanwise 
collocation point locations, the program reads the spanwise locations at which the blade 
is divided into segments for discretization. Collocation points are located at the center 
of each segment. The integer, ISEL, specifies whether a piecewise constant or 
piecewise quadratic representation is to be used for the computation. Following the 
reading and writing of input data, some arrays that will be required later are set up. 
Arrays FX1 and FX2 contain the average values of the factors 



^ R 1 ~ R 0> Vl - x 2 ) 
2Rj A j 


over each chordwise segment. These values will be used in the near field calculation. 
Array GF relates the endpoint and midpoint values for the spanwise segments in the 
piecewise quadratic representation. If p-(j =1,. . j6) represent the ends of the 5 
spanwise segments (Pj = R q , Pg = Rj) and r j(j = 1,. . 5 ) the midpoints, then 

5 

g(pj) = £ (GF) . g(r.) (i = 2, 3,4, 5) 

1 j=l 13 J 


(2-5) These steps are generally the same as the corresponding steps in Program 
ASYMP1. 

(6) For the current azimuth interval, the slopes and intercepts for the linear 
approximations to x. (T, ), y^( T^) and r, (T^) are defined as XBI, XBS, YBS, RBI, RBS. 
The far field contribution is computed using Subroutine FFINT. For computing the 
common part and near field contributions (which are both dependent on the local 
spanwise dipole strength) it is necessary to divide that part of the interval (PB1, PB2) 
which has the trajectory within the blade span into subintervals such that each 
subinterval has a trajectory wholly within one spanwise segment. This is done by calling 
Subroutine SUBIVL. The common part and near field contributions are then computed 
respectively by calling Subroutines CPINT and NFINT, summing the contributions over 
each subinterval. 

(7-12) Comments on these steps are the same as for Program ASYMP1. 


Subroutine SUBIVL. - 


Input: RBl(r bl ), RB2(r b2 ), PBlO^), PB2(H' b2 ), R(p.,j = 1,. . £) 

Output: I, PI, P2 

Comments: As explained in Step 6 for the main program, the azimuth interval must be 
subdivided so that the trajectory segments within each subinterval lie completely within 
one spanwise segment. This is done by comparing the endpoints for each segment (p., 
p. + j for the j segment) successively with the spanwise coordinates at the ends of the 
interval, r b j and r b2 . The terms I, PI and P2 are all arrays of dimension 5, 
corresponding to the u number of spanwise segments. For the segment J, I(J) = 1 or 0 
depending on whether a portion of the trajectory does or does not lie within that 
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segment. If 1(3) = 1, the azimuth positions corresponding to the ends of that portion of 
the trajectory segment are stored in P1(J) and P2(3). 

Subroutine NFINT. - 

Input: PI, P2, X, ISEL 

XI, XS, YS, RI, RS - through common block MAIN1 

Output: Tl, T2, T3 

Comments: This subroutine computes various integral terms required for the near field 
contribution over the azimuth interval (PI, P2), with X being the chordwise location of 
a segment over which the surface pressure has been averaged. The terms Tl, T2, T3 
correspond to the terms n , n^, x\^ in Appendix C of reference 3. For the piecewise 
constant representation, ISfiL = 0, only Tl is used. 

Subroutine CPINT . - 

Input: PI, P2, ISEL 

XI, XS, YS, RI, RS - through common block MAIN1 

Output: Tl, T2, T3 

Comments: The terms required to set up the common part contribution over the 

azimuth interval (PI, P2) are calculated. The terms Tl, T2, T3 correspond to Cp C 2 , 
C, in Appendix C of reference 3. For the case ISEL = 0, only the term Tl is relevant. 
The common part is singular at the collocation point, although the complete pressure 
field is regular due to the singularities in the far field and common part cancelling out 
in the limit. Since this has been established, when P2 = 0 (corresponding to the 
collocation point) the routine sets all the terms to zero. 

Subroutine FFINT . - 

Input: PI, P2, R, ISEL 

XI, XS, YS, RI, RS - through common block MAINl 

Output: Tl, T2, T3 

Comments: This routine calculates the terms required for the far field contribution 
over the azimuth interval (PI, P2) with R being a boundary of one of the spanwise 
segments (p ., j = 1, . . j6). The expressions for Tl, T2, T3 are derived in Appendix C of 
reference 3 ) As pointed out above, the singularity in the far field gets cancelled out in 
the limit. However, unlike the common part, there is a finite residue left over after 
cancellation. 

In addition to the subroutines listed in this section, both programs make use of 
Subroutine GELIM to solve the system of simultaneous equations. This routine uses 
direct Gaussian elimination with pivoting and details can be found in the Langley 
Computer Programming Manual. 
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DESCRIPTION OF INPUT 


The input data required for the programs can be divided into two parts. The first 
consists of data pertaining to rotor geometry, flight condition, blade motion and some 
additional parameters relevant to the computational scheme. These are assumed to be 
a part of the INPUT file; that is, they must be entered after the job control statements 
when the job is submitted for execution. The second part consists of airfoil data (if 
required) and these are assumed to reside in a file called AFDATA which must therefore 
be available (if the data are to be used) when the job is executed. The format for the 
airfoil data is described separately following descriptions of the basic input for the two 
programs. 


Program ASYMPI 

The READ statements for the first part of the input are given below, followed by 
explanations of the data items. 


READ (5,*) R , AR, NB, TW, MU, ALR, CT, MINF 

READ (5,*) Nl, N2, N3, W 

IF (Nl. EQ. 1) READ (5,*) THC 

IF (N2. EQ. 0) READ (5,*) GAMA 

IF (N2. EQ. 1) READ (5,*) A0 

IF (N3. EQ. 1) READ (5,*) A1 

IF (N4. EQ. 1) READ (5,*) B1 

READ (5,*) (RB0 (I), I = 1, NSP) 

READ (5,*) DP1D, DP2D, UTMIN, NAFD 


Each READ statement corresponds to a line of data input. It may be noted that all the 
above statements specify free format for the data entry so that the different items in a 
single line of data can be entered in any convenient format, separated by commas. 


R0 = 

root radius/tip radius 

AR = 

aspect ratio 

NB = 

number of blades (integer) 

TW = 

built-in linear twist ( e , ref. 2) 

= 

pitch angle at root minus pitch angle at tip 

MU = 

forward speed/tip speed (floating point) 

ALR = 

inclination of tip path plane to flight path, in degrees 
(forward tilt positive) 

a • a 

CT = 

y y 

rotor thrust coefficient = T/p (ttRj ) (flRj) 
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MINF = 
N1,N2,N3,N4 


THC = 
GAMA = 

A 0 = 

A1 = 
B1 = 
Note: 
RB0 = 

DP1D, DP2D 

UTMIN = 

Note: 
NAFD = 


Mach number corresponding to forward speed (note: this is used only 
for interpolating from airfoil data; if airfoil data are not used, this 
item is not needed and can be set to zero) 

= integers associated with the four blade motion parameters (THC, 
A 0, Al, B1 respectively) - if a blade motion parameter is to be 
specified in the data, the corresponding integer is set to 1; if it is 
to be calculated by the program, the integer is set to 0. 

pitch angle at blade root in degrees 

coefficient representing blade flapping inertia 

2 it (air density) (chord) R, /(mass moment of inertia of blade about 
flapping hinge), to be specified if the coning angle is to be calculated 

blade coning angle 

first harmonic longitudinal flapping coefficient 
first harmonic lateral flapping coefficient 
Flapping angle = A0 - Al cos - B1 sin 

array (of dimension 5) containing the spanwise locations of the 
collocation points, as fractions of the tip radius, e.g., 0.3, 0.5, 0.75, 
0.85, 0.95 

= normal and reduced azimuth intervals to be used for the 
numerical integration (step 6 of the main program) in degrees, 
e.g., 15.0, 5.0 

minimum value of local onset velocity at a collocation point for using 
the normal velocity boundary condition (see step 3), as a fraction of 
the tip speed 

Local onset velocity at (r^, ) = (r^/Rj + ** sin ^bo^ 

integer related to use of airfoil data: 

NAFD = 1 - airfoil data used 
NAFD = 0 - airfoil data not used 


Program ASYMP2 

Much of the input is identical to that for ASYMP1, as can be seen from the READ 
statements below. 


READ (5,*) R0 , AR, NB, TW, MU, ALR, CT, MINF 

READ (5,*) Nl, N2, N3, N4 

IF (Nl. EQ. 1) READ (5,*) THC 

IF (N2. EQ. 0) READ ( 5 ,*) GAMA 

IF (NZ EQ. 1) READ (5,*) A0 

IF (N3. EQ. 1) READ (5,*) Al 
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IF (N4. EQ. I) READ (5,*) B1 
READ (5,*) (R (I), I = 2, 6) 

READ (5,*) DP1D, DP2D, UTMIN, ISEL, NAFD 

R = array (of dimension 6) containing the spanwise locations (as fractions 
of the tip radius) marking the division into 5 spanwise segments 

(R (1) = Rtf, R(6) = 1.0) 

e.g., R (I), I = 2, 6 — ► 0.5, 0.7, 0.8, 0.9, 1.0 

ISEL = integer specifying choice of piecewise representation: 

ISEL = 0 - piecewise constant representation 
ISEL = 1 - piecewise quadratic representation 

DP1D and DP2D are the normal and reduced azimuth intervals, in degrees, over which 
the induced velocity contributions are summed. Typical values are 30.0 and 10.0, 
respectively. 

Airfoil Data 

The data used in this case (if NAFD = 1) consist of lift coefficients, over a range 
of incidences and Mach numbers. 

NXL = number of Mach numbers in the table 

NZL = number of incidences in the table 

MCL = array containing the NXL Mach numbers (ascending order) 

ACL = array containing the NZL incidences (ascending order) 

CL = two-dimensional array (of dimension at least (NZL, NXL) ) containing 

the lift coefficients 

The format for data entry is as below. 

Line 1: NXL, NZL (30X, 212 format) 

Col 31 - 32 NXL 

33 - 34 NZL 

Lines 2(a), (b),. . . MCL (I), I = 1, NXL (7X, 9F7.0 format) 
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Col 

<** 

00 

MCL (1) 



15-21 

• 

MCL (2) 



• 

64 1 70 

MCL (9) 


If NXL is greater than 9, additional lines are 
NXL entries have been made. 

entered with the same format until all 

Lines 3(a), (b) . 

. . ACL (1), 

CL (1, I), I = 

1, NXL (F7.0, 9F7.0 format) 

Col 

1 - 7 

ACL (1) 



8-14 

CL (1, 1) 



15 - 21 
• 

CL (1, 2) 



• 

64-70 

CL (1, 9) 



If NXL is greater than 9, additional lines are entered as below (7X, 9F7.0 format). 


Col 8-4 

CL (1, 10) 

15 - 21 

CL (1, 11) 

64-70 

CL (1, 18) 


Lines 4, 5, 6 . . . are identical in format to line 3 and contain the data with 

ACL (2), ACL (3) . . . ACL (NZL). 

DESCRIPTION OF OUTPUT 

The presentation of output is basically the same for Programs ASYMP1 and 
ASYMP2. Each page of output is formatted to fit within letter paper size (11" x 8.5"). 
Given below is a page-by-page description of the output. 

Page 1 contains the input data and some auxiliary parameters, as well as the 
locations of the collocation points at which the normal velocity condition has been 
replaced by the zero lift condition. Page 2 contains the basic solution for the 
collocation coefficients used to express the variation of the dipole strength function 
g(r. » T^). These are 55 in number, corresponding to the coefficients in equation (9) of 
reierence 2 for ASYMP1 and to the coefficients in equation (E3) of reference 3 for 
ASYMP2. Also presented on this page are the blade motion parameters (0 > a Q , a^, bj) 
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and the computed values of the rotor thrust coefficient and the moment coefficients 
about 2 the yotor X-and Y - axes. Page 3 contains the distribution of sectional lift 
/(p ft Rj ) in tabular form, at 5 radial and 24 azimuth locations. Page 4 contains a 
similar table of sectional lift*R^/th^ust per blade. Pages 5 and 6 contain tables of 
sectional pitching moment/( p ft 21 Rj ) and the center of pressure locations. Page 7 
presents the variation with azimuth of the total blade lift, the moment due to lift about 
the hub and the radial location of the center of- lift. Pages 8-15 present the 
distribution of surface pressure differential/( pu Rj ) at 5 spanwise and 10 chordwise 
locations, for every 15 degrees of azimuth. 

EXAMPLES OF JOB ENTRY, INPUT DATA AND OUTPUT 

The current versions of ASYMP1 and ASYMP2, written in Fortran IV, are intended 
to be run on the Cyber network at the Langley Research Center inasmuch as they use 
subroutine GELIM which is to be accessed from the subroutine library FTNMLIB. With 
appropriate changes, therefore, they can be run on any other system with a Fortran IV 
compiler. In order to compile and execute the programs, the basic sequence of job 
control statements would be as follows. 


GET, ASYMP1 (or ASYMP2). 

GET, AFDATA = airfoil datafile name, (if airfoil data are used) 
MAP, OFF. (if load map is not required) 

FTN, I = ASYMP1 (or ASYMP2), L = 0. 

ATTACH, FTNMLIB/UN = LIBRARY. 

LIBRARY, FTNMLIB. 

LGO. 


The control statements must be followed by the relevant input data. Examples are 
given below for two conditions: (A) data for program ASYMP1, applied to Case 1, 
y = 0.29 (see p. 16 of ref. 2) and (B) data for program ASYMP2, applied to Case 2, 
y = 0.29 using the piecewise quadratic representation. In each case the input data is 
followed by the corresponding output listing. 

Example A 

Line 1: 0.17, 5.43, 2, 0., 0.29, 6.7, 0.00394, 0. 

Line 2: 0, 1, 0, 0 

Line 3: 0. 

Line 4: 0.3, 0.5, 0.75, 0.85, 0.95 

Line 5: 15.0, 5.0, 0.1, 0 

In the first line MINF has been set to zero since airfoil data are not going to be 
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used. The second line specifies that 0 q , and are to be computed while a will be 
input. Since this example involves a teetering rotor with no coning angle, the tftird line 
specifies a = 0. 



ROOT RADIUS/ UP RADIUS- .17000 
aspect ratio- s.43coo 
NUM dEK .IP BLADES* 2 

L1NE4* r -y 1ST ( ^rjor TO TIP)- 0.00000 DEGREES 
FORWARD SPE-J/TIP SPEED- .29000 

RGTJR INCIDENCEIFOR-ARJ TILT POSITIVE)- 3.70000 DEGREES 

FRcrSTREAM MACH NUPBER- 0.00000 

THRUST COEFFICIENT- ,00394 

CONING ANGLE* 0.00000 DEGREES 

TOTAL INFLOW RATIO- .04070 

MINIMUM UT- . IGGJGI ZERO LIFT CONDITION APPLIED acLCw THIS VALUE) 
*0*fUL AZIMUTH INTERVAL* 15.00000 DEGREES 
RECUCcO AZIMUTH INTERVAL* 3.00000 DEGREES 
4IEFGIL J.ATA TABLES NOT USED 


*« 

.3 00 

PSI- 

2 2 9. C 91 DEGREES 

JT- 

.051 

ZERO 

LIFT 

CONDITION 

APPLIED 

3* ' 

. 3 00 

PSI- 

2fcI.aI3DEGREES 

n 

)— 

.013 

ZERO 

LIFT 

CONDITION 

APPLIED 

? « 

.300 

PSI* 

2 94. 5 4-j'JcSXcES 

UT- 

.036 

ZERO 

LIFT 

CONDITION 

APPLIED 
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SOLUUJN FOR CutFFICIt'ITS 


( Go ( I ) * I = i * s SP J 

20-xic- J2 -.A37<3 c- 02 -.7A2SE-u2 -.1972E-C2 -.A996E-03 
(GC(I*J>* J*1>NHM)*1*1*N3P) 


2 A2c c-uj 

-.AA19E-03 -.I377E-03 

-.1375E-G3 

• £ A loE-0 A 

2011 E-0 2 

-• o356c-02 -.ARGIE-03 

. 13A3E-02 

- » 2 S 7CE-03 

*tl99E-05 

-.1021E-02 -.3516E-03 

-.A87AE-G3 

- . 9 597E-0 A 

291;E-03 

-.7772E-C3 13AAE-03 

-.19A9E-C3 

• 3 1 87E-u3 

10UE-03 

-.B'lllE— 03 «llA2c~D3 

.3015E-03 

- . Ao£2£-0 3 


(GS(I*J)*J»l* TrM) . I»1*N3P) 


* 2 9o7t-0 3 

• 1 2 *9 E — 0 3 

.19E7E-DA 

- . 1 27 8E-G3 

-.12616-03 

2 556E-92 

•oi03c-C3 

331AE-D2 

. 3A99E-03 

.77832-03 

-•AA3it-C3 

. 116oE-G2 

— • 15A8E— 93 

-.139AE-02 

-. 2G96E-03 

. 92338-03 

3151E-C3 

-. Lol3£-03 

• 2038E-CA 

-• A07CE-0 3 

• jio i E— 0 3 

« 26A9E-C3 

-• A 271E-93 

-.2 01 3E-05 

.38776-03 


PITCH ANGLE AT 3LAOE ROOT* 7 . 693A6D6 GR EES 
CGMNG ANGLE* 0.6D0G3 DEGREES 
FLAPPING COEFFICIENT* Al* A. 27211 DEGREES 

FLAPPING COEFFICIENT* 31- .03698 DEGREES 

COMPUTED THRUST COEFFICIENT* .39A0E-02 

COMPUTED MOMENT COEFFICIENT A 3 OUT ROTOR X-AXIS»-.3092E-16 


computed moment coefficient ajout rotor y-axis* *33i3E-i7 



r a * l r- 1 


SECTIONAL UF T/ ( ?.-i J*{JMEGA**2)*(P 1**3)) 


k/9 1 1 

• 30006*03 

.53336*03 

• 75 JOt +C0 

. 35006*00 

• 9500E+00 

PS1 






C.v 

.22*96-0? 

.60226-02 

,ieie6-oi 

.20*56-01 

• 1 620E-01 

15 .C 

. i9o96-^2 

• 772*6-02 

• 1 6 3*E-0 1 

.20*76-01 

.15396-01 

30.0 

.196*6-0? 

. 66006-02 

. 1 565F-01 

. 1727E-CI 

.13196-01 

*5.0 

.1*566-02 

.55126-02 

.11096-01 

.12136-01 

.9265E-02 

ot.O 

• *5 3 2i.-w3 

. 35716-02 

.75716-02 

.6*536-02 

. 67 3 3E-02 

7 5.0 

-.76562-03 

.1*236-02 

.72*96-02 

. b*36£-02 

.71136-02 

9C.0 

-.15226-02 

• *0 76E-03 

,96506-02 

.11306-01 

•9*29E-02 

105.0 

-.12*76-u2 

.13 7 36-02 

.12576-01 

.1*256-01 

.11596-01 

12C .u 

-• 2 ?C it-C* 

. 5*736-32 

.1*516-01 

.15936-01 

.12706-01 

135 ,0 

.15626-02 

. 65316-02 

.15556-01 

.16926-01 

,13*3E-01 

liG.0 

.27316-02 

.951*6-02 

. 16636—0 1 

. 1311E-01 

. 1 * *8 E-01 

165.0 

.26636-12 

. 550*6-02 

.172*6-01 

. 191*6-01 

. 1 5**E-01 

18C.0 

.19066-0? 

.6*516-02 

• 1672E-01 

. 18976-01 

. 1 5*2E-01 

195. j 

• 59C7E-03 

. 56536-02 

.15016-01 

.17**6-01 

.1*326-01 

310.0 

- • 16 1 7E-03 

.*39*6-02 

.12896-01 

.15*86-01 

.12936-01 

32 5.0 

-.1C7CE-03 

.29976-02 

.11166-01 

.13926-01 

.11866-01 

2*C.O 

. 25e 66-03 

.19756-02 

.10056-01 

.12906-01 

.1118E-Q1 

255.0 

.19356-03 

.1867E-02 

.95286-02 

.12*06-01 

.1088E-01 

270.0 

-• 2695E-03 

.2*566-02 

.98056-02 

.127*6-01 

. 1121E-01 

265.0 

-. 39256-03 

. 31356-02 

.109*6-01 

.13976-01 

.12206-01 

3CC.0 

.38376-03 

• 3823E-32 

.12**6-01 

• 15*1E-01 

.13286-01 

315.0 

.16676-02 

.*9656-02 

.13706-01 

.16336-01 

. 1386E-01 

330.0 

. 2fcGCE-J2 

• 65*96—02 

. 1*856-01 

. 17206— Cl 

.1*266-01 

3*5.0 

.26606-02 

•7790E-O2 

.16*76-01 

. 13696-01 

• 1515E-C1 
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TAbte 2 - itcrii^Ai LiFr*si/THRUsr pep slade 


K/01 1 

PSI 

.3CGCE+00 

. 5000E + 00 

.7500E+CC 

. 88O0E+0O 

c.o 

•3635E+03 

. 129b 2+01 

. 29372 + 01 

• 3304E +C1 

15.3 

• 3162c 400 

. 1200E + 01 

.2963E+C1 

.33032+01 

3 C . 0 

.3011E+00 

.10556+01 

.2528E+Q1 

• £ 7902 +01 

A 5 . 0 

• 23r2c *00 

. 890 7E + 00 

.1792E+01 

.19502+01 

6 C • 0 

• 7322E-0X 

.593 1 E+00 

.12232+01 

. 1 3662 +01 

75.0 

-.12696+00 

.23092+00 

.11716+01 

.13712+01 

9C.0 

- • 24 3 96 +00 

.55362—01 

.15592+01 

.18252+01 

105.0 

- • 20 1 fcb +00 

.31396+00 

. 2C31E+C1 

• 2 j036 + 01 

X2C.0 

-• 4t 6 6E-02 

.8843E+00 

.2345E+01 

.23736+01 

135.0 

.25246+00 

.13676+01 

.25296+01 

.27336+01 

ISC. j 

. 441cE *'JJ 

.1537E+01 

. 2087E+C1 

.29272+01 

165.0 

• 4b 0 9t +00 

.13742+01 

.27656+01 

.30926+01 

150.0 

• 308OE+O0 

.11232+01 

.27026+01 

.30645+01 

195.0 

.95446-31 

.91342+00 

•2425E+C1 

.2819E+01 

21C.0 

— . 293eE-01 

. 71032+00 

•20826+01 

.25022+01 

225.0 

-. 172 9E-01 

. 49 4 3 E+00 

,ie04E+01 

.22502+01 

240.0 

.4179E-01 

.31916+00 

.15246+01 

.20846+01 

255.0 

• 3207E-01 

.3C17E+00 

. 1 540E + 0 1 

.20046+01 

270. 0 

-.435 5E-0X 

. 39O8E+00 

.15846+01 

• 2C58E +01 

265.0 

-. 63426—01 

. 5065E+00 

. 1769E+01 

.22576+01 

300.0 

.62G0E-01 

.61 352+00 

• 2C10E+01 

• 249QE +01 

315.0 

.2726c* 00 

. 83232+00 

. 22136+01 

.26466+01 

330.0 

, 4202L +03 

.10532+01 

.24006+01 

•2779E+C1 

345.0 

.42986+00 

.12606+01 

.25602+01 

.30206+01 


• 95 OOE ♦ 00 

• 26 1 7E + 0 1 
.2567E+01 
.. 2 1 31c +01 
.X4P7E+01 
. loeeE+oi 

• 1 1 49E +0 1 
.15236+31 
.19736+01 
.2052E+G1 
.2171E+01 
. 23 4GE + 3 1 
.24956+01 
•2491E+01 
•2314E+01 
.2089E+01 

• 191 7E + 01 
. 1907E+01 
.1758E+01 
•1811E+01 

• 1971E+01 
.2146E+01 

• 22 43E+01 
.23076+01 
.2449E+01 


1 A IL L 3 


SECTIONAL PITCHING 'nh£\T/(SHQ*{QMEGA* + 2)*( ftl»*4 ) ) 


(AaJUT CUAPT65-C 30*01 


7/7U 

. 300CE +0 0 

.50006*00 .750CE+00 

. 85C06 *00 

.95006+00 

PSI 





C • J 

• 33516-04 

• 2i92£-05 -.55386-04 

- • 1GC26-03 

- . 22 78E-03 

15. J 

.It 72c-04 

-.76306-05 - . 72 3 3E-04 

-.11556-03 

-.23756-03 

j C • 0 

. 7 2 4 c £ -9 6 

-.17306-04 -. 7620E-C4 

-.11246-03 

-.21546-03 

A 5.0 

-.70966-05 

-.33036-0* - • 6 51&E-04 

- .9294c-04 

- . 17 3 8E-03 

oO . 3 

- . 57 1 9t-05 

-.27236-04 -.43866-04 

-.57946-04 

-.13526-03 

75.0 

• 16326-05 

-. 73766-05 -.35886-04 

- . 49076-04 

- • 1 12CE-03 

7 1 . J 

. 12 8 2E— 09 

.17726-04 -.24536-04 

-.34696-04 

-.77396-04 

10 5.0 

• 26u 7c -0 4 

.31536-04 -.30366-05 

- . 1921E-04 

-.90116-04 

1 c 0 « 3 

• 3649E-C.4 

.29036-04 .12656-04 

- . 2C376-C5 

-.61186-04 

13 5.C 

.32792-09 

.20776-04 .27016-04 

.85616-05 

-. 5242E-04 

I 5 0 . W 

.lbC 3t-j9 

. 1620E-04 .25336-0 4 

• 5792E-05 

-.6087E-04 

lc 5 . 0 

-. b96 IE-05 

.12076-04 .11476-04 

-.32236-05 

-.90766-04 

loo. 0 

-. 29 3 3E-09 

.26376-03 -.59726-05 

-.25596-04 

-.1013E-03 

195.0 

-.33456-04 

-• 17 77 E-04 -.19756-04 

-.4064E-04 

-•1175E-03 

210.0 

-. 4466E-04 

-.31336-04 -.29986-04 

-.52216-04 

-.12986-03 

225.0 

-. 526CE-09 

- . 32 116-04 - . 3742E-04 

-.59776-04 

-.1381E-03 

24C.0 

-• 50d7E-04 

-.23396-04 -.3402E-04 

-.61536-04 

-.14126-03 

255.0 

-.3203E-C4 

-.14096-04 -.32296-04 

-•5747E-04 

-.14116-03 

27C.0 

- . 101 9E-U5 

-.36*76-05 -.2C75E-C4 

-.51396-04 

-.1433E-03 

265.0 

. 2796E-09 

.10*76-04 -.12386-04 

— . 4766E-04 

-.1499E-03 

30C.O 

.9311E-J9 

.25546-04 -.10336-04 

-.47466-04 

-.1580E-03 

315.0 

• 475 3E-04 

.32586-04 - , 1 348E-04 

- • 50646-0* 

-.16626-03 

330.0 

.97332-0* 

.26786-04 -.21096-04 

— • 39656-04 

-.1800E-03 

39 5.0 

.*3836-09 

•143oc-04 -.35336-04 

-.77476-04 

-.20416-03 
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table 4 


C6+Tc* Of P^sSSUnc LOCATION FkG+< LEADING 6 '3G F ( F 9 AC T X uN OF CHGPO) 


ft/Pi: 

.3CCCc+00 

.50006+00 

• 7500c + CC. 

• t 5C06 +00 

.95006+00 

PsI 






0. 3 

. 39 7 5L to 0 

.25226+00 

. 23016+00 

.21796+00 

.13726+00 

15.0 

•30666+00 

.24336+00 

.2241c+0J 

.21316+00 

.15226+00 

30. 0 

.2:252 + vu 

. 23046 + 00 

. 2 IS It + jO 

.20746+03 

•1432E+30 

45. J 

.216X2 + 00 

. 21446+00 

. 21166+00 

• 20vjlc+Cj 

.12736+00 

to.O 

• 16 7 4c +u 3 

.2J14E+00 

,207o6+'0C 

.19746 + uO 

.11 3cE + 00 

7 4.0 

•2364L+C0 

. 21626+00 

, 21766+ OC 

• c!2 26 +C0 

.14706+00 

90. J 

.1449t+JJ 

.53766+00 

.2334E+CO 

.22996+00 

.1924E+00 

104.3 

.11332+00 

. 354:6+00 

•2453F+0C 

•24126+00 

.20486+30 

1 2 0 . D 

-.77536+01 

. 26 4 7c ♦ J 0 

.23676+00 

.24926+CC 

.21356+00 

13 5.0 

. 3d 7 36 + 00 

. 26546+00 

.26136+00 

.25336+00 

.22456+00 

13 0 . 0 

. 26640+00 

.26116+00 

•26U26+CC 

•25216+00 

.22256+00 

lo 5 . 0 

.2 3510+00 

.25936+00 

,25446+00 

.24726+03 

.2153E+33 

1 d 0 . 0 

• l655t +00 

. 25 J26+00 

.24776+CC 

. 24126+00 

.20706+00 

195.0 

-.14260+00 

.22926+00 

. 24146 + 00 

.23436+00 

. 19636+00 

210.0 

•ld59E+0l 

.20336+00 

,23496+00 

.22795+03 

. 1343E+00 

22 5 .0 

.34660+01 

. 17996+00 

.22816+00 

.22196+00 

•17386+0 0 

240.0 

-.1O37E+01 

. 17086 + 00 

. 22466+00 

.21836+00 

.1674E+00 

255.0 

-.80560+00 

. 20066+00 

.2273E+00 

.21976+00 

.16516+00 

270.0 

.27470+00 

. 24036+00 

• 2 3606 +C0 

.22366+00 

.1664E+0O 

265.0 

-.20300+00 

. 27136+00 

• 2426E+C0 

.22776+00 

. lb 966 +00 

300.0 

.99510+00 

. 29366*00 

.2446E+00 

.2298E+C0 

.17226+03 

315.0 

.43436+00 

.29296+00 

.24366+00 

.22986+00 

.17176+00 

330. D 

.36916+03 

. c 7 7 3 6 +00 

.24076+00 

. 22736+CO 

•lo75E+00 

345.0 

.35760+00 

.26236+00 

.23606+00 

.22296+00 

.16196+00 



TAtLt 5 - TOTAL 3LADE lIFT. MOMENT' ABOUT HUB AMD RADIAL CENTER OF LIFT 


TOTAL 3LA06 LIFT/{SH0*( | jM6GA**2)*(Rl**4)) 

TOTAL dLAOc LIFT/THkUST R2k 8LAOE 
"G!*E>iT a IOUT HU3/(RHD*(GMcGA**2)*Ul**5) ) 

RADIAL CENTER OF LIFT/R1 

PSI TjTAL BLADc LIFT (MOMENT CENTER 





A30UT HUB 

OF LIFT 

0.0 

• 6359F-02 

.19312+01 

, 63 76E-02 

.71982+00 

X 5 . 0 

. 5657E-C2 

.13992+01 

• 6281E-02 

.72562+00 

3 L. J 

. 74242-02 

• 1 2 JOE + 01 

.53532-02 

. 72106+00 

45.0 

.54585-02 

. 83192+00 

• 3385E-C2 

.71192+00 

6 C • 

• 36422-C2 

.55356+00 

.25532-02 

. 72832+00 

7v#y 

.23 58L-02 

.93135+00 

.23002-02 

. 30492+00 

9L .0 

. J3792-C2 

. 356 IE +00 

.23622-02 

.84682+00 

1 C 5 • 0 

. A 907E-02 

.77372+00 

. 33762-02 

.80632+00 

12C. 3 

. 69332-02 

. 105-E+01 

.53702-02 

.75122+00 

135 .6 

.73 aOE-C? 

.12732+01 

. 56362-02 

' .715 22+00 

ly C |0 

. 36S0E-O2 

.19032+01 

.61132-02 

• 7043E+00 

165.0 

. 3721E-02 

• 1 + 092 + 01 

. 6215E-C2 

.71262+00 

18C .0 

. 8036E-02 

.12986+01 

.58752-02 

.73115+00 

195.0 

• 6911E-C2 

.11175+01 

. 5202E-C2 

.75262+00 

210.0 

.57692-02 

.93135+00 

. 5559E-02 

.77202+00 

225.0 

. 48995-02 

. 73982+00 

. 3933E-02 

.78422+00 

240.0 

. 9362E-C2 

.70975+00 

. 3436E-02 

.78792+00 

255.0 

. 9175E-02 

.6756E+00 

. 3291E-02 

.78822+00 

27C.0 

. 43592-02 

• 7 095c +00 

. 3437E-02 

• 7634E+00 

285.0 

. 9925E-C2 

.79596+00 

• 3359E-02 

.. .78126 + 00 

300 .0 

• 5752E-C2 

.92952+00 

. 4373E-02 

• 7c 0 2 2 +00 

315.0 

■ . C352E-C2 

. 1075E+01 

,43752-02 

.73282+00 

33C.0 

.75272-02 

.12162+01 

.53702-02 

.71 i5t+00 

345.0 

• 8 339E-02 

.13575+01 

.59272-02 

•71122+JO 
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TAoLE o -SURFACE PRESSURE 0 1 FF £R ENTI AL / < RHO* { GM t G A + + 2 ) * ( R 1 **2 )) 


AZIPJtri ANGLE* O.G DEGREES 


P/liU 

. 30CCS+0D 

. 5003E+00 

.75006+00 

.85006+00 

.95006+00 

x/C 






. 03000 

.2681E-01 

. 15 3 36 + 00 

. 35 156 + 00 

,51166+CO 

.38686+00 

. lGCCO 

.20886-01 

. 9950 £-01 

.23366+00 

.2 7726+CO 

.2519E+00 

. 20000 

•1698E-01 

• 6o 7 3E-01 

.13576+00 

•17&?£+CG 

.15-y3E. + 00 

• 3 w v w- 0 

.15 1 9E— 01 

.51216-01 

.11586+00 

.12916+00 

.99996-01 

• ^ 0 \j 3 Cj 

.13906-01 

.51206-01 

.39256-01 

.98326-01 

.08876-01 

.50030 

.U716-01 

. 33 73 ;-01 

.70556-Cu 

.76556-01 

.56356-01 

.70000 

• lOOBE-ul 

.22096-01 

.52966-01 

• 55366-ul 

.17756-01 

.90000 

.61C26-02 

.11266-01 

. 20056-01 

.19226-01 

.15036-02 

# Q *; 0 C J 

.43986-J2 

.77556-02 

.13556-01 

.12536 -Cl 

-.52205-03 

.99000 

.2016E-02 

. 3516E-02 

.57556-02 

.51S3E-C2 

- . 9029 E-0 3 

AZIMUTH ANGLE 

* 15.0 0E3RE 

ES 




k/Ri : 

. 3C0CE +03 

. 50036+00 

.750C6+00 

.35C06+C0 

• 9 E + 00 

X/C 






.05000 

. 88336-01 

• 1366 = 00 

• 3ol 36+00 

.51826+00 

. 33 39E + 00 

.10000 

.20796-01 

.93706-01 

.25526+00 

. 73076 +oO 

.25935+00 

• ZuOO'u 

.16505-01 

.ol33t-0i 

. 1 3696+00 

. 17736+00 

. 15536 + 00 

• 00000 

.i23ct-ul 

.56726-01 

.11556+00 

.12636+00 

• 9 75 2E-0 1 

.50000 

.Ulit-Ol 

.369dE-Cl 

.38656-01 

. 97606-01 

.66336-01 

■ .30000 

.97396-02 

.297 3E-01 

• 6950E-G1 

.75976-01 

.5550E-01 

• 7 0 o 3 u 

. 78 2Ct-u2 

.13956-01 

.51356-01 

• 52D6E-C1 

.15726-01 

.90000 

.51666-02 

. 95376-02 

. Io71£-01 

• 1 7obc — C 1 

.32156-03 

.95000 

• 29 8 <J t*“02 

. 55876-02 

.12526-01 

.11506-01 

-.13556-02 

.99003 

.13606-02 

.28536-02 

. 5257E-02 

.56996-02 

-.1259E-02 

AZIMUTH AMGLt 

* 30.0 0EGR2 

ES 




R/Rl * 

.30006+00 

.50006+00 

.75006+00 

.65006+00 

.95006+00 

X/C 






. C5GC0 

• 3292E-01 

. 1265S+00 

.31396+00 

•35d5E+00 

.32576+00 

.IOC 00 

.2281E-01 

• 8537E-01 

.21126+00 

.2397E+00 

.21056+00 

.20000 

.15326-01 

. 55556-01 

.1355E+00 

.15066+00 

.12266+00 

.30000 

• 11 756-01 

. 5106E-01 

• 9792E-01 

.10816+00 

.80196-01 

.50000 

.95066-02 

. 3179E-01 

.75636-01 

. 8 1C&E-C1 

• 53 52E-01 

.50000 

• 7661E-02 

• 2503E-01 

.57736-01 

.01566-01 

.35816-01 

. 70CC0 

. 598 7E— 02 

• 1517E-01 

. 33556-01 

.33966-01 

. 1068E-01 

.90000 

, 2565E-02 

.7177E-J2 

.15636-01 

.13666-01 

- . 1 5 52E-02 

.95000 

. 1766E-02 

.53706-02 

• 9fc 12E-02 

.36636-02 

-.23856-02 

.99000 

.7 998E-G3 

.21226-02 

. 5027E-02 

.35826-C2 

- . 1 5 30E-02 


AZIMUTH ANGLE* 55.0 OLG^iES 


5/kls 

.30006+00 

. 5 30 06 + 00 

. 75Q0E+00 

.85G0E+C0 

. 9500E+00 

X/C 






.G5CC0 

.28516-31 

.1113=500 

.22556+00 

.23676+00 

•2376E+00 

•iOwOO 

. 1 92 6 £-01 

. 75716-01 

. 15176 +C0 

• 17C86+G0 

.15256+00 

• ^ u 2 %j\j 

.12216-01 

. 57256-01 

, 9565E-01 

.1O62E+00 

. 8719E-01 

. 3CC00 

. 3d'J5t-02 

. 3510E-C1 

• 68836-01 

•7526E-01 

• 5569E-01 

.50v0t' 

.56276-02 

.25756-01 

.5x796-01 

,55656-CL 

• 3 591E-01 

.50000 

. 50556-02 

. 19 7 0 5-01 

. J950E-01 

. 5153E-C1 

.22126-01 

• 7 0 0 C 0 

. 2659t-02 

.11156-01 

.22156-01 

.22006-01 

.57096-02 

• *C G \j o 

. 12786-02 

.57 90E-02 

. 92956-02 

.82966-02 

-.31106-02 

• 9 "j i* \j *j 

.o5fc3E-03 

• 3 1 5 j 6-02 

.6035E-02 

.5153E-02 

-. 3237E-02 

.79000 

. 36536-31 

. 1 3 2 7 6-02 

• 2 506E-02 

• 2L21E-G2 

-.1877E-02 

AZIMUTH ANGLE 

* 60.0 DEG 5 

F.ES 




s 

• 3C0C t +00 

.53036530 

.75006+00 

• £ “OOfc ♦ GO 

.95006+00 

X/C 






. Ci 3 C C 0 

.97656-32 

.77096-01 

. 1562E+C0 

.loCOE+OO 

. 17576 + 00 

• 1 C C u 0 

.03126-02 

.51296-01 

.1G53E+C0 

. 1195E +00 

• 1 123E + 00 

• '300 0 

. 3r Oot -c 2 

.31356-01 

.05526-01 

.72536-01 

• 0357E-UI 

. 3CCG3 

• 2265E-32 

.22536-01 

. 56636-01 

.52106-01 

. 3987E-01 

• Iwvtu 

. 15176-32 

. 1 j5 3E-01 

• 3565E-01 

. 3e2o£-01 

• 2508E-01 

•50Cc0 

. o3 9 8E— 0 3 

. 12506-01 

• 2635E-01 

.26375-01 

.15816-01 

• 7 w V/ C 0 

.16726-03 

. 65oU5-02 

.15566-01 

.15736-01 

• 2082E-02 

• 9 »j C 0 o 

. 2569fc-05 

. 2530E-02 

.59166-02 

.55205-02 

-.31536-02 

. 95CGC 

. 12376-05 

. 1 59 7E-C2 

.33156-02 

• 33396-02 

-.29796-02 

.99000 

.22606-35 

.6571 £—03 

.15736-02 

. 1303E-02 

-. 1637E-02 


AZIMUTH ANGLE* 75. C Q6G*£tS 


k/Us 

• 3C00E +00 

.50006+00 

.75006+00 

. 8 5006 +C0 

.950CE+00 

X/C 






. C 3 3 0 0 

— . 15666—01 

.28256-01 

,15395+CC 

.17366+00 

•1725E+00 

. lUGGC 

-. 1075E-J1 

. 18896-01 

.98096-01 

.11636+00 

.1115E+00 

. 2C0U0 

-.7196E-02 

.11836-01 

• 6236E-01 

.7351E-C1 

• 6567E-01 

• 30o0 0 

-.55586-02 

. 85305-02 

.55356-01 

.52966-01 

.52076-01 

. 50000 

-.55696-02 

.62 77E-02 

. 35556-01 

.39976-01 

.27376-01 

.50300 

-.36826-02 

.57556-02 

. 267CE-0 1 

. 3060E-C1 

.17966-01 

.70000 

-.23516-02 

• 266 5 E-02 

. 1 557E-01 

.1728E-01 

.55156-02 

. 900 00 

-.10506-02 

. 12196-02 

.6765E-C2 

• 7270c-02 

- • 5 9036-03 

.95000 

-.65356-03 

.6392E-03 

. 5537E-02 

. 57C26-02 

-. 1C16E-02 

.99000 

- . 256 8E-03 

.37966-03 

• 1855E-02 

.19356-02 

-. 6767E-03 



AZIMUTH UGLE* 90. 0 DEGREES 


X / S I : 

. 30CCE +00 

.53006+00 

.75006+00 

.85006+00 

.950CE+00 

X/C 






• 05C UC- 

-.33936-01 

-.70376-03 

.19596+00 


. 2085E+00 

• 1 C L C 

-.22266-01 

.55136-03 

.12696+00 

.19926+00 

. 1371E + 00 

• 2 ocoo 

-, 1395E— 01 

.15936-02 

.82396-01 

.96696-01 

. 8295E-0 1 

. 30000 

-. i OCOE-Ol 

.21796-02 

.61536-01 

.71796-01 

.56916-01 

.•*0000 

-. 79996-o2 

.29 90E-C2 

.**0296-01 

.55936-01 

• 9C96E-C 1 

. aOCOO 

-.55666-02 

.28676-02 

.3«97E-01 

.99326-01 

.23826-01 

.70000 

-.30166-02 

.27116-02 

.23976-01 

• 2 7ilE-ol 

• 1339£-01 

.90000 

-.10276-C2 

.20396-02 

. 1123E-C1 

.12616-01 

. 3937E-02 

• 9'ioCO 

-. 575 7E-03 

. 16136-02 

. 755OE-02 

• 6909E-02 

.2207E-02 

• 4 1 3 0 \j 0 

-.16976-03 

.63016-03 

.32036-02 

.35526-02 

.73926-03 

AZIMUTH ANGLE 

*105.0 06 3 R 

66 S 




K/U s 

. 30C0E +00 

•53006+00 

.75006+00 

.65026+00 

.95006+00 

X/C 






• 050uO 

-. 3361E-01 

. 22536-01 

.23296+00 

. 26926+00 

• 2929E+00 

. 10000 

-, 2 ifc 6E-G1 

.17936-01 

.18996+00 

.1335E+C0 

.16176+00 

.20000 

— . 12 9 * 6- J l 

.191 3 E-01 

. 10626 + L-0 

•1210E+C0 

.10116+00 

.30000 

-.60296-02 

.12895-01 

. 8C836-01 

.91526-01 

.72136-01 

• 9CCC0 

-. 52 95E-02 

. 11626-01 

.69626-01 

•7267c -Ll 

. 53786-01 

. 5CC00 

-. 3267E-02 

.10726-31 

. 52586-01 

.5c73E-Cl 

.90566-01 

. 70000 

-.72196-03 

• 87596-32 

. 39056-01 

• 3797E-C1 

• 2203E-01 

.90000 

.56166-03 

• 55986-52 

. 1691E-01 

.18266-01 

.87396-02 

• *5000 

. 60 306-0 3 

.90826-02 

.11506-01 

• 123**6-01 

• 5 3 1 6E — 02 

• ^ y 0 u 0 

.36896-03 

.13936-02 

. 99516-02 

.52736-02 

. 2211E-02 

AZIMUTH ANOLE 

» 1 2 0 • 0 0 c G 8 £ c 3 




9/*l« 

. 3000E+00 

. 50036+00 

.75006+00 

' .35006+00 

. 95006+00 

X/C 






.05000 

-.15586-01 

. 57986-01 

.25976+00 

.2920E+00 

.25786+00 

. 10000 

-.65076-02 

. 62996-01 

.17976+00 

• 201Cc +00 

. 1737E+00 

.20000 

-.2659E-02 

.99106-01 

.12116+00 

.13936+00 

.1112E+U0 

.30000 

-.10976-03 

.355^6-01 

.93516-01 

.10236+00 

.31626-01 

.90000 

.15736-02 

.30016-01 

.75816-01 

• 32676-C1 

.62796-01 

.50000 

• 2662E-02 

.25796-01 

.62566-01 


• 9903E-01 

.70000 

.30866-02 

.18996-01 

.91736-01 

.99286-01 

.28316-01 

• 9C0QC 

• 310CE-02 

.10296-01 

.21986-01 

.22216-01 

.12396-01 

.95000 

.29126-02 

.72506-02 

.19786-01 

.15136-01 

•7960E-02 

.99000 

• 1168E-02 

• .32926-02 

.69506-02 

. 6529E-G2 

.32036-02 



AZIH'JTH Atv3lc*135.0 0339c=S 


*/<l : 
*/C 

. c 3 ;oc 

a 30006+00 

.50006+00 

.75006+00 

.85006+00 

.9500E+00 

• 1 4 4 fc E —0 i 

.14336+00 

•2742E+G0 

. 3u596 +00 

.26966+00 

. ioc*oo 

.12G5k-01 

. 10 356 + 00 

• 1 v06t ♦00 

.21136+00 

. 1827E+00 

• 200 GO 

.10706-CL 

.70936-01 

• 1295E+00 

.14216+00 

. 1184E+30 

. 300 CO 

.10146-01 

. 55726-01 

. 1C066+C0 

.10956+00 

.38066-01 

• + COGO 

.96olc-02 

.45376-01 

• 3236E-01 

• 86536-01 

.69766-01 

.SOCOO 

.91356-02 

.33416-01 

. 66476-01 

.72936-01 

. 5451E-01 

. 70COO 

. 77946-02 

. 26346-01 

.4 53 56-01 

•4632E-C1 

.3298E-01 

. +3 OOG 

.50696-02 

.13936-01 

.24326-01 

. 2456E -01 

• 1 A49E-01 

• 5 j C 0 0 

. 37366-02 

• 96536-02 

.16846-01 

. lbolc-01 

• 93 1 3E-02 

.97v,eo 

A Z 1 HUT H ANGLE 

. 17536-02 .42426-02 

*150.0 DEGREES 

. 7407E-02 

.72516-02 

.37176-02 

h /9 i t 

. 3000c +00 

. 50006+00 

.75006+00 

.85006+00 

. 95CCE + 00 

A/C 

• 0 3 v» 0 0 

.+2656-01 

• 1 66 7 6 + 00 

. 29226+00 

.32876+00 

.29236+00 

.IOC JC 

• JG*» 6 1 — G 1 

. 1 1 60c + 00 

.2C3CE+00 

.22696+00 

.19786+00 

• 20000 

.21736-01 

.78936-01 

.13786+00 

.15236+00 

.12796+00 

. 30000 

.17616-01 

.61496-01 

. 1C716+CO 

.11726+00 

• 94 90E-01 

a + Ju 00 

• 1 *• 9 G t —o l 

.5 (256-01 

.87306-01 

.99536-01 

•7385E-01 

• 5 'j Vrf 0 0 

.12796-01 

.41776-01 

.72456-01 

• 773iE-01 

. 53 35E-01 

• 7 0 C 0 0 

.92116-02 

.26246-01 

.4887E-01 

.51346-01 

.35016-01 

• ^2000 

• 5 1 0 8 6 -02 

.14 7 oc-Ql 

• 2 5546-01 

.25976-01 

. 15 19E-01 

. 05000 

. 3c62L- 02 

. 1>j20E-01 

• A 7 656-0 1 

. 17746-C1 

.97246-02 

•V4GG3 

AZIMJTH ANGLc 

. 16666—02 .4473 6—02 

* L o5 . 0 Cc'j^tES 

.77706-02 

.76736-02 

.33626-02 

? /su 

. 30006+00 

. 50006+00 

.75C0E+00 

.85006+00 

.95006+00 

x/C 

•ooocc 

.53966-01 

.14926+00 

. 30906+00 

.35296+00 

. 3 1 65E+C0 

• lUww’v 

.36926-01 

• 1035$ *00 

.21376+00 

.24276+00 

.21326+00 

. 20u0 J 

.24266-01 

. 70296-01 

.1+366+00 

. 16136+00 

.13646+00 

• JuCGO 

• 1 9c 16-0 1 

. 6+546-01 

. 1 108 E + 00 

1 1236c+G0 

• 9996E-01 

• +0C CO 

.14296-01 

.44366-01 

.89636-01 

.99136-01 

.7 676E-01 

. 500C0 

. 1 1 3 8c— J1 

. 36716-01 

.73766-01 

.30666-01 

, 5978E-01 

. 7JC00 

.70106-02 

. 24636-01 

.48916-01 

.52536-01 

.3475E-01 

a 40000 

.33006-02 

. 12356-01 

.25066-01 

.26106-01 

, 14 54E-01 

• 95C0C 

.2221E-02 

.69126-02 

.17246-01 

.17736-01 

• 9209E-02 

• 99C00 

.95556-03 

.39346-02 

.75346-02 

.76316-02 

.36286-02 


2 5 



4 L I i'i 'J T H 4 'J Ci L £ « 1 3 0 • 0 0 1 G 9 £ c S 


* / < 1 s 

• 300 Cc ♦ CO 

. 30d0E+00 

.7500E+C0 

. 35006+00 

. 950CE + 00 

* /C 

• C 5 G 0 0 

.43475-01 

.12506+00 

. 3C64E+C0 

.35546+00 

.3234E+00 

• 1 0 C w j 

.2c62;-01 

. do30 6-01 

.21085+00 

• Z 441 E + 0o 

• 21o7E+00 

. iOuOO 

.1727t-01 

.57736-01 

.1405E+00 

« 16 1 3E ♦ GO 

.1369E+00 

. iGo.'C 

. 1 1 7 2£ -01 

• 4414E-J1 

.1072E+OC 

•1222E+CG 

.98996-01 

• ^ 0 0 0 J 

.eitAE-G2 

• 3 ? 3 it E — 0 1 

• 65866-C1 

.47066-01 

.74856-01 

. ioj j; 

.5i04£-C2 

.23306-01 

. 65 4 4c -C 1 

.78436-01 

.57296-01 

.7(.„C0 

• 21 79t-0 2 

. id 73E-01 

,4i396-Gl 

• 4442E-C1 

.31906-01 

• V 0 L' jo 

. 23911-0 3 

.95175-02 

.22716-01 

. 2421E-01 

• 1249E-01 

.9iOCo 

• ^.6 7 Cc — L t 

. 35706-32 

. lsilE-Oi 

. I6o36-0l 

.77146-02 

AZIMUTH 4\G 

-.32372-04 
U*145.0 DiOf 

.2? ooc-oa 
EES 

.67376-02 

• 69 30 1 -02 

.29496-02 

K/U s 

. 30006+00 

• 5030c +00 

• 7 50 C E +CC 

. 85 006 +00 

.950CE+00 

x/C 

• uH Jv 

.2331t-01 

. 10356+30 

• 3 6 GO E + C0 

.33416+00 

• 3099E+ JO 

• 1 O Um w 

.14076-01 

. 73326-01 

.19196+CO 

• Zc 7s»c +00 

.20645+00 

• 2 vj o C 0 

. 6612E-02 

.47956-01 

. 12676+00 

. 14925+=0 

. 12396+CO 

.300 GO 

. 231 76 -02 

■ 3554 6-31 

.95776-31 

.11206+00 

• 4i66E-01 

.4oOC3 

.3.645C-33 

. 27 54 1-01 

.75915-01 

• 8c0bc-01 

. 68 346-01 

• p Oo C J 

-.127*6-02 

.21656-31 

.61186-01 

. 7C406-C1 

. 51 3CE-01 

. 7CC0g 

-.30096-0? 

.12996-01 

.38336-01 

.43776-01 

. 2704E-01 

• ^OoO'J 

-. 2664t-G2 

. 30126-02 

• 1 b 3 7t-0 L 

.2Co3t-0l 

. 94 76E-0 2 

• 9 j o 0 j 

-.2C20t-c2 

.43336-02 

. 1 2 336-01 

• 13306-01 

.55266-02 

.99000 
A Z IMUTi AN 3 

- . 9*. 2 36-03 
L6*210.0 063* 

• 1 74 36-92 
5 6 S 

.55=36-02 

. 66+9E-C2 

. 1444E-02 


.30006+30 

.5003E+30 

.75006+00 

.86006+00 

. 95006+00 

x/c 

• 0 5 G C 0 

. 1 260E-0 1 

.90746-01 

.24506+00 

.30266+00 

.28926+00 

•lOOcO 

.C164E-02 

.50336-01 

.1672E+00 

•20556+C0 

.19146+00 

. 2 jOCG 

.50316-03 

.33246-01 

.10446+00 

• 1 3 3 J c +00 

.11796+00 

. 3CG JO 

-.23636-02 

.27306-01 

.31336-01 

.49015-01 

.32716-01 

.40000 

-.4?59E-32 

. 2025E-01 

.64156-01 

.77036-01 

. 6C 36 £-01 

• 5 C G t»0 

— • 54 60E-0 2 

. 15136-01 

.51116-01 

• 60 66 E *"01 

.44286-01 

• 7 G \j 0 0 

-.63216-02 

.79426-02 

.31606-01 

. 36326-01 

.21726-01 

• 9 00 GO 

-•4b39£-32 

.29936-02 

.15006-01 

• 1677E-01 

• 6369E-02 

.95000 

-.34266-02 

.13836-02 

.101 IE-01 

• 1109E-C1 

• 3315E-02 

.99000 

-. 1573E-02 

. 769 9E-03 

• 43476-02 

.46526-02 

. 949CE-C3 


A Z I MU T ri ANGL6»225.C 0 £ G 9 c 6 $ 


?./n i 

• 300 vc 3 

.50J05*00 

.75006+00 

.35006+00 

> 9500E + 00 

A /C 






> C 3 V 0 u 

.1736E-01 

.65256-01 

.21556+00 

.27696+00 

. 2721E+00 

• IOC 00 

.9027E-02 

• 93 23 £-01 

. 1 9716+00 

. 1C73E+00 

.17916+00 

• 20CC0 

.1969E-C2 

. 2o98E-01 

. 9593E-01 

. 12066 +00 

.10896+00 

• 3 CO 0 u 

-.17086-02 

. 15296-01 

. 7G73E-01 

. 33806-01 

.75156-01 

• 4:o*.o 

-.90586-02 

.1292E-01 

.59376-01 

.6393E-G1 

. 5375E-01 

• 5 C 0 0 J 

-.561 bfc-02 

.90706-02 


.53996-01 

. 33 91E-01 

• 7 JCC J 

-. oS 6 16-02 

. 3339E-02 

• 2 598E-01 

.31536-01 

• 1733E-01 

.90000 

-.52336-02 

.92166-03 

. 1 1696—01 

.13356-01 

.39756-02 

.95000 

-.39136-02 

• 96 9 56-03 

•7926E-02 

.90936-02 

• 1675E-C? 

• 9 9 C G 0 

-.13296-32 

. 131 9 £-03 

• 3375E-02 

.37996-02 

.25166-03 

AZIMUTH A 1GL 

£»290.0 DEGREES 




p/m: 

. 30COc *03 

. 50006+00 

.75006+00 

.35006+00 

.95006+00 

A/C 






. 050^0 

. 2 37 Ac-01 

. 9 3 1 1 £-01 

.19696+00 

• 25376 + 00 

.25996+00 

. 10000 

. 1 3626-0 1 

.23516-01 

.13366+00 

.17936+00 

.1705 E+00 

• CuwOd 

. 5 3 A 31-02 

• 1790E-01 

.66396-01 

.11216+00 

• lu2 76 +00 

. 3 0 0 C J 

• i 0 6 9 b — 0 2 

.11316-01 

. 63736-61 

.32216-01 

.70116-01 

• +0 000 

-.169 St -02 

• 5l99t-j2 

.*.9206-01 

.63036-01 

• 9991E-01 

.50000 

-.35196-02 

. 5-.956-02 

•335 16— Oi 

.99026-01 

. 3962E-0 1 

• 70000 

- . 537 96-32 

• 19806-02 

. 22316-01 

.23936-01 

.19616-01 

• 9C0GG 

-•95G9t-02 

.21516-03 

.10196-01 

.12236-01 

.26006-02 

.95000 

-. 3992t-02 

, 99 3 2c-09 

. o790t-02 

.79306-02 

. 77106-03 

.99000 

— • i a 3 6 l— 02 

• 529 9E-05 

• 2 6^ t>t~C Z 

• 32£9£-w2 

-. 103CE-03 


AZIMUTH 493L£*25'./.0 Oi^-:S 


k/ k 1 : 

. 30006 +00 

. 50006+00 

.75006+00 

.65006+00 

• 9500E + 00 

X/C 






•CpuCO 

,1P?1E-01 

. 37536-01 

.16976+00 

.296UE+C0 

.2599E+00 

.10000 

. ci 99o£-02 

• 23996-01 

.12566+00 

.16776+00 

.16686+00 

.20000 

.33936-02 

.16216-01 

.31696-01 

.10736+00 

• 1C03E + 00 

• 3000C 

.12106-02 

• 1 L 596-01 

.60366—01 

. 7919E-C1 

.63276-01 

. 90000 

-. 3799E-C3 

.65296-02 

.97026-01 

• 6C79E-C1 

.9796E-01 

. 50000 

-•1892E-02 

.62976-02 

.37016-01 

.97336-01 

. 3395E-01 

. 70000 . 

-.32336-02 

.30276-0 2 

.22196-01 

.27586-01 

.13896-01 

. 900 00 

-.26526-02 

.10036-02 

.10086-01 

.11906-01 

.21656-02 

.95000 

-.2195E-02 

. 62 36E-03 

.67026-02 

.77306-02 

.97796-03 

.99000 

-.10996-02 

.25906-03 

.23976-02 

« 3189E-C2 

- . 2906E-03 
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AZIMJTH ANGLE*270.0 D6‘39£SS 


R/*n 

• 3C00E + 00 

. 5000 £ +00 

.7500E+00 

• 6 b OoE + CO 

• 9500E +00 

x/c 






• C5ooC 

-.59006-02 

• 4472E-C1 

•1851E+0C 

.25176+00 

• 262 1 E +00 

.13060 

-• 355 2E-02 

.31126-01 

.12666+00 

. 1706E +00 

• 1720E+00 

.20000 

-.19096-02 

.2101 E-01 

.83256-01 

.11026+00 

.103eE+00 

.30000 

-.13016-02 

• 16056-01 

.62526-01 

.61426-01 

• 710CE-01 

.4t>0G0 

-.1055E-02 

. 12 74E-0 L 

• 49 L 36-0 1 

• o29 36-01 

. 5019E-01 

• 5CCC0 

-.971 IE-03 

• 1017E-U1 

• 3 920E-01 

.47326-01 

• 3528E-01 

. 7CCCC 

-. 9 33 oE-03 

.62276-02 

.24256-01 

. 2 9 c.36-01 

• 1496E— 01 

• 70000 

-.63516-03 

. 23695-92 

. 1 149E-01 

• 1292E-01 

.2566E-C2 

. 9?00Q 

-.4764E-03 

.19236-02 

. 77506-02 

. 34 71E-C2 

.6641E-03 

.07000 

-.2I46E-03 

.33636-03 

.33396-02 

•3524E-02 

- • 2071E-03 

AZIMUTH ANGLE 

* 2 0 5 • G 0 r G k 

EES 




H / ■» l * 

.3orot+oo 

. 50306 + 00 

.75006+00 

. 85 OOt +00 

. 950CE+00 

X/C 






• OsOC'J 

-.1952E-C1 

.519 JE-01 

. 2O25E+0C 

.27266+00 

• 2833E +00 

• iocog 

-. Ill IE-01 

. 36996-01 

.13916+00 

.13526+00 

•1S62E+O0 

. 20GuO 

-.456 JE-02 

. 261 7 E-01 

. V2316-C1 

.12036+00 

•1128E+00 

. 30000 

- . 15 5 7E-G 2 

.20975-01 

.09996-C1 

.89376-01 

•7752E-01 

. 400C0 

• 1 7 1 Oc-Q 3 

.17476-01 

.55575-01 

.69516-01 

• 5 5 15E-01 

• 5 C 0 G j 

.12106-02 

.14696-01 

• 4432E-01 

•5437c -Cl 

. 3911E-01 

.70000 

.2091E-02 

.99906-02 

.28426-01 

.33086-01 

. 1712E-01 

• 90000 

• 17 6 25-02 

. 51755-32 

.13866-31 

.14395-01 

. 3403E-02 

• 5ucO 

• 13 bdE-G2 

. 35 736-02 

.04336-02 

. 49036-02 

.11696-02 

.77000 

. O730E-0 3 

. 15736-02 

.40976-02 

.41576-02 

- . 2 3 84E-04 

AZIMUTH ANGLt 

* 300 • 0 DEG* 

6ES 




R/Rl: 

. 30UGt *00 

.50336+00 

. 75006+00 

.35006+00 

, 9500E + 00 

X/C 






.03000 

-.U16E-01 

.58536-01 

,22926+00 

. 29 89 E +00 

.3055E+00 

. lCGuO 

-.43766-02 

. 42496-01 

. 15766+00 

.20336+00 

• 2009E + 00 

. 20000 

. 1265E-02 

.31036-01 

.10436+00 

.1324E+00 

.12176+00 

.30000 

.37606-02 

. 2565E-01 

.79736-01 

.-.861E-C1 

. 8372E-01 

.40000 

.5445E-02 

. 21985-01 

.63526-01 

. 7695E-01 

. 5961E-01 

.50000 

.6214E-02 

. 190 OE-Ul 

.51436-01 

.6097E-C1 

• 42 34E-01 

. 70C00 

.627 7E-02 

.13676-01 

• 3289E-G 1 

. 3710E-C 1 

• 1876E-01 

.70000 

. 434 2E-02 

.75326-02 

. 16196-01 

•1704E-01 

. 4078E-02 

.75000 

.32166-02 

. 5290E-02 

.11036-01 

• 1131E-01 

.16296-02 

.77000 

.15106-02 

.23656-02 

.47846-02 

. 4766E -02 

.18 69E-03 
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AZIMUTH AIGU* Jii.0 063M66S 


X / U « 

• BC OCc +00 

. 50036+00 

.7500E+00 

.85006*00 

. 950CE + 00 

X/C 






• 0 5 C 0 0 

.11336-01 

. 763 7E-01 

• 253&E+CO 

.31736*00 

•3198E+00 

• lOCwC 

• 11 396—01 

.5515E-01 

• 17A2E *0U 

• 21 olc *00 

• 209AE + 00 

• 200 CO 

•1219E-01 

• A00AE-01 

.1156E+00 

.1A06E+00 

.12656*00 

.30001 

. 1 2bcfc-31 

. 3293E-01 

.6775E-01 

• 1 O^ot ♦CO 

.8561E-01 

• ^00 00 

.12716-31 

.23 126-01 

•6982E-C1 

.eibi5-oi 

,61 3AE-01 

• 500 Go 

.12396-01 

. 2 A3 16-01 

. 56A7E-01 

, 6Ao3E— Cl 

• A331E-01 

• 7 C v C C 

.106 It “01 

. 17556-01 

• 360AE-01 

• 3 93A c-Cl 

. 1891E-01 

.90003 

. 85A dfc -02 

.97 505—02 

• 1769E-01 

• 1 & 10 E -01 

• A03EE-02 

. 95003 

. 97966—02 

. 6377E-02 

.1203E-01 

• 1203E-C1 

. lfclSE-02 

. 9 9 0 C 0 

.21916-02 

. 303 35-32 

• 52076-02 

.50756-02 

. 20 55E-03 

AZIMUTH A.mOL£« 

Biu • 0 Ot G3 c 

ES 




»/* L s 

• 300 Ct ♦ 00 

. 50005*30 

. 750CE+00 

.85006*00 

. 9500E+00 

X/C 






.03000 

• 23186-01 

.10756*00 

. 2776E*C0 

.23826*00 

.3313E+0O 

• 1 v. fj C 0 

. Z2926-01 

. 76186-01 

. 1901c*00 

.22815*03 

.2171E+00 

• 2 C v Cl 

.19Bl=-01 

.53536-01 

. 1255E+0G 

. 1 A 776 *00 

• 1303E + 00 

.301,00 

• 16 A 5E—0 1 

• A2 33E-01 

.9A78E-01 

.10966+00 

.88506-01 

• <rCCC3 

,17356-ul 

• B3dJc*"0i 

.76035-01 

» 8 A97c -Cl 

.52206-01 

.50000 

•I 0 I 8 E-CI 

. 30352-01 

. 6039E-C1 

•to99S-vl 

. A3 A1E-01 

.70000 

.13116-01 

. 21226-01 

. 38ZCE-01 

. A039E-C 1 

• 1 8 21E-01 

. 9oC 00 

.792 lt-02 

. 11A75-01 

. lc59E-01 

. 1 8A06-01 

• 3305E-02 

.95003 

.56 7 AS — J2 

. 83196-02 

. Ic626-0i 

.I220E-C1 

.10A6E-02 

.99000 

,25726-02 

• 3562E-C2 

. 5 A 5 AE-02 

.51376-02 

- . 1 7 5AE-0 A 


AZIMUTH aMGL6*3A5.0 DEGK6ES 


K/9H 

• 300CE *00 

. 5 J00E + 00 

.75006*00 

.83006*00 

.95006*00 

X/C 






• 350 CO 

. j0A3£-Oi 

.13536+00 

.31236*00 

, 370A6+C0 

• 3573E+00 

.10CUO 

.2A1AE-01 

. 9AA3E-J1 

• 2131E+00 

• 250AE+C0 

.23335*00 

. 2C00C 

• 201 8E-01 

. 6A 566-01 

• 1395E+C0 

. IfelOc + C0 

,13916+00 

. 3CCo0 

. le A1E-01 

. 50 A IE-01 

, 1 OA AE +00 

.llfaAE+OO 

• 9387E-01 

• A 0 0 0 0 

. 1707E-0 1 

. A122E-01 

. 8201E-01 

.91216-01 

• 6525E-01 

. sOCCu 

. 1 1 7 7E-01 

, 3 A 2 56-01 

. 65A76-01 

. 7133E-C1 

. AA95E-01 

• 7000C 

• 1266E-01 

.23106-01 

.A07AE-01 

. A225E-C1 

. 1 791E-01 

.90000 

.76536-02 

.120oE-01 

« 195GE-01 

• 188 Ae-Cl 

.23786-02 

.95003 

. 5A97E-02 

.83556-02 

. 1318E-01 

. 12A0E-01 

.2A10E-03 

.99000 

• 2 50 AE -02 

.36355-02 

• 5o77E-02 

•5187E-C2 

- , A990E— 03 
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Example B 


Line Is 0.16, 17.2, 4, 7., 0.29, 6.1, 0.0057, 0 
Line 2: 0, 0, 0, 0 

Line 3: 9.6 

Line 4: 0.5, 0.7, 0.8, 0.9, 1.0 

Line 5: 30.0, 10.0, 0.1, 1,0 


30 


ROOT RADUS/TIP RAOI'JS* SO/R1 * .16000 

ASPECT vATIO* 17.20000 
N 0 M 3 c R u F OLAOtS* 4 

LINEAR TWIST ( R GOT TO TIP) * 7.00003 DEGREES 

FORWARD SPE c O/TIP SPEED* .2*000 

RuTOk INCIDENCE (FQR-ARO TIlT POSITIVE) * 6.100C0 DEGREES 

FREESTktAM 1 AC H NUMBER* 0.00000 

THRUST COEFFICIENT* .DO 573 

FLAPPING INERTIA COEFFICIENT- R. 60000 

TOTAL I<FLDW RATIO- • JVG7 D 

MINIMUM OT* • lCC'OoI ZcH J LIFT CONDITION APPLIED 3tL0W THIS VALUE) 

NORMAL A7IMJTH SPACING- 3C.C00'J0 DEGREES 
RF.OuCtO AZIMUTH SPACING- 1 J.DOoDO DEGREES 

PItCE-ISc QUADRATIC APPROXIMATION OF SPANkISE DIPOLE STRENGTH VARIATION 
********************* 4 «*****«*******************+***** ***************** 

AIPFDIL DATA T43LES NO I USED 


k = 

. j30 

°S I * 261. 316 DEGREES 

*JT» 

• <H3 ZtR'J 

LIFT 

CONDITION 

APPLIED 

R* 

.300 

PSI- 2V4.5450EGREES 

UT- 

• 0 66 Z c PO 

LIFT 

CONDITION 

APPLIED 


31 



SOLUTION FOR COEFFICIENTS 


<GG(I), 1*1, N$P) 

-•1030E-01 -.3620S-01 -.5546S-01 -.6553E-01 -.61856-31 
t ( GC f I » J ) jp J«1,NHM),I*1,NSPJ 


.2 8&2E-02 
• 3532E-02 
-.7379E-03 
46206-02 
-.99E4 c -02 


-.15086-02 

-.1*796-01 

-.17486-01 

19866-01 

-.1SS4E-01 


19176- )3 
-.39o7£-)3 
-• 62176-02 
-.59946-32 
-•7231 £-32 


.62736-04 
-.20306-0 2 
-• 3336E-02 
-• 47916-02 
-.53936-32 


• 1 3136-02 

• 3 57 4E-03 
-• 3 5 59E-3 2 
-.29866-02 
-• 9otfcE-03 


( C OS < I»JJ, J*l,.\Hf1),I.l,NSP> 

-.10356-31 • 1 31 5 t-C 2 -.2362E-32 -.5932E-L3 -.487GE-03 


— • 9 391fc-02 


•1593E-02 . 65 53E- J3 


.225dE-G2 -.82596-03 


-. 342EE-03 
.98406-32 
. 13396-01 


•5335c -02 
• 5338E-02 
• 6391E-02 


• 30466-03 

• 2 3 68 E-02 
•40116-02 


•3133E-C2 
.3525E-C2 
. 2664E-02 


• 2 0346-32 
• 1 809E-02 
.63772-03 


FITCH ANGLE AT BLADE ROOT* 15.36391 DEGREES 
CONING ANGLE* 5.82645 DEGREES 
FLAPPING COEFFICIENT, Al* 6.28243 DEGREES 

FLAPPING COEFFICIENT, 31* 1.79994 DEGREES 

COMPUTED THRUST COEFFICIENT* .570JE-02 

COMPUTED MOMENT COEFFICIENT A 3 OUT ROTOR X-AX IS «-. 19886-16 
COMPUTED MOMENT COEFFICIENT ABOUT ROTOR Y-AXIS* .3755E-16 



T A 3l£ 1 - iFCTIOMAL . L I F T/ ( RHO* ( JMt 5A**2 ) * ( *1**3 ) ) 


R/U i 

. 33U06+G0 

. 50006*00 

.75006+00 

. 65005 +00 

PH 





0.0 

. I369 c -G2 

. 71636-02 

• 13266—01 

.15996-01 

15.0 

. 21 95E-02 

. 67 966 — 02 

.10686-01 

.12696-01 

30.0 

.27336-02 

. 59736-02 

• 7 to5E-02 

.86335-02 

*9.0 

.26596-02 

. 53316-02 

• 6193E-02 

.66386-02 

60. o 

. 23075-v 2 

. 53J5E-C2 

.61926-02 

.62575-02 

7 5.0 

. 25 1 IE-02 

. 55996-02 

. 6 37 6 1 — C2 

• fc4 286 -02 

9 0.0 

. 28 1 4t-02 

. 5399E-02 

.63336-02 

.64636-02 

105.0 

.33836-02 

. 59296-02 

. 6850E-02 

•C&376-02 

12C.0 

. 3o02E-02 

. 62 1 3S-02 

• 38 53E-02 

• 625C6-02 

13?. 0 

.35196-02 

. 7 lOuc-02 

. 109OE-C1 

• 1G085-C1 

15 0.0 

.337 1E-G2 

.62956-02 

. 1 1 36E-01 

.11446-01 

165.0 

.31006-02 

• o96o£-02 

. 111CE-C1 

.11996-01 

lec.c 

.29336-02 

.89336-02 

. 1050E-C1 

.12036-01 

198.0 

•1399E-32 

.71076-02 

• lCOlE-Ol 

.11926-01 

210.0 

.96 66E-03 

. 56256-02 

.92996-02 

.11356-01 

225.0 

• 5559E-09 

. 55 5 9E-02 

.31046-02 

.10416-01 

2*i C. 0 

• 5068E-09 

. 39 06E-02 

.69996-02 

.96386-02 

255.0 

.51836-09 

. 3070E-02 

.65746-02 

.93856-02 

270.0 

-.7825E-09 

. 2929E-02 

.65616-02 

.92846-02 

285.0 

-.11266-03 

.2313E-02 

.63796-02 

.89746-02 

300.0 

• 1002E-03 

.30336-02 

.64066-02 

.91006-02 

315.0 

. 38 2CE-0 3 

.93366-02 

• 7803E-02 

.10816-01 

33C.D 

. 5* 9 fc E —0 3 

. 5891 £—02 

. 107CE— 01 

.139*6-01 

345.0 

.77736-03 

.686 76-02 

. 13256-C1 

.16436-01 


• 95006 + 00 

. 1593E-C1 
.12 386-01 
. 3 190E-02 
.5290E-02 
.43646-02 
. 4689E-02 

• 5 267E-02 
.5779E-02 
.65046-02 
.76705-02 
• 9094E-02 
• 1C33E-01 
• 11C1E-01 
• 11026-01 
• 1 059E-0 1 
.1015E-01 
.9987E-02 
• 9998S-C 2 
• 9878E-02 

• 97376-02 
.1032E-01 
•12266-01 
.1499E-01 
. 1 673E-01 
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IA8LE 2 - SECTIONAL L IF T*R 1/ TH RUST ? ti 3LAQE 


H / \ l s 

PSI 

0.0 

15.0 

30.0 
At. 3 
oO.O 

75.0 
9 0.5 

10 5.0 
12C.J 

135.0 
15C.J 
les.o 

13C.0 

195.0 
210. C 

225.0 

290.0 

255.0 

270.0 

205. 0 

300.0 

315.0 
33C.O 
3*. 5.0 


• 33C06 +00 

.3J96E+C3 
. 9906c +0J 
.61C5E+OJ 
. .59296+50 
,51696+00 
. 5 16 it +C0 
.6286E+00 
.75576+00 
' .80966+05 
.78606+00 
.7529E+C0 
.69256+00 
,59366+00 
.31296+00 
.10976+00 
• 1292E-01 
.11376-01 
.1158E-01 
-.17986-01 
-.25196-01 
,22396-01 
.6533c-0l 
.12286+00 
. 17366+00 


.60006+00 

.15006+01 
.15076+01 
. 13 39E + 01 
.12026+01 
.11356+01 
. 12516+01 
.13066*01 
.13236+01 
. 13356+01 
. 15366+01 
• 1353 E + 01 
.20016+01 
. 13956+01 
.15386+01 
. 12566+01 
. 10186+01 
, 8502E+0C 
. 68 5 3E + C0 
. 59256+00 
.5177E+00 
.67756+00 
. 97976+00 
. 1 3056+01 
. 1t39E+01 


. 7 5006 +00 

.29636+01 
.2386E+01 
.1712c +01 
. 1383F+C1 
•13836+01 
.19266+01 
.19156*01 
.15396+01 
. 19116+01 
•2J36E+01 
.25376+01 
.2979E+01 
.23956+^1 
.22376+01 
. 2076E+CI 
• 18106 + 01 
.1563E+01 
.196eE+01 
• 1966E + 01 
• 1925E+01 
.19316+01 
.17936+01 
. 2 39 Of + 0 L 
. 29 59 c +0 1 


• 6 5 00 c + 00 

.35726+01 
. 282*6 + 01 
.19736+01 
. 19836+01 
. 1293E+C1 
. 1936t +01 
.19996+01 
.15336+01 
.18936+01 
.22536+01 

• 255os +01 
.26796+01 

• ct 99E + 01 

• 2b696 +01 
.2535E+C1 
.23256+01 
•21536+cl 
.20966+01 
. 20796 +C1 
.20056+01 
.20336+01 
.29196+01 
.3115E+01 
.36676+01 


.95006+00 

.35376+01 
.27656+01 
•133CE+01 
.11826+01 
.97986+00 
. 10 97E + 01 
. 11 77 E +01 
.12916+01 
.1953E+01 
.17136+01 
.2O31E+0L 
. 2308 E +01 
.2959E+J1 
. 2962E +01 
.23666+01 
• 2267E +01 
•2231E+C1 
.2233E+01 
• 2206E+01 
.21756+01 
.23056+01 
.27366+01 
• 3399E+01 
. 3737E+01 



TABLE 3 - SECTIONAL PITCHING H3H5N T/ { RHQ* ( OMcGA*** ) * ( 3 1** A ) ) 


PSI 

0.0 

15.0 
3C.9 
A 5.0 

60.0 

75.0 

90. 0 
li/6 .0 
120.0 

135.0 
15 0.0 
165.3 

130.0 

195.0 

210.0 

225.0 

240.0 

255.0 

270.0 

285.0 
3C0.3 

315.0 
33 ci • 0 

345.0 


( ABOUT CUA9TSk 
: .33006+00 

.-CHOkO) 

• 60 00 6 + 00 .750CE + 00 

.85006+00 

.95005+00 

-.2121E-05 

-.35336-06 

• b£ 69b-0o 

.12846-05 

.19395-05 

-.13936-05 

. 39396-06 

. 1 387E-05 

.20436-05 

.27105-05 

- . 7 2 306—36 

.91816-06 

• 1830E-05 

.2436E-C5 

.30456-05 

-.314ig-36 

. 13256-05 

.17716-05 

.22685-05 

.27546-05 

-. 26346-06 

.62276-06 

.11266-05 

• 14fctt-C5 

.17976-05 

-.62706-06 

-.25306-06 -.56156-07 

.79756-07 

.21576-06 

-. 12266-05 

-. 14656-05 - , 1 599E-C5 

- • 168 36-07 

-.17776-05 

- • 1 876E-05 

-.27476-05 -.32356-05 

-.3; ifa£-C5 

-• 38625-05 

-.23506-05 

-.38406-05 -.46686-05 

- . 52206 -C 5 

-.57716-05 

-.24586-05 

-.45056-05 -.56436-05 

-.84026-05 

-. 7160E-05 

-.21036-05 

-.46126-05 -.60066-05 

— • 69366—05 

-.78655-05 

-.13176-05 

-.41596-05 -.57376-05 

-.67906-05 

-.78435-05 

-.247CE-06 

-. 32716-05 -.49516-05 

-.60716-05 

- . 71 91 E-05 

.87046-06 

-.21726-05 -.38636-05 

-.49906-05 

-.61175-05 

.17756-05 

-.1122E-05 -.27316-05 

- . 3804E-05 

-•4S77E-05 

.224 8E-05 

-. 34 85E-06 -.17916-05 

- . 27 53E-05 

- . 3714E-05 

.21696-05 

• 69 53E-08 -.11946-05 

-.19956-05 

-.27966-05 

. 154 8E-05 

-.74836-07 - . 9 764E-06 

-.15776-05 

-.21795-05 

.52406-06 

-.49206-06 -. 1056E-05 

-.14336-05 

-.18095-05 

-.67146-06 

-.10546-05 - • 1 266E-05 

-.14C8E-G5 

- . 1 5 50E-05 

-.17726-05 

-.15336-05 -.14086-05 

- . 132 IE-05 

- . 12 35E-05 

-.25466-08 

-.17545-05 -. 1 314E-Q5 

-.10215-05 

-. 72746-06 

-. 26528-05 

-.15978-05 -.90238-00 

— .43876-06 

. 2 54CE— 0 1 

-.26746-05 

-.10836-05 -.20736-06 

.38025-05 

• 9676E— 06 
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TABLE A 


CENTER OF PRSiS'jRE LOCATION FROM LEADING EOG 5 t F R AC T I GN OF CHORD) 


rc/R i t 

PSI 
0 #0 

15.0 

30.0 

45.0 
eC.O 
7 : . o 
90. D 

105.3 

12C.3 

135.0 

150.0 

163.0 

lec.o 

195. C 

210.0 

225.0 
2 AC , 0 

255.0 

270.0 

235.0 

300.0 

315.0 
330. C 
3 A 5 , 0 


. 33C06 +00 

• 21 6 25 +uU 
.237Ct+uO 
• 2A46E +00 
. 2A76E +C0 
.24735+00 

• 2 A A AE +00 
.24115+00 

• 23SfcE +C 0 

• 23fctt +C0 

• 23575 +30 

• 2 3 726 +GC 
.2A13c+O0 
•24795+00 
. 2627E +00 
.3275E+00 
•1078E+01 
.1123E+01 

• 36 1 5E +00 
.1-129E+C0 
• 3721E+00 

- . 112 IE +00 

• 11 3 5c +00 
.1A37E+00 
,17951+03 


.60306+00 

. 24936+00 
•25126+00 

• 2 6 3 15 + 00 
. 2539E+03 
. 252 A E + OC 
. 24905+00 
« 24496+00 
• 2405 E+Ou 
.23746+00 
.237JE+OD 
. 23 355 + 3'w' 
. 2A05E+00 
.2A21E+00 
.2 A37E+00 
. 2A595+00 
.2A3A5+00 
. 2500E+00 
. 2A95E+00 
• 2459E+00 
. 2A07E+00 

• 2396E + 00 
.24135+00 
• 2AAAE+00 

• 2 + 6 3 E + 00 


•75CCE+CG 

. 2510E+00 
.2527E+CO 

• 2 349E +00 
.2559E+OG 
.25375+00 
.24965+00 

• 2AA3fc + C(J 
. 2ACAE+00 
.23335+00 
.23905+00 
.23925*00 
. 239 AE +00 
. 2A93E + C 0 
. 2A21E+00 
•2AACE+00 
•2A55E+00 
. 2 A o£ 5 +00 
.2A70E+0C 
.2A67E+00 
.2A59E+00 
.2A555+00 
. 2A665 + 00 
.2A83E+00 
.2A975+00 


.55005+00 

. 23165 +00 
.25335+00 
.25576+00 
.25705+00 
. 23A3E +00 

• 25C3E +00 
. 2 + 4 76 + 00 
.239A5+00 
.23705+00 
.23705+00 
.23765+00 
. 2 3 8 A5 +03 
.23975+00 
.2A1A5+00 

• 2 A315 +C0 

• 24 + bE +00 
.2+536+00 
. 2A665 + 00 
. 2a686 +00 
. 24b3t +00 
.24706+00 
. 2 Ac 16 +00 
.24946+00 
. 2505E+00 


• 9500E+00 

.2525E+GO 
.2545E+00 
.23765+00 
.26075+00 
.25845+00 
. 25095+00 
.24315+00 
.23626+00 
.2316E+CO 
.2309E+00 
.23236+00 
. 23A5E +00 
• 23666 + 00 
.2336E+00 
.2A06E+C0 
.24255+00 
.2AA3E+00 
.2A55E+00 
. 2A63E +00 
.24675+00 
.2A76E+00 
.24385+00 
•250CE+00 
.2512E+00 



TABLE 5 


TOTAL BLADE LIFT. MOMENT A30UT HUB AND RADIAL CENTER OF LIFT 


TOTAL oL A 0 £ LIFT/(-tiC*(u.'iEviA**2)*(5l**'»)I 
TOTAL olVDc LIFI/THRUST °ER BLADE 
MOMENT ABOUT HUB / t RHO* ( OIL GA+* 2) * ( R 1** 5 )) 

racial center of lift/m 

PSl TCIAL BLADE LIFT MOMENT CtNTER 





ABOUT HUB 

OF LIFT 

c • c 

. 6342E-02 

. 14176401 

• 47 85 t-C2 

.75456+00 

15.0 

. 55556-02 

.12432+31 

,4034c-02 

,725Cc*O0 

3C.C 

.44642-02 

.99712400 

.30672-02 

.68716+00 

45.0 

. 36o26-02 

. 22252400 

.24391-02 

• 6625E+00 

0 0.0 

. 34: BE-02 

.77132400 

.22 892-02 

• 66 30E +00 

7 5.0 

• 35s5E-C2 

• 9uD72 40 0 

. 2335E-C2 

.6652E+00 

90.0 

. 36256-02 

.30462400 

.25122-02 

.65672+00 

105.0 

,4 1466-02 

. 9262E 400 

.27032-02 

.65 192+00 

110.0 

•4o49c— 02 

. 103 Sc 4 01 

.30 7-32-02 

. 66 lc E +00 

135.0 

. 5 3102-02 

.11362401 

. 3590E-C? 

.67512+00 

150.0 

• 5386S-U2 

. 13152401 

.40332-02 

•6653E+00 

155.0 

. 60S9E-02 

.1350E4J1 

. 42 1 3E-02 

• 69 26E + 0 0 

100. u 

. 5304E-02 

. 1 296E 401 

• 4106E-02 

.70756+00 

195.0 

. 515dE-02 

.115 IE 401 

• 37 80E-G2 

.73376+00 

210.0 

. 44016-02 

•9B32E400 

. 3 360E-02 : 

.76336+00 

225.0 

. 3794E-02 

.84762400 

• 2970E-02 

.78272+00 

24C.0 

. 3419E-02 

. 763 8E 400 

• 2703E-02 

• 7907E+00 

255.0 

. 3192E-02 

.71292400 

. 2555E-02 

.8C04E+00 

27C.0 

. 2499c -02 

.66992400 

• 2442E-02 

. 8142E+00 

2d 5 . 0 

. 2906E-02 

.64962400 

.23722-02 

.81556+00 

300.0 

.31932-02 

.71332400 

.25436-02 

.7963E+00 

315.0 

.40582-02 

.90642400 

• 3 1 59 E— 02 

.77366*00 

330.0 

. 52 85E-02 

.11302401 

.40882-02 

.77352+00 

345.0 

.624 IE-02 

.13942401 

• 4301E-02 

.76936+00 
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TABLE o 


Su+FACc PRESSURE 0 1 FF 6 RENT I AL/ ( RHO* < OMEGA** a ) * ( R 1**2 ) ) 


AZIMUTH ANGLE- 0.0 02GR52S 


S/Rl l 

.33002+00 

• 6330E+00 

. 750CE + 00 

• 6 b 00c +00 

• 95CCE+00 

X/L 






. C500G 

. 64s 3E-0 l 

. 40572+00 

.7561E+00 

. 4 1 33E +00 

• 90b9£ +00 

.IOCCO 

.49232-01 

.27 4 +E + 00 

.52C0E+00 

• 1 278c +00 

.5229E+00 

.20000 

. 3467c-01 

. 16655+00 

,34612+00 

.41742+00 

. 4136E+00 

.30000 

• 2 7c 62-01 

. 14272+00 

• 2639E +C0 

. 31202+CO 

.31462+00 

• =*0000 

• 23 4 5E-G 1 

. 1145E+00 

•2113E+00 

• <: 5 4 32 +00 

•2513E+0Q 

.50000 

. 2gC 56-ol 

.93742-01 

.1722E+C0 

.20712+00 

•2043E+C0 

. 70000 

.14312-01 

.61572-01 

•1124E+00 

• 1 3 S c +00 

•1327E+00 

.9C0C0 

. 7892E-02 

. 31 45E-C1 

. 6 706E-01 

.68292-01 

. 67C0E-C l 

.95000 

. 55 355-02 

.21662-01 

. 39242-01 

.4p9AE-01 

•4o02E-01 

.99000 

•2+612-02 

.94962-02 

.17182-01 

.2C54E-01 

. 2U13E-01 

AZIMUTH ANGLE 

• 15. C OtGR 

2 2 S 




R / S 1 s 

. 33CCE +00 

. 60 oOE+OO 

.75002+00 

.25002+00 

.95005+00 

A / V* 

• C5C GO 

. I196E +00 

. 3348E+00 

.61222+00 

,728oE+00 

.71356+00 

. IOCCO 

.9321E-01 

. 26462+00 

.42052+00 

• vOOlc+OO 

.4893E+00 

• 2C00G 

• 5oo65-01 

• 1 7 5 1 2 + 00 

.27912+00 

.23172+03 

•3239E+00 

• 3 0 0 0 v 

. 44 1 9E-0 1 

.13422+00 

•2123E+OC 

.25202+00 

• 245bE + 00 

.40000 

• 361fcE-91 

. U74E + 00 

. 16432 +C0 

. 2U10E+00 

. 1955E +00 

• 30 C 00 

. 30 1 2E-cl 

.37522-01 

. 1378E+C0 

.16322+00 

.1535E+00 

.70000 

. 204 9E-C1 

,570(32-01 

. 3943E-Q1 

• 1C 572 +00 

.1C22E+00 

•^CCOj 

. 1C83E-01 

. 23952-01 

. 45142-01 

.5324E-01 

.5128E-01 

..95000 

•7523E-02 

. 19402-01 

. 31u0E-01 

.36542-01 

.35162-01 

.99000 

. 3320E-02 

.87132-02 

. 1356E-01 

.15972-01 

. 1 5 36E-0 1 

AZIMUTH ANGLE 

- 30.0 OEGR 

2 2 S 




R/Rl : 

.330GE+00 

. 60002+00 

. 7 5 JOE + OO 

.65002+00 

• 9500E+00 

X/C 






.05000 

• 1526E+00 

• 3423E+00 

.4423E+00 

•511QE+00 

.4767E+00 

.10000 

.10552+00 

•2354E+00 

. 3033E+00 

*35012+09 

. 3262E +00 

.20000 

• 7094E-01 

.15616+00 

. 20C6E+00 

.2313E+00 

•2148E+00 

.30000 

.54652-01 

•1185E+00 

. 152CE+00 

.17512+00 

.1621E+00 

.40000 

.4420E-01 

. 94 5 4 E-01 

.1209E+C0 

.1391E+00 

.12842+00 

.50000 

.36406-01 

. 7o98E-01 

• 9 796E-01 

.11252+00 

.1035E+00 

.70000 

•2423E-01 

• 4982E-01 

• 63112-01 

. 7231E-01 

.66062-01 

.90000 

.12 54E-01 

. 25 1 IE-01 

.3161E-01 

.3613E-01 

.32772-01 

.95000 

• S6s 92— 02 

. 1723E-01 

•21o76-01 

.2474E-01 

• 224 IE-01 

.94000 

. 3B10E-C2 

.75342—02 

. 946GE-C2 

. lOdOE-Cl 

.97636-02 


AZIMUTH AMGL E* 45.0 0CGR6ES 


J/«ll 

• 33CCb + 00 

. 60005+00 

. 750C6+C0 

.65005+00 

. 95 OOE + 00 

X/C 






• 05006 

• 14966 +00 

. 30985+00 

.3585E+00 

.38566+00 

.31096+00 

• 1 C u C 0 

• 1032c +00 

. 21245+00 

• 2456E+0C 

. 26406+00 

. 21 22E+00 

. 20C00 

. 69 06E-01 

.14075+00 

.16225+00 

.17406+00 

. 1391E+00 

.30CCC 

• 529 5E-01 

. 10685+00 

.1227E+00 

.1314E+CC 

. 1 C44E+00 

• 40300 

.42626-01 

. bsllE-ol 

.9746E-01 

.10426+00 

.82296-01 

• c sj L 0 

.8493E-01 

.69055-01 

.78945-01 

.84UE-CI 

.66016-01 

. 70000 

.2305E-1U 

.44635-01 

.50625-01 

. 5 3 oOc-Cl 

• •» 1 67E-0 1 

• 4 0 G 0 

.Ue2E-01 

.22435-01 

. 2527E-01 

.26755-01 

•2043E-01 

. 9a00u 

.8153E-02 

.13395-01 

.17315-01 

. 162UE-01 

.13926-01 

• 4 00 

. 337 7E-02 

• 67 2 4E-02 

.75516-0 2 

.7977E-02 

• 60 53E-02 

AZi.UTH ANGLE 

* oO .0 0EGR6 

ES 




r m « 

• i3CCx+0J 

. 60005+00 

. 75006+00 

.85006+00 

• 95G0E+00 

" / » 
• O-iwOO 

.1300c 400 

. 3038E+00 

. 355GE+00 

.36106+00 

• 2547E+00 

• 1000k 

• 1 ^6 8l— 0 l 

. 2037E+OC 

• 2 4435 + 00 

.2475E +C0 

.17 41E + 00 

• 3 G u gj 

.600 36-01 

.13365+00 

.1 6196+CC 

.16376+00 

• 1145E+00 

. 30000 

. 4o03E-01 

. 10 555 + 00 

• 1 229E +00 

.12416+00 

.86 31E-01 

• 4 o 0 o 0 

•370oE-0i 

• 3422 E-01 

. 9797E-01 

.96746-01 

.68266-01 

.50000 

.50380-01 

.63505-01 

. 79516-01 

. 6C00E-01 

. 5497E-01 

. 70o00 

. 2005E-01 

• 445 3E-J1 

.51426-01 

. 5156E-C1 

.34966-01 

.->0000 

• 1029E-01 

.22436-01 

• 2 5866-0 i 

.25836-01 

. 173CE-01 

.95000 

. 7094E-02 

.15445-01 

.17755-01 

.17716-01 

.11826-01 

.99000 

.311 3E-02 

.67545-02 

. 7755E-02 

.77336-02 

.51476-02 

AZIMUTn A9GL6 

* 75.0 05GREE3 




s/ais 

• 33006+C0 

. 6000 c + OO 

.7500E+G0 

. 8 50CE +00 

.95006+00 

x/c 






.35000 

.12906+00 

.31725+00 

.36225+00 

.36556+00 

• 2672E + 00 

•1CCC0 

. 691 6E-01 

• 2134E+00 

• 2493E+00 

.25156+00 

.18386+00 

. 2C000 

.3997E-01 

.14595+00 

.16626+00 

.Ie76c+C0 

.12236+00 

• 30*00 

.46215-01 

. 11155+00 

• 127CE+00 

.12506+00 

.93296-01 

.40000 

. 3 73 8E-01 

.89595-01 

.10196+00 

.10266+00 

• 74 68E-01 

• 5C0U0 

. 30 7 9E-C 1 

.73265-01 

.83196-01 

.83716-01 

• 6089E-01 

. 70C00 

.2051E-01 

.48116-01 

. 5450E-01 

. 54756-Ci 

. 39 74E-01 

• "i U C J *’ 

.10625-01 

. 24575-01 

.27766-01 

.27866-01 

• 20 17E-01 

.95u00 

.73405-02 

. 1692E-01 

.19116-01 

. 1917E— 01 

. 1387E-01 

• 99C00 

• 3226E-02 

.74186-02 

•8373E-02 

.83966-02 

• 6074E-C2 
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azimuth angle* ;o.o degrees 


R/U» 

. 33 G Lb +CC 

• 60 006+00 

. 7500E+00 

. 35006 +Cu 

.9:006+00 

X/C 

< C5 C DG 

•1553E+00 

. 326:6+00 

•3539E+C0 

.36096+00 

. 2927E + 00 

. 1G0CC 

.1C77E+00 

.22576+00 

. 2446E+0C 

.24956+00 

■ 20266 + 00 

.20000 

• 72 3 4E-01 

. 1517E+00 

.1644E+C0 

.16786+00 

.13666+00 

• 30GC0 

. 564 3E-01 

.11696+00 

• 1266E+C J 

.12926+00 

.1G55E+C0 

• **^u 30 

.45696-01 

.944t>6-01 

.1024E+0C 

.10456+00 

• 3549E-01 

• 5CCC0 

.3799E-01 

• 7776E-Ci 

.34266-01 

. 6e05£-Gl 

• 70 5 — 01 

.70000 

.25566-01 

.51726-01 

. 5tCot-01 

• 5 7 2 8 c “C 1 

• 4718E-01 

. 9000C 

. 1336E-01 

.26756-01 

.29006-01 

.29646-01 

. 24 5 3E-01 

.95000 

.92566-02 

. 18436-01 

. 20046-G1 

.20496-01 

. 1697E-01 

.99000 

AZIMUTH ANGLe 

.40766-02 • 61226-02 

*105.0 DEGREES 

. 38G56-02 

.90036-^2 

.74^56- j 2 

k/91 : 

• 33C0c +C0 

• 60 00 1 ♦ 00 

.75006+00 

. 8500E+00 

.950CE+00 

X/C 

. 05C 00 

.16526+00 

. 32646+00 

•37956+00 

.37316+00 

.31 39E +00 

. 10COO 

.12676+00 

. 2264 E + 00 

•26326+C0 

. 26256+00 

.21956+00 

. 20C00 

. S 7 4 0 E -0 1 

. 15 3 3 E + 00 

. 17326+OC 

. 178CE+00 

,14906+00 

. 30000 

.67986-01 

.11396+00 

.13836+00 

.13836+00 

.11636+00 

,<.0000 

.55486-01 

• 9: 74E-01 

•1125E+00 

.11286+00 

.95316-01 

• 56*^00 

.46106-01 

.89166-01 

.93276-01 

.93:86-01 

.794 8E-01 

.70000 

.31236-01 

• 5401E-01 

.62876-01 

.63256-01 

. 54 20E-01 

.90000 

.16446-01 

. 23 2 9 £-01 

.32936-01 

. 33226-01 

.29706-01 

.9:000 

.11406-01 

. 19506-01 

• 2232c-G l 

• 23C4E-C1 

. 1 9 9 46-01 

.99000 

AZIMUTH ANGLE 

.50266-02 
*120.0 OEGR 

. 89356-02 
EES 

• 100 56-G1 

.10156-01 

•6804E-02 

R/ R 1 > 

. 3300E +00 

. 90006+00 

.75006+00 

• 65C06 + C0 

.95006+00 

x/c 

.05000 

.1959E+00 

. 33706+00 

.46366+00 

.44936+00 

. 3480E+00 

. 10CC0 

.13636+00 

. 2353E + 00 

.32556+00 

• 3 iZ 6 E ♦ 00 

. 2432E+00 

. 20C00 

.92906-01 

. 1605E+00 

. 22 ICE +0C 

.21282+00 

.16716+00 

.30000 

.72496-01 

. 1251E + CO 

.17136+00 

. 1660E+C0 

.13146+00 

.40000 

. 5935E-01 

. 1023 E + 00 

.14026+00 

. 13596 +00 

. 1084E+00 

.50000 

.49466-01 

.85156-01 

.11656+00 

.11316+00 

.9094E-01 

.70000 

. 33 6 9E-01 

• 5789E-01 

.78866-01 

.76976-01 

• 62 76E-01 

.90000 

•1762E-G1 

.30576-01 

.41486-01 

• 4C69E-01 

• 33 60E-01 

.95000 

•1230E-O1 

. 2123E-01 

.28786-01 

.28256-01 

• 2341E-01 

.99000 

.54656-02 

.93666-02 

.12696-01 

.12476-01 

. 10 35E-C1 



A2I«UTri AN GL t » 1 3s • 0 OEGriScS 


•i/Ui 

.33C0E+00 

. O303E+00 

• 750CE+00 

. fi 50CE + 00 

.95006+00 

X/C 






.05000 

•looeE+oo 

•33676+00 

.57335+00 

.55926+00 

.6092E+09 

• 1 0 0 C w 

• i 3295 +00 

• 2o TO 6 + 00 

.39826+00 

.36216+00 

. 28 62E +00 

.20003 

.90666-01 

• 18 32E + 00 

.27035+00 

• 2602E+G0 

. 19696+00 

< 30C 00 

« 7C6 66 -jI 

• 1+28E+00 

,21016+uO 

.202TE+00 

.15516+00 

• ~t C C 0 0 

, 6 eiOE-GL 

. lioTt+OG 

.17196+00 

.16606 +00 

.12816+00 

• 3 GC'GO 

. 68696-01 

•973&E-31 

• 1 626E + CC 

• 1 3u3c + 00 

. 1 J76E+00 

. 7CGOO 

.33126-01 

.65255-31 

• 9o36E— Ol 

.96116-01 

.76676-01 

. 90000 

• 1756E-01 

.35026-01 

. 506 7E-0 1 

.69766-01 

. 3995E— 01 

• /5v'O0 

.12216-01 

.29326-01 

.33156-01 

. 3655E-01 

. Z785E-01 

.99CC0 

. 53906—02 

.10736-01 

. 1550E-C1 

. 1 5266-01 

.12 3 2E-01 


AZIMUTH AMilt*l3C.C D6G66eS 


2/9 11 

• 3'iCGE+uC 

. 60006+00 

. 7 500t + 00 

. 85006 +00 

.95CCE+00 

X/C 






. C5C-CC 

.13376+00 

.65626+00 

.62316+00 

.62636+0 0 

.6675E+00 

• 1 1/ L 0 C 

.12786+00 

.31586+00 

.93276+00 

.63516+00 

.36056+00 

.2 JCCO 

• 36986-01 

. 21 63E+00 

•293oE+0C 

.29536+00 

. 23 37E +00 

.30000 

.67316 -u 1 

. I357€*0o 

. 2232E+OC 

•22C16+CO 

.18376+00 

. 6 0 C 0 0 

.55676-01 

. 136 1E+00 

. 16616 + CC 

.16816+00 

.1516E+C0 

. 50C ou 

• 6o 1 96-o 1 

.11316+00 

.15656+00 

• 13066 +00 

.12696+00 

. 70C 00 

.31616-01 

.76506-01 

.10656+00 

• 1 066 E + 00 

.87516-01 

.90000 

.16 596 — jl 

.62326-01 

.39925-01 

. 56166-01 

• 6(j795-01 

. 93C CO 

.11326-01 

.27975-01 

. 36 05E-C1 

• 3G98E-G1 

.32596-01 

. 99000 

.50866-02 

. 12336-Oi 

. 1 6 79E-0 1 

.17196-01 

. 1661E-01 


AZIMJTH ANGLE*lo5.0 DEGREES 


P/91 i 

.33006+00 

• 60006 + 00 

. 7500E+00 

. 650C6+00 

. 9500E + 00 

X/C 






• C 3000 

.17136+00 

.69366+00 

.60936+00 

» 6561 E +00 

.55796+00 

. 1 v C 0 

.11876+00 

.36266+00 

• 6 2 30 E + 00' 

.6559E+00 

. 3690E+00 

.20000 

.30266-01 

.23136+00 

.23&9E+00 

•309dE+03 

.2560E+00 

. 30C 00 

.62166-01 

• 1798E+0G 

.2229E+C0 

.26106+00 

. 2033E+00 

.50000 

.30536-01 

. 15536+00 

• 18 17E + 0Q 

.19686+00 

.17116+00 

.50000 

.61826-01 

.12126+00 

.1803E+00 

.16366+00 

. 163 IE +00 

.70000 

•2611E-C1 

.81696-01 

.10196+00 

• 11096+00 

.93C6E-01 

.90000 

. 166 96-0 1 

.62732-01 

.53566-01 

.58396-01 

•5215E-01 

.95000 

.10176-01 

. 29656-01 

.37136-01 

.60526-01 

.36286-01 

.99CGC 

. 6680E-02 

. 13066-01 

.16366-01 

•1787E-01 

. 1603E-01 


41 



A L IMUT 4 ANGLfc* 130.0 DEGREES 


«/«l» 

. 33CG6+00 

.60006+00 

• 7 50CE + 00 

. 25002+00 

.95006+00 

x/c 






• C 3 C 0 u 

.13736+00 

.4o93E+00 

.5730E+00 

• 663sJc+0j 

• 599EE+00 

• iOGOO 

< 94o6E-0l 

• 3254E+00 

.40096+00 

.46036 +00 

.41676+00 

.2COOO 

. 63 3 3E-01 

•2198E+00 

•2715E+00 

.31246+00 

. 2840E +00 

.30000 

• 45 536-01 

. 170CE+00 

.21062+00 

.2425E+G0 

. 22 1 oc + 00 

• -tO o 0 C 

.39046-01 

. 13306+00 

• 1715E+G0 

. 1976 c +00 

.131AE+C0 

•50030 

. 319oc-0 l 

.11416+00 

• 1421E+00 

• 1640E +C0 

. 1 5 12E +00 

• 7 0 0 0 c 

.21086-01 

.7550E-01 

.95806-01 

•1107E+O0 

.1O3CE+00 

.90000 

. 10806-01 

. 393 3 £-31 

•5C19E-01 

« 5811E-C1 

• 5449E— 01 

. 95:00 

. 74966-02 

.27612-01 

• 3473E-01 

.40296-01 

.3736E-01 

•99000 

.3266E-32 

. 12 15E-01 

• 1 5 32E-01 

.17766-01 

. 1671E-01 

azimuth angle 

•195.0 CESR 

EES 




R/*l* 

, 3300E +00 

.6300E+00 

.7500E+0C 

. 6500c +00 

.95C0E+00 

X/C 






. .OSOOu 

.8270E-01 

. 3957E+00 

.5546E+0C 

.65906+00 

. 50366+00 

• iDoOO 

. 56 3 6c -01 

.2733E+00 

.3842E+0C 

• 45o72 +00 

.41936+00 

• 'c. 0 C C Q 

.3633c -01 

.13446+00 

.2594E+00 

.30666+00 

.23486+00 

.30000 

• 27 5 62-01 

.14226+00 

.2O07E+OC 

.23912+00 

.22156+00 

■ • 40000 

• 21642-0+ 

.11526+00 

.16296+00 

. 19 4? E + 00 

.13082+00 

• .50000 

,1730c -01 

, 9496E-01 

.1347E+CC 

. IsOSE+OO 

.1502E+00 

. 700 jO 

.10842-01 

. 63336-01 

•90316-01 

. 1C60E +00 

•1018E+00 

.90000 

. 52 71E-02 

. 3239E-01 

.47062-01 

. 5tt3E-01 

• 53 56E-01 

• 95000 

. 3 36 5E-02 

.22742-01 

. 325C2-C1 

• 3903E— 01 

• j? 16 t “0 1 

• 9 v C 3 3 

. 1856E-02 

. 10006-01 

. 1435E-01 

.17216-01 

.16366-01 

AZIMUTH ANGLE 

•clC.O DEGREES 




K/R1« 

.33002+00 

• t JO 0 c + 00 

• 7 5GC t + 00 

.65002+00 

, 9500E +00 

X/C 






• 0 5 0 0 w 

.33256-01 

. 3L54E + 00 

. 5179E+00 

• 63 09 1 + GO 

. 5838E+00 

• IOC 00 

.21742-01 

. 21732+00 

.3582E+0C 

.4 366E +00 

.40496+00 

.20000 

.12986-01 

. 1452E + 00 

.2411E+00 

.2943S+00 

. 2741E+C0 

. 30C00 

.97562-02 

•1L24E+0C 

.18606+00 

• 2272c +00 

i 21256+00 

.40000 

•6C92E-02 

• 9063 £-01 

•1505E+00 

.1642E+C0 

. 1730E+00 

• 5CCOO 

. 42 1 6E-02 

. 7452E-01 

.1241E+00 

.15206+00 

.14336+00 

.70000 

.17 6 7E-02 

.49416-01 

• 8275E-01 

• 1016E +C0 

. 9654E-01 

.90000 

.39452-03 

• 2549E-01 

.42916-01 

.52846-01 

.50556-01 

.93000 

, 164 5E-03 

.17 5 9E-01 

.29676-01 

. 3655E-01 

. 35O3E-01 

.99000 

• 50366-04 

• 7726E-02 

.13046-01 

. 1608E-C1 

• 1 543E-01 


\Z I MOTH 4N8L2-225.0 DEGREES 


R/R IS 
X/C 
• CiC 3C 

.3300200 

. &G002+00 

.75002+00 

.85002+00 

.95002+00 

.11532-01 

. 25 7 6 £*00 

.45382+00 

.58112+00 

.56292+00 

• 1 OGU a 

.6494E-02 

• 17762+00 

.31336+00 

.40172+00 

.38982+00 

.2000 0 

.2409E-02 

.11972+00 

.21052+00 

.27022+00 

.26302+00 

• 30 COO 

• 3737E-03 

. 9G89E-01 

. 1620E +00 

.20612+00 

.20332+00 

• 4GC CO 

-.97622-03 

. 7306E-01 

.13062+00 

.16632+00 

.16502+00 

• 5C009 

-.1575E-G2 

.59932-01 

• 1 0762 + 00 

.13662+00 

. 1363E +00 

. 7 0 l' 0 0 

-.23542-02 

. 39346-01 

.71412-01 

.9228c-01 

. 91 29E-01 

. 90CC0 

-.18392-02 

.20132-01 

.36872-01 

.47772-01 

. 4754E-01 

, 4>C 00 

-.13702-02 

.13972-01 

.25472-01 

.33012-01 

.32902-01 

• 5 90 jO 

-.64122-03 

.60942-02 

.11192-01 

. 14512-01 

• 1448E-01 

AZIMUTH ANGLE 

•240.0 DtGRE 

£5 





i< / l • 


. 60002+00 

. 75002+00 

.8500E +C0 

. 9500E+00 

X/C 



• 



• w ? 0 'j 'j 

• 109 7£ — 'j i 

. 21 o 32 + 00 

• 3*326+00 

.54022+00 

. 5571E+00 

. 10000 

.01582-02 

. 14 99£ + 00 

.27142+00 

.37312+00 

• 3852E+00 

• 2C uOu 

.22536-0’ 

• 9924o-01 

• 1 6 Z 0 L +00 

.250‘te +00 

• 2 592E + 00 

• 30000 

• 305 82-03 

.75792-01 

• 1 39 8 1 +0C 

.19262+00 

• 1998E+00 

• 4 jv. CO 

-.96932-03 

. C076E-O1 

.11272+00 

.15542 +0j 

. 1616E+00 

. 50000 

-. Io526-02 

.49612-01 

.92512-01 

.12732+03 

. 13322+00 

. 70 O0C 

-.22952-02 

.32472-01 

. 6123E-01 

.84762-01 

. 3374E-01 

• 90 COG 

- • 17 9 6c-u2 

. 16532-01 

. 31522-01 

.43722-01 

.45962-01 

. 75000 

-.1335E-02 

.11336-01 

.21752-01 

. 3019C -01 

. 3+78E-01 

• 9900 J 

-.02222-03 

. 4995E-02 

.95498-02 

.13262-01 

. 13972-01 

AZIMUTH ANG 

L E» 25 5 • 0 OEGR 

EES 




R/Rl : 

. 3300E+G0 

• 6030 c +00 

.75006+00 

.85006+00 

. 9500E +00 

X/C 






.05000 

.87092-02 

.17422+00 

.36992+00 

• 52 74E +00 

. 5o00E+00 

. luC 00 

.50022-02 

.11992+00 

.2 552E+00 

.36402+00 

.38682+00 

.20000 

.20122-02 

.80022-01 

.17106+00 

. 24402+00 

.25972+00 

.30000 

. 52 7 2E-03 

.61166-01 

.13122+00 

.16742+00 

. 1998E+00 

.40000 

-.38706-03 

.49092-01 

. IG57E+0C 

. 15112 +C0 

• 1613E+00 

.50000 

-.97712-03 

• 4010E-01 

.86742-01 

.12402+00 

.1327E+00 

.70000 

-.15052-02 

. 2a30£-01 

• 5733E-C1 

.32072-01 

.88062-01 

.90000 

-.12072-02 

.1341E-01 

.29472-01 

.42242-01 

.45466-01 

.95000 

-.9067E-03 

• 92 33E-02 

• 2033E-01 

.29152-01 

• 3139E-01 

.99000 

-•4238E-93 

. 4046E-02 

• 8923E-02 

.12802-01 

. 1 3 79E-01 
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AZIMUTH ANGL£*<£7G«G DtG^tES 


A/U s 

• 33006+00 

. 6)006+00 

.75006+00 

.65006+00 

•95OCE+00 

x/c 






• Q5C 00 

-.2995E-C2 

• 136 2E +00 

.36396+00 

.52226+00 

. 5595E+30 

.10C3G 

-.20532-02 

.99036-01 

.25956+00 

.26036+00 

.38266+0 0 

.20000 

-.18166-02 

.63116-01 

. 17066 + 00 

.29196+00 

.25636+00 

. 30000 

-.17296-02 

. 93 52E-01 

.13106+00 

.13536+00 

.19736+00 

.9CGG0 

-.16601-02 

.39166-01 

• 105bE +CG 

.19936+00 

•1591E+00 

a 00000 

-.15796-02 

. 32136-01 

• 3 oo 5E-C 1 

. 1225E +00 

.13076+00 

. 700C0 

-.13276-02 

. 21396-01 

.57226-01 

.81036-01 

• 0657E— 01 

.90000 

-. 8242E-03 

.11016-01 

. 29+96-01 

. 9 lc 66-01 

. 99 6C 6-01 

.93000 

-.59296-03 

.76006-02 

.20396-01 

.26756-01 

.30786-01 

.99000 


• 3 i 3 36-02 

.89316-02 

• 1 t 626-01 

.13526-01 

AZMUTH ANGLE 

*285.0 O60R6 

6 S 




R/91S 

.33006+00 

. 60006+00 

.75006+00 

.65006+00 

.95006+00 

X/C 






• C5CG0 

-.38956-02 

. 12736+00 

.35776+00 

.50976+00 

.5A756+00 

. IOoOC 

-.50926-02 

.33616-01 

•297CE+C0 

.39326+00 

• 3773E+00 

• .2 0000 

-.32216-02 

.59976-01 

.16536+00 

.23296+00 

.25326+00 

• 3CC v-G 

-.20226-02 

.96996-01 

.127 Ac +00 

.17916+00 

.1999E+C0 

.90000 

-.12706-32 

. 37336-01 

.10286+00 

.1A+9E+C0 

.1567E+C0 

• .50000 

-.75096-03 

.31396-01 

• 89 506-01 

.11356+00 

.12866+00 

. 70000 

-.11536-03 

.21 106-01 

. 56036-01 

. 78356-01 

.350AE-01 

.90000 

.13226-03 

.11056-01 

.23396-01 

• A0306-01 

•9379E-01 

.■»5C00 

.12396-03 

. 76596-02 

,19996-jI 

,27606-01 

.30166-01 

. 9900 J 

.65796-09 

.33716-02 

.37596-02 

.12206-01 

.13256-01 

AZIMUTH ANGLE 

*300.0 G6G96 

6 S 




K/Ul 

.33006+00 

. 60006+00 

. 750CE+00 

• 65006 + C0 

. 9500E +00 

X/C 






. U 5 0 0 u 

-.90396-03 

• 1O6O&+00 

. 3533E+00 

.51226+00 

. 56186+00 

. lOGuO 

.51356-03 

. 1156E+00 

. 2A73E+0U 

.253AE+CO 

.A012E+00 

.2 DC 00 

.16566-02 

.73916-01 

.16696+00 

•2367E+00 

.2635E+00 

.30000 

.25796-02 

.60396-01 

.1280E+00 

• 1 8 16c 4*00 

.2059E+00 

.90000 

.29906-02 

.99636-01 

.10396+00 

.19636+00 

. 1657E + C0 

.50000 

.31936-02 

. 91136-01 

.8501E-01 

.12006+00 

. 13596+00 

. 700u0 

.30856-02 

.27326-01 

. 5699E-01 

.79326-01 

• 3 963E-01 

• 90C00 

.20756-02 

. 19b06— 01 

.29196-01 

.90776-01 

.95996-01 

.95000 

.15156-02 

. 1012E-01 

.20126-01 

.28126-01 

• 3171E-01 

.99000 

.69926-03 

.99626-02 

.88906-02 

.12396-01 

. 1 391E— 01 


A Z I *i U T -• A * o L z ■ 3 1 5 • 0 0 c G ft E £ S 


*/<i : 

. 330Ct +00 

. 60006+00 

.75006+00 

,e£006+00 

. 95 CCt + 00 

X/0 






. C ; C 0 0 

• 1 2 2 36-0 l 

. 29276+OC 

.93856+00 

.61026+00 

. 6939E +00 

• 1 U Vj v 

.10C5E-01 

. J.631E + jO 

•3026E+CG 

.92066+00 

.97 8CE + 00 

.200 JO 

.36726-02 

.11 36E+00 

• 20 29E + 00 

• c813c +G0 

. 31936+00 

• 3CC0 0 

• e9 3 76-02 

. 6 79 OS-01 

• 1558E+CC 

.21556+00 

•2993E+0O 

.<.0000 

. 8O9t>E-02 

.71906-01 

.12566+00 

.17336+00 

.19636+00 

• ?C'C 00 

• 76 9 3 fc^O 2 

• 5909E-01 

.IU316+00 

.19206+00 

.1606E+0C 

. 70C jO 

.o963E-J2 

• 3963E-01 

.58296-01 

.93506-01 

.10556+00 

• 9 0 \j 0 0 

. 90 166-02 

. 20636-01 

•3512E-QI 

• 97906-01 

• 5399E-01 

• 9 i> C C 0 

.26396-02 

.19326-01 

.29296-01 

.33026-01 

• 3716E-01 

• 99000 

.1 J096-02 

. 6303E-02 

.10696-01 

• A 8 e “C 1 

. 1629E-01 

AZIMUTH MOLE* 

33C.0 0F69 

EES 




S/U t 

• 3 3 C C c ♦ 0 0 

• oO 0 OE + GG 

• 7 i «j 0 E + G C 

.85006+00 

.95006+00 

x/c 







. 2 L 6 1 £ — c 1 

.32396+00 

.30956+00 

.79076+00 

. 8521E + CO 

. 10 COL 

.16016-01 

• 2253E+ JO 

• 4 lo7t + 00 

.59956+00 

•5369E+00 

.2OJ00 

.13116-01 

.15136+00 

•2785E+00 

•36336+00 

. 39096+00 

. 3060 J 

.11676-01 

• li6‘3fc*CC 

•21336+ OC 

.27786+00 

• 29866 +00 

• ^ 0 0 0 0 

.11016-01 

.99516-01 

.17156+00 

.22306+00 

• 2399E +00 

.50000 

•lC21c~Gl 

. 77056-01 

. 1 909E +00 

• 1 322E ♦■CO 

.1959E+C0 

.701-00 

.62 79E-J2 

• 51 36E-01 

.92996-01 

.11956+00 

.1279E+03 

, 9 JC 00 

.60276-02 

.23366-01 

. 9722E-01 

.60996-01 

.65135-01 

• 9 ? 0 G u 

. 3oCC6-02 

. 13 66C-C1 

• 326 1E-C 1 

.92006-01 

.99826-01 

,990oJ 

.16266-02 

•316CE-02 

• 1 930E-C 1 

.18916-01 

.19896-01 

AZIMUTH AHoL 6* 

395.0 DECS 

6 £ S 




k/Sl S 

.33006+00 

. 60006+00 

.7500E+00 

. 85006+00 

• 9 500E +00 

X/C 






• G 5 G 0 C 

. 3A21E-C1 

•3361 fe + OC 

.7519E+00 

. 9392E + 00 

.9592E+00 

. 10CCC 

.25266-01 

. 2669E+0C 

.51766+00 

. 6927E +00 

• 6561E+C0 

.20000 

. 19126-01 

.1736E+00 

. 39536+00 

.92816+00 

.93666+00 

• 3 0 \j C 0 

.163 5E-01 

• 1371E+00 

.26386+00 

.32635+00 

.3328E+00 

• iOC 00 

.19516-01 

. HOoE + OO 

.2116E+CC 

• 261 3 1 ♦ GO 

.26636+00 

• 5 0 0 G J 

.12996-01 

. 9067E-01 

. 17296+00 

.21366+00 

.21706+00 

.70000 

• 9998E-02 

. 59 97E-01 

.11336+00 

.13966+00 

.19166+00 

.90000 

.58626-02 

• 30 3 ^ c~Cl 

.57756-01 

.70986-01 

.71 SCE-Ol 

.95000 

. 91 59E-0 2 

.21286-01 

.39766-01 

.98836-01 

.99376-01 

.99000 

• 1 868E-02 

.9392E-02 

. 1792E-01 

.21396-01 

• 2 1 6 1E-Q1 
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Sample Airfoil Data 


Printed below is a sample print of typical airfoil data to illustrate the previously 
described input format. 
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0 . 

.2 

.3 

1127 

.4 

.5 

.6 

.7 

• 75 

.3 

- 39 . 

.9 

- 1.18 

1 . 

- 1.18 

- 1.18 

- 1.18 

- 1.18 

- 1.18 

- 1.18 

- 1.18 

- 1.18 

- 21 . 

- 1.18 

-.3 

- 1.18 

-.8 

-.81 

-.83 

-.85 

-.85 

-.85 

-.71 

-.63 

- 16.5 

-.64 

- 1.007 

-.64 

- 1.007 

-.944 

-.960 

-.965 

-.965 

-.965 

-.795 

-.760 

- 15 . 

-.70 

- 1.19 

-.70 

- 1.19 

- 1.09 

- 1.055 

-.99 

-.98 

-.98 

— .8 3 

-.79 

- 14 . 

-.72 

- 1.333 

-.72 

- 1.333 

- 1.22 

- 1.096 

- 1 . 

-.97 

-.97 

-.84 

-.805 

- 13 . 

-.73 

- 1.334 

-.73 

- 1.334 

- 1.28 

- 1.12 

- 1 . 

-.96 

-.96 

-.35 

-.815 

- 12 . 

-.735 

- 1.255 

-.735 

- 1.255 

- 1.26 

- 1.13 

- 1 . 

-.947 

-.94 

-.85 

-.82 

- 11 . 

-.74 

- 1.161 

-.74 

- 1.161 

- 1.19 

- 1.12 

-.994 

-.930 

-.923 

-.85 

-.81 

• 

0 

H 

1 

-.74 

- 1.055 

-.74 

- 1.055 

- 1.01 

- 1.082 

-.985 

-.910 

-.900 

-.845 

-.805 

- 8 . 

-.73 

-.344 

-.73 

-.844 

-.88 

-.907 

-.922 

-.87 

-.84 

-.82 

-.77 

- 6 . 

-.695 

-.633 

-.695 

-.633 

-.66 

-.684 

-.741 

-.77 

-.75 

-.77 

-.72 

- 4 . 

-.593 

-.46 

-.593 

-.46 

-.46 

-.4805 

-.52 

-.5657 

-.6401 

-.7244 

-.53 

- 2 . 

-.396 

-.25 

-.396 

-.25 

-.25 

-.2399 

-.27 

-.282 

-.3274 

-.3605 

-.3469 

3 . 

-.1399 

0 . 

-.1399 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

2 . 

0 . 

.23 

0 . 

.23 

.23 

.2403 

.25 

.271 

.3245 

.3417 

.3676 

4 . 

.1245 

.45 

.1245 

.45 

.45 

.4916 

.51 

.5731 

.6252 

.6368 

.5428 

6 . 

.45 

.70 

.45 

.70 

.70 

.7224 

.75 

.7967 

.7367 

.7169 

.6430 

8 • 

.593 

.89 

.593 

.89 

.69 

.9416 

.89 

.8711 

.7513 

.7071 

.70 

10 . 

.695 

1.08 

.695 

1.08 

1.08 

1.0333 

.93 

.8802 

.8600 

.845 

.805 

11 . 

.73 

1.161 

.73 

1.161 

1.19 

1.1200 

.994 

.93 

.923 

.85 

.81 

12 . 

.74 

1.255 

.74 

1.255 

1.26 

1.1305 

1.0000 

;947 

.940 

.85 

.82 

13 . 

.74 

1.334 

.74 

1.334 

1.28 

1.12 

1 . 

.96 

.96 

.85 

.815 

14 . 

.735 

1.33 

.735 

1.33 

1.22 

1.096 

1 . 

.97 

.97 

.84 

.805 

15 . 

.73 

1.19 

.73 

1.19 

1.09 

1.055 

.99 

.98 

.98 

.83 

.79 

16.5 

.72 

1.007 

.72 

1.007 

.944 

.96 

.965 

.965 

.965 

.795 

.76 

21 . 

.7 

.8 

.7 

.8 

.81 

.33 

.85 

.85 

.85 

.71 

.68 

39 . 

.64 

1.18 

.64 

1.18 

1.18 

1.18 

1.18 

1.18 

1.18 

1.18 

1.18 


1.13 

1.18 
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APPENDIX A 


LISTING OF PROGRAM ASYMP1 



oooo o n n n 


PROGRAM MAlN{lNP'JT»aurPJT,TAPE5*lNPUT»TA?£6«0UTPUT, 
IAFQATA.TAPE1-AFDATA) 

C 

C CALCULATION OF THE UNSTEADY AIRLOADS ON A HELICOPTER ROTOR BLADE IN 
C FORWARD FLIGHT. THE METHOD USES AN ACCELERATION POTENTIAL DESCRIPTION 
C OF THE FLOW FIELD AND A MATCHED ASYMPTOTIC EXPANSION TECHNIQUE TO 
C OBTAIN A SOLUTION CORRECT TO 0 ( I / ( AR *AR ) ) . 

C REF. I - TH. VAN HOLTEN. REPORT VTH-1 89, TECHNI S CHE HQGE SCHOOL* DELFT » 

C NETHERLANDS. 

C REF. 2 - G. A. PIERCE ET AL, NASA CR 165742. MAY 1981. 

C FOR IDENTIFICATION OF PROGRAM STEPS. REFER USER'S MANUAL. 

C 

REAL "U.LAM.MCL.MINF.MLOC 

DIMENSION R30< 5 ). XOUT < ID )» PBOOt 24 ), A( 59,39 >. 8(59,1). 
1IPIVOT{59),PMIN(ZO).FI1(5),FI2(5>,GO(5),GC(5,5),GS(5,5),F1(24), 
2F2(24)»F3(24),FLTl(24),FLr2(24)»FMT(24),C0L(24)»GST(24»5),WK(59), 
3DGS2(24,5),FL1(24»5)»FL2(24,5),FM(24,5),XCP(24,5), 

4P0UT( 5, 10),PDIP(24,5) 

DIMENSION CL(5C,20),MCL(2D)»ACL(50) 

COMMON/MAIN1/NL.PI 

Yl(XL»XL1,XL2,YL1» YL2)»YL1MXL-XL1)*{ YL2-YL1)/(XL2-XL1) 

DATA NSP»NHM»NAZ,NC0r,NC0FP,QMAX/5,5,ll,59,6G»0. 5/ 

DATA X OUT/. 05, .1, .2, .3, .4, .5, .7. .9, .95, .99/ 

P I *4 . * A T AN ( I , ) 

PROGRAM STEP 1. 

READ AND WRITE INPUT DATA ANO ASSOCIATED QUANTITIES. 

REA0(5, *)R0, AR.NB, TW.MU. ALR.CT.MINF 
R £ AO ( 5 » *>M»N2,N3»N4 
IFINl • EQ • 1 ) R E AD ( 5. * ) THC 
IF ( N2 .EQ. G ) R b AD ( 5, * ) G AM A 
I F ( N 2 .EQ. 1) REAO( 5.*) AO 
IF(N3 .EQ. 1 ) RE ADC 5,* ) Al 
I F ( N4 .EQ. 1)READ(5.*)31 
READ( 5,*) (RBOCI ),I»1,NSP) 

READ! 5»*)0P1D,DP2D»UTMIN»NAF0 

NAFD-0 -- AIRFOIL TABLES NOT USED. 

NAFO-1 — AIRFOIL TABLES USED. 

IFtNAFD .EQ. 0)00 TO 8 
RE AQ( 1, 1 ) NXL. NZ L 

1 FORMAT ( BOX, 2 12 ) 

RE AO ( 1, 2) (MCL(I),I*1»NXL) 

2 FQRMAT(7X,9F7.0) 

NL1-NXL/9 
N L2*NL 1 + 1 
DO 3 I-l.NZL 
DO 3 J* 1, NL2 
Jl-( 

J2« J*9 

I F ( J 1 , GT. N XL ) GO TO 3 
IF ( J2 .GT. NXL)J2»NXL 

I F ( J .EQ. 1)READ(1,4) AC L ( I ) » ( CL ( I , J3 ) , J 3- J 1 , J 2) 

4 F0RMATCF7.0.9F7.0) 

I F ( J .GT. 1)READ(1,5) (CL( I. J3)» J3-J1.J 2) 

5 F0RMAT(7X»9F7.0) 

3 CONTINUE 
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8 CONTINUE 

DP1-0P1D*PI/130. 

CP2-OP2O*PI/180. 

WRITE !o»lu) 

10 FORMAT ( 1HI ) 

WRITE!S»20)RO*AR,NB»TW,NU» ALR»NINF 

20 FORMAT ( //6X*"R0QT RAOIUS/TIP RADIUS-", F10. 5//6X, "ASPECT RATIO*", 
1F10.5//6X, "NUMBER OF BL ADE S*"» 1 2/ /6X, 

2"LINE«R TwISTIPOOT TO T I P ) -"» F 10. 5, IX, "DEGR EES" / / bX, 

3"FQRw ARO iPEEO/TIP 3PEE0»",F10.5//bX, 

A" ROTOR INCIOENC£!FORWARO TILT POS I T I VE ) ■", F 10. 5, 1 X, "DEGR EE S "/ /6X, 
5"FR6ESTR£AM MACH NUMBER-", F 10 . 5 I 
TW-TW+PI/190. 

ALR«ALR*PI/180. 

I F < N 1 .EQ. 0 ) THC *0 . 

I F ( N 2 .EQ. 01A0-0. 

I F ( N3 .EQ. OJAl-O. 

I F I N4 .El. 0)91-0. 

WRITE(6,30)CT 

30 FOR MAT( /6X, "THRUST COEFFICIENT-", F10. 5) 

IF(N1 . EQ. 1 ) WR I TE ( b> AO ) THC 

AC F05MAT( /bX, "PITCH ANGLE AT BLADE R DOT- ", F10. 5 , IX , "DEGP E E S" ) 

I F ( N2 .EQ. 0) WR ITE I 6, 50 ) GAMA 
50 FORMATI/bX, "FLAPPING INERTIA COEF F IC IE N T»" , F 1C. 5 ) 

I F ( N 2 .EQ. 1) WR I T£( b, 60) AO 
60 FGR*AT!/6X» "CONING AN GL E -" , F10. 5, IX, "DEGR E ES " ) 

IF ( N3 .EQ. 1 ) WR I T E ( 6, 7 C ) 4 1 

70 FORMAT! /6X, "FLAPPING COEFFICIENT, Al- ", F 10 . 5, IX, "DEGR EE S" ) 

I F ( NA .El. 1 ) WR ITE ( 6, 3D > 31 

80 FORMAT! /6X» "FLAPPING C JEFF IC I ENT, 31-" , F 10 . 5, 1 X , "OEGR EE S " ) 
THC-ThC+PI/180. 

AO- AG+P 1/190. 

Al»Ai*PI/160. 

31-31*PI/190. 

LAM-MU*ALR*SaRT!.5*<-MU*MU+-SQRT(MU*MU*MU*MU*CT*CT)> ) 

WKl TE ! 6, 90) LAM, UTMIN, 0P10, JP20 

90 F0RMAH/5X, "TOTAL INFLOW R A TI 0-"» F10 . 5/ /bX, "MINI MUM UT«",F10.5, 

1" 1 Z E RO LIFT CONDITION APPLIEO BELOW THIS VALUE)"//6X, 

2"NQRM AL AZIMUTH IN TER VAL » ", F 10. 5, IX, "DEGRE E S"//6X, 

3"RED'JCc D AZIMUTH I NTERV AL - ", FID . 5 , IX, "DEGREES ") 

SLCR-1. 

CL0-0. 

IF ( NAFO .EQ. 0) WRITE! 6, 91) 

91 FORMAT! /6X, "AIRFOIL DATA TABLES NOT USED") 

IF ( NAFO .EQ. 1) WRITE! 6, 92 ) 

92 FOR,' 1ATI/6X, "AIRFOIL DATA TABLES USEO") 

CALCULATE QUANTITIES NEEDED FOR TRAJECTORY SEGMENT ADJACENT TO THE 
COLLOCATION POINT. 

ETA1P«ALJG!!3.+SQRT!5.))/2.) 

CFl-CJiH! .5*ETA1P)/SINH< .5*ETA1P) 

CF2*. 5+CF1-ETA1P 
CF3».25*CF1-SINH!£TAIP) 

EXl-EXPi-ETAlP) 

CFA-. 5-EX1 

CF5--. T A1P-, 5*EXl*.2b*EXl*EX10./B. 

L» 1 
C 

C PROGRAM STEP 2. 



o o o o o o o o r> o o o r: c> o 


C START OUTER LOOP FOR COLLOCATION. THIS LOOP SETS THE CURRENT AZIMUTH 
C STATION. 

C 

00 100 J»l,NAZ 
P30>2.*J*PI/NAZ 
P 2 w U C J) «360.*J/NAZ 
CP1-C0SIP3O) 

SPI-SIN(PBO) 

C S 2»C0S(?.*P30) 

SP2«SIN(2.*PB0) 

START INNER LOOP FOR COLLOCATION. THIS LOOP SETS THE CURRENT RADIAL 
LOCATION. 

DO 100 1*1, NSP 
DO 110 M»l> NCQF 
A(L»N)«0. 

110 CONTINUE 
B ( L, 1 ) • 0. 

Z30«2.*(RBO(I)-.5*(i.tRO))/(l.-RO) 

SaZ«S3RT(l.-Z8C*ZB0) 

UT-RBO! I)+MU*SP1 

PROGRAM STEP 3. 

TEST THE TANGENTIAL VELOCITY AT THE C0LL0CAT1JN POINT, TO DECIDE 
WHETHER NORMAL VELOCITY BOUNDARY CONDITION OR ZERO LIFT CONDITION 
SHOULD BE APPLIED. 

IF( UT .GT. U T M I N ) G 0 TO 120 
WRIT£(o,l30)R30(I),PB00(J),UT 

130 FORMAT! /oX,"R*" ,FB. 3, 1X»"PSI*"» F8 . 3, "QEGRE E S", IX, "J T* ", 

IFb. 3, IX, "ZERO LIFT CONDITION APPLIED") 

GO TO l<tO 

120 DP=(1.-R0)/(2.*AR*UT) 

PLIM»PB0-(2.+R30( I ) *C P l > / S OR T < MU*MU + L AM*L AM ) 

IFINAFO .EO. 0)G0 TO 133 

PROGRAM STEP 4. 

CALCULATE LIFT CURVE SLOPE FROM DATA TABLES FOR THE CURRENT 
COLLOCATION POINT, USING THE LOCAL INCIDENCE AND MACH NJMBER. 

MLQC*UT*rtINF/MU 

ALOC*THC-TW*(RBO(I)-RG) / ( l.-RO) +B l*CPl-Al* SP1 
I-CMU*A0*CP1«-LAM)/UT 
AL0C*ALQC*160.ZPI 

CALL TABSCH(MCL,NXL,ML0C,IMCL1, IMCL2,INT) 

IFIINT .£0. -1) IMCL1»IMCL2»1 

I F ( INT .EQ. 1) rMCLl«IMCL2«NXL 

CALL TABSCHI ACL , NZL, ALOC , IACL1, I ACL2, INT) 

IF ( INT .EQ. 0 ) GO TO 131 
I F < INT .EO. — I ) GO TJ 132 
IACL1-NXL-1 
I ACL2*NXL 
GO TO 131 
132 I AC L 1* 1 
IACL2-2 

131 SLC1MCU IACL2,IMCL1)-CL( IACL1, IMCL1) ) / ( AC L ( I AC L 2 ) 

1-ACL ( I ACL 1) ) 

CL01*CL( IACL1,IMCL1)-SLC1*ACL(IACL1) 

SLC1-SLC1*180./PI 


SLC2*(CL( UCL2, IMCL2J-CU UCLl, IMCL2 )) / (AC L ( I AC L2 ) 

1-AC L ( I ACL 1 ) ) 

CL02*CL( IACLl, IMCL2 ) -SLC2* ACL (I AC LI ) 

SLC2»SLC2*130./PI 

IF(IMCL1 .EO. IMCL2)SLCR*SLC1/(2.*PI> 

IFdMCLl .NE. IMCL2)SLCR*YL{ML3C»rtCL(IMCLl),MCL(IMCL2)» 
1SLC1,SLC2)/(2.*PI) 

IFdMCLl .EO. IMCL2)CLJ*CL01 

IFdMCLl .NE. IMCL2)CL0*YL{HL0C»MCL(IMCL1)»MCL(IMCL2)» 
1CL01»CL02) 

133 CONTINUE 
C 

C CALCULATE THE CONTRIBUTION TO THE INDUCED VELOCITY FROM THE 
C TRAJECTORY SEGMENT AOJACENT TO THE COLLOCATION POINT, AND AOD TO 
C COEFFICIENT MATRIX ELEMENT. ALTERNATIVELY. SET UP THE ZEPO LIFT 
C CONDITION. 

C 

I AO DO 150 11*1, NSP 

IF(I1 .GT. 1 ) GO TO 160 
IF(uT .LE. UTMINJGO TO 170 
I F ( SLC R .EO. 0. ) GO TO 151 

xl»-AR*(l,tZB0)*CFl/UTtrtU*CPl*CF2/(UT*UT)-(2,*AR*(l.+ZB0)- 
1FUN3(R0,AR,RB0( I), 1) ) /UT 

R2*( 1 .-KO ) ♦ ( i.*Z30)*CF2/(2.*UT*UT)-(l.-R0)*MU*CPl*CF3/(2.*AR* 
1UT*UT*UT) 

W3*( l.-R J)*( l.-RO)*( l.*Z30)*CF3/( S . ♦ AR*UT* UT*UT ) 
«L«Ri-ASMl.*Z30)*d.-5LCR)/(SLCR*UT) 

GO TO 180 

151 W1*4R*( l.+ZBO) 

W2*«i3*U. 

GO TO 160 

170 fcl»AR*( 1.+Zfl0)-FUN3(R0»AR,R30(I)» 1) 

W2*0. 

W 3* 0. 

GO TO 180 
160 NL-Il-l 

I F ( UT .LE. UTMINJGO TO 19G 
P 0S*PNM ( 1 1—1, 0, Z60) 

P1S*PNM( I 1-1 , 1, Z90) 

P0*FUN3(R0, AR,RBG(I)» 3) 

I F ( SLCR .EO. 0. JGO TO 152 

»1»-AR*S0Z*P1S*CF1/(2.*PI*UT)-MU*CP1*NL*(NL+1.)*P0S+CF2/(2.*PI 
L*UT*UT) + MU*iMU*CPl*CPl*NL*CNL + l.)*PlS*CF3/(A.*PI*AR*S0Z*UT 
2*UT*UT)-(AP*SQZ*PlS»P0/2. )/ C P I * UT ) 
W2*(l.-R0)*SGZ*PlS*CF2/(A.*PI*JT*UT>+( l.-RO )*MU*CP1*NL* ( NL + 1 . ) 
i*P0S*CF3/(4.*PI*AR*UT*UT*UT) 

W3«( l.-RO )*< l.-RO) *S0Z*PIS*CF3/(16.*PI*AR*UT*UT*UT) 
«1»W1-AR*S3Z*P1S*(1.-SLCR)/(2.*PI*SLCR*UT) 

GO TO 180 

152 Wl*AR*SQZ*PlS/(2.*PI) 

V2-W3-0. 

GO TO 160 

19 0 Rl*( AR*S0Z*PNM(I1— l,l,ZBO)*FUN3(RO,AR,RBOd),3) )/(2.*PI J 
W 2 » 0 • 

W 3*0 . 

160 M*(U-1 J + NAZ + 1 

A(L»M J*A(L»H)+Wl 
UO 150 12*1, NHM 
M*M + 1 

A(L,M)»A(L»M) + (Wld2*I2*V|3)*CQS(I2*P60)-I2*W2*SIN(I2*P30) 
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M*fUl 

A(L» 1>«A(L*.M)M W1*I2*I2*W3)*SIN(I2*PB0) + I2*,«2*CQS(I2*PB0) 

150 CONTINUE 

I F { UT .or. UTMINI30 TO 203 
M«n>P*n4Z*1 

A(L»,1)«A(L»MI+( l.-Ra)*(2.*rlU*CPl>/<4.*AR)'Ml.-R0)*(l.-R0) 

1M-1. ) /(lfc.*AR*AR» 

M-MH 

A(L,i).A(L,.1)+(l.-Rj) *(2.*iU*SPl+R80( I) )/( 4.*AR) 

M*M4i 

A(L#M)»A(L#M> + (l.-80J'M-2.MU*SP2-2.*R8O{I)+CPl>/C4.*AR)t<l.-R0> 
1*(1.-R0)*(2.»SP1>/ Uo.*AR*AR» 

M-M + l 

A(L#.i)«A{L»M)>( 1«-R0 ) + ( 2.*NU*CP 2—2.*R80( I)*SP1)/(4.*AR)+(1»-R0) 
1*(1.-R0)*(-2.*CP1)/(16.*AR*AR) 

BIL.ll ■3<L#1)-TW*(2.*MU*RG*CPI-NU*MU*SP2-4.*HU*R80<I>*CP1)/(4.*AR) 
l-< l.-R0)*(-R0 + 4.**33< m-4.*MU*5Pl)*TW/ ( 16.*Aft*AR ) 

L *L ♦ 1 

oa ra loo 

200 CONTINUE 

IFISLCR .h£. 0. ) GO TO 201 
3 (L. 1 >*3<L. l)-UT*UT*CL0/{ 2.*PI) 

L*LU 
GO TO 100 

2U1 8<L»l)*3(L,i>+UT*CLO/(2.*PI*SLCRi 

C 

C PROGRAM STEP 5. 

C START CALCULATION OF THE INDUCED VELOCITY CONTRIBUTION THAT REQUIRES 
C INTEGRATION WITH AZIMUTH. 

C START LJJP FOR NUMBER OF BLADES. 

C 

DO 210 I3L*1»N8 
D3*2.*PIMI8L-1)/NB 

CALL SUBROUTINE TO DETERMINE AZIMUTH POSITIONS AT RHICH TRAJECTORY 
IS CLOSE TO A BLADE. 

I F ( DP 1 .Nt. DP2JCALL ON I N ( N U. L AM» D8, R 30 (I) » P80, P L IN* Dr. AX, 1(1 1 N , 
1P.MIN) 

Jl-1 
Kl«0 


SET LJRER AND UPPER LIMITS FOR AZIMUTH SUB-INTERVAL. 
PB2-PB0 

IFII3L .EQ. 1 ) P Bl* P 30-DP 
IFII3L .GT. 1)P81«P30-0P1 
220 CONTINUE 

Gl* • 5* ( P32+P81 ) 

G 2* . 5* ( P 3 2-P B 1 ) 

PROGRAM STEP 6. 

START LOOP FUR 5-POINT 3AUSS-CHE8YSCHEV INTEGRATION. 
00 230 U»1.5 

G«GUG2*C0S( (2.+Il-l.)*PI/10. ) 

F ACT* SORT { (G-P3 1> *( P3 2-G » ) 

00 230 M* 1» NCQF P 
Ml» ( M— 1 J/NAZ+l 
N2«N-<«i-l)*NAZ 
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IF(M1 .GT. NSP ) GO TO 240 

M3«M2/2 

M4»M2-M3*2 

I F { M 1 .GT. 1 ) GO TO 250 

I F ( M 2 .£0. UFN*FUN2tR0»AR»TW»MU>LAM»R80(I)»P80»DP»03»G»l) 

GO TO 260 
250 NL-M1-1 

IFIM2 .EO. 1)FN*FUN2(R0>AR>TW>MU*LAM.RBG(I)»P80»DP»Q3»G»2) 

2tO IF{«2 .EO. 1) A<L,M)*AU>M) + rN*FACT*PI/5. 

I F ( M2 .GT. 1 .AND. M4 .60. 0 ) A ( L» M ) « A ( L» M> +FN*COS M3* ( G + DB >> 
l*F4CT*PI/5. 

I F ( 1 2 ,C-T. 1 .AND. .M4 .GT. 0) A( Li M ) » A( l, fi ) *FN*$lN ( M 3* ( G*09 ) ) 
l*FACT*PI/5. 

GO TO 230 

2 AC I2*M+2-(NSP*NAZ) 

FN*FUS2(R0» AR, TW.1U,LAM,R30m»P30*0P» OB,G» I2J 
IF(M • LE . NCOF) AlL.MJ « A C L # M ) +FN*F ACT*P I / 5 . 

I F ( M .cJl. NC0FP)3(L»l)*8(L,l)+FN*FACT*PI/5. 

23C CONTINUE 

PROGRAM STEP 7. 

t NO OF LOOP FOR GA US S-CHE 3 YSCHE V INTEGRATION. 

OEOEFINE UPPER AND LORE R LIMITS FOR THE NEXT AZIMUTH SU3-InTERVAL . 
TEST TO SEE IF THE AZI1UTH LIMIT HAS BEEN REACHED. ALSO TEST TO SEE 
IF THE TRAJECTORY IN THE NEXT SEGMENT IS CLOSE TO A BL40E, IN WHICH 


CASE REDUCED 

SPACING 

MUST 

BE USED. 



PB2-P31 







I F t P32 . LE 

. PL I M ) GO TO 

210 




I F ( QP 1 .EO 

. DP2JG0 

TO 

270 




I F ( J I .GT. 

I MIN ) G J 

TO 

270 




I F ( K1 .EO. 

0 .AND. 

PB 2 

. LE . 

(PMIN(J1)*PI/6.))G0 

TO 

290 

I F ( K 1 .EO. 

1 .AND. 

P 32 

.GT. 

(PMINI JD-PI/6. ) )G0 

TO 

230 

IFIK1 .EO. 

1 .AND. 

P 3 2 

.LE. 

( PlINt JU-PI/6. ))G0 

TO 

290 


I F ( ?8 2 .GT. (P30-PI/6.) .AND. 1 3L .EO. 1 ) P 3 1«PB 1-DP2 
I F ( P3 2 . LE . (PBO-PI/6.) .OR. IBL .GT. 1 ) PB 1»P B1-0P1 
I F C K 1 .EO. 0 .ANO. PB1 .LT. ( PMIN ( J 1 ) + 0P l » ) PB1»P MIN ( J 1 > *DP1 
GO TO 220 
270 P81-P31-0P1 

GO TO 220 
280 P 31*P 3 1-OP 2 
< 1*1 

GO TO 220 
290 P B1 *P 81-DP1 
< 1*0 
Jl-Jl+1 
GO TO 220 
210 CONTINUE 

PROGRAM STEP 8. 

ENO OF LOOP FOR NUMBER OF 3LA0ES. 

SET UP THE COEFFICIENT MATRIX ELEMENTS CORRESPONDING TO THE 
4 AUXILIARY UNKNOWNS. AND THE RIGHT HAND SIOE. 

M »N S P +N AZ + 1 

A(L,m-A(L»M)-UT-(l.-R0)*(2.*MU*CPl)*CF4/(2.*4R*UTJ- 

1(1.-R0)*(1.-RG)*(-2.*MU*SP1)*CF5Z(4.*AR*AR*UT*UT) 

M*M*1 

A(L..M)«A(L»MJ +MU*CPI-(1.-R0)*(2.*MU*SP1+RB0(I) ) *C F4/ < 2 . * AR*UT ) 
l-( l.-RO)*(l.-RO>*( 3.*MU*C?1)*CF5/(4.*4R*AR*UT*UT) 
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m«-hi 

A(L,M>*A(L,MK.5*U*{1.-C?2DR90(I)*SP1-(1.-R0)*<-2.*MU+SP2 
l“*-**'*30tI)*CPl)*CF4/(2.*AR*UT)“{l.-R0)*(l.-R0) + (- 4.*MU*CP2 
2*2.*R8j(I)*SPl-2.*MJ*CPl*CPl)+CF5/(4.*AR*AR*UT*UT) 

1 

A(L»M)*A(L»M)-R80(I)*CPl-.5*MU*SP2”(l.-R0)*(t:.*MU*CP2”2.*R30(I) 
i*iPl)+wF4/(2.*AR*JT)-(l.-R0)*(l.-P0)*(-4.*MU*SP2-2.*H30(I) 
2*CP1-2.*MU*CP1*SP1)*CF5/(4.*AR*AR*UT*UT) 
3(L>n*J(L,l)-MU+ALH-T,M'(R3 0<n-ROI*UT/<l.-RG)-TW*{-2.*MU*RO + ''p 1 
l,rtJ*wU*SP2*A.*.MJ*S.3j( IJ*CP1)*CFA/ (2.*AR*UT)-(1.-R0)*TW* 
2(2.*MJ*R9+SP1+2.*M'J*MU*CP2-4.*'1U*R80( I ) tSPI *4 . ♦MU + MU 
3*CPX*CPi I ♦CF5/( 4.*AR*AH*UT*UT) 

L*L *1 

ICO CONTINUE 
O' 

C PROGRAM STEP 9. 

C. ENi) of collocation LOOP. 

£ CALCULATE S 3 ANWIS£ INTEGRALS NEEOEO IN SUBSEQUENT STEPS. 

OJ 300 I * 1 » N S P 

Fum«o. 

F 1 2 ( I ) *0 . 

300 CONTINUE 

Gi».5*( l.+RO) 

G2».5*< I.-RO) 

00 310 1-1,10 

G«Gi+G2*C0S( ( 2. *1-1. H-PI/20.) 

F AC T* S JR T ( (G-RC)M l.-G) ) 

•JO 310 11*1, NSP 
NL* I 1—1 

IF ( II . GT . 1 ) GO TO 320 

FIKIll-FIK I 1) +FACT*FUN3( RO, AS, G# 1 ) *? I / IQ. 

FI2(U)»FI2( II) +FACT*FJN3(R0, AR , G, 2 ) *P I / 10 . 

GO TO 310 

32 0 FIl(Il)-FIKll) ♦FACT*FUN3(SO,AS,G,3)*PI/10. 

FI2CI1) -FI2( UJ+FACT*FJN3(R0, AR,G, Al+PI/10. 

310 CONTINUE 

C PROGRAM STEP 10. 

: SET UP THE EXTRA A EQUATIONS NEEOcO TO CLOSE THE STSTEM. 

DO 330 1-1,4 
DO 340 M* 1, NCOF 
A(L,M)-0. 

340 CONTINUE 
3 ( L, 1 ) * 0. 

IF(I .EQ. I .AND. N1 ,SQ. 1JG0 TO 350 

IF(I .EQ. 2 .AND. N2 .EQ. 1)G0 TO 360 

IF(I .£0. 3 .AND. N 3 .EQ. DGO TO 370 

IF(I .EQ. 4 .AND. N4 .EQ. DGO TO 380 

I F ( I .LE. 2)11*1 
1 F ( I .GT. 2)11- 1-1 
I F ( I .GT. DGO TO 390 


C THE FOLLOWING EQUATION SETS THE TOTAL BLADE L IF T, AVERAGE 0 OVER 
C THE AZIMUTH, EQUAL TO THE THRUST COEFFICIENT. 

M-l 


A(L»M)«A(L»M1 + (1.-S0)*< 1.-R0)-(1.-RQ)*FIH1)/AR 
M-1+NAZ 
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A(L.i1)»A<L*M)'Ml.-R0)m.-R0)/<3.*PI) 

DO *00 I2-2.NSP 
M*(I2-1)*NAZ*U 

4 (L*?i)-A(L». M ) + { l»-RO) *FI1< I2)/( 2.*PI*AR) 

4GC CONTINUE 

M*NSP*N AZ+ l 

A(L*M)«A(L»M)-( l,-RO)*( l.-R0)*<l.-R0)+<l.-R0)/<16.*AR*AR*AR) 

M-M + l 

ML»1)>UliH| +< l.-RO) *<l.-RO)*< l.-RO)* <l.*ROJ/(8.*AR*AR) 

3CL*1 )*3(L>1) *CT/NB-< 1 ,-RO ) ♦ ( 1 . -RO ) * ( 1 .-RO ) ♦ ( 2. *RC) *T W/ < 16.*AR*AR* 
1AR) 

L-L+l 
GO TO 330 

THE NEXT THREE EQUATIONS (I»2*3*4) REPRESENT MOMENT EQUILIBRIUM 
A30UT THE HUB (ZEROTH HARMONIC* FIRST HARMONIC COSINE AND FIRST 
HARMONIC SINE COMPONENTS). 

39 C M»I1 

AU.#;il»A(L>N)-FI2m + AR*{ ,5*( l.-RC*RO)+( l.-RO)* (l.-RO )/£>.) 

M-M+NAZ 

A<L*M)«A(L*M)*AR*(1.-RJ*R0)/(&.*PI) 

M»M+NAZ 

A(L*M)*A(L»M)+AR*(1.-R0) *( l.-RO) /CIO. *PI) 

00 410 1 2* 2* N SP 
M»(I2-1)*NAZ*I1 
A<L,M)«A(L,M)*FI2<I2)/(2.*PI) 

410 CONTINUE 

M*NSP*N AZ+1 

I F ( I • 1 3 • 2)A(L»M»*A<L»M)-(l.-RO)*(l.-RG>*(l.-RO)*(l.«-RO)/(32.* 
1AR*AR) 

I F ( I .EQ. 3)A(L»M)«A(L*i)+MU*(l.-R0)*(l. -R 0*R0 ) /(4.*AR) 

M »M *1 

1 F { I .EO. 2)A(L*fl)*A(L*M)*:(l.-R0)*(l.-R0*P0*R0)/(12.*AR) 

1-2. /GAMA 

I F ( I .EQ. 4)A<L*M)«4{L*M)*C1.-R0)*(1.-R0*R0)+MU/<4.*AR) 

M*M + 1 

I F ( I .EQ. 3)A(L*H) m A(L*M)-(l,-RO)*(l,-RO+RO*f<OI/C6.*AR) 

I F ( I .EQ. 4) A(l#H)-4(L*D ♦ < l.-RO )*( l.-RO )+( l.-RO*RO > / (16.*AR*AR) 
M*M + 1 

I F ( I .EQ. 3)A(L»M)«4(L*H>-( 1 . -RO ) * ( 1 ,-RO ) * ( 1 ,-RO* RO > / (16 . *4R* AR ) 

I F ( I . EQ. 4)A(L*M)«A(L* 1 ) -( 1 . -RO ) * ( l.-RO*RO*RO) / < 6.* AR ) 

IF ( I .EQ. 2)B(L*l)»3(L*l)-(l«-RQ)*(l.-RO)*(l.-RO)*(l.— R0)+TW/(12. 
1*AR*4K)+R0*(1.-R0)*(1.-R'J)*(1.*R0)*TW/{32.*AR*AR) 

I F ( I .EQ. 3)B(L*1>«3( L*1)*MU*TW*{-RO*(1.-RQ*RO)/4.*(1»-RO*RO*RO)/ 
13. ) / AR 

I F ( I .EQ. 4) B ( l* 1 ) »3 ( L» 1 1 -HU*TW* ( 1 •— RO) * ( 1 •— RO*RO ) / ( 9.* AR* AR ) 

L«L + 1 
GO TO 333 

THE FOLLOWING EQUATION SETS THE COLLECTIVE PITCH EQUAL TO THE 
GIVEN VALUE. 

350 M«NSP*NAZ+i 
A(L*M)»1. 

6 (L* 1 ) »THC 
L-L + l 
GO TO 330 
C 

v. THE FOLLOWING EQUATION SETS THE CJNInG ANGLE EQUAL TO THE GIVEN 
C VALUE. 
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c 

360 rt-NSP*NAZ+2 
A(L»H)"1. 
3(L,1)»A0 
L-L *1 
GO TO 330 


c ?s E rSe“?;is G vA;t;tr iaN ssTs rHE c,cuc ,i,ch 

c 

370 H-NSP*NAZ+3 
A ( L , M > « 1 . 

3 ( L » 1 ) » A 1 
L-L + l 
GO TO 330 
C 

c r^E following equation sets the c/clic pitch coefficient, 
to the GIVEN VALUE. 

c 

36 C M»NSP*NAZ+4 
A(L»N)»1. 

3 ( L, 1 ) * 31 
L-L + l 

330 CONTINUE 


A1 » £ OU AL 


31, E3UAL 


w 

C FRJGRAN STEP 11. 

C SOLVE THE STSTEN OF SIMULTANEOUS EQUATIONS AND PRINT THE SOLUTION. 

CALL GELIN(NC0F,NC0F, A, 1,3, I P I VOT, 0, wK , I ER R ) 

IF! IE RR .20. 1)30 TO 420 
L»1 

DO 430 I1»1,NSP 
Got ii )»a ( l, i) 

L-L + l 

DO 430 12-1, NHM 
GC(I1,I2)»3(L,1) 

L-L-t-1 

GS(I1,I2)»3(L,1) 

L-L + l 

430 CONTINUE 

THC -8 (NSP*NAZ*1, l) 

THCD«THC*130./PI 

A0-8(NSP*N4Z*2»1) 

ADD -A 0*180. /PI 
A1»3(NSP*NAZ*3, 1) 

AID- Al*180./ PI 
S1-3(NSP*NAZ+4»1) 

810-31*180. /PI 

WR 1TE ( 6, 10 ) 

WRITE(6* 440 ) 

440 F0RMAT(//6X»" SOLUTION FOR COEFF IC IENTS H /6X , 25 < 14- ) / /6 X, 

1 H (G0(I),I-1,NSP)") 

WRITE (6, w50) (GO (I >, I-l, NSP) 

450 F0RMAT(/6X,5(E10.4,1X) ) 

WRITE (6,460) 

460 FORMAT! //5X» "( GC( I, J ) , J »1,NHM), I- 1,NSP )"// ) 

00 470 I-1,NSP 

WRITE! 6, 450) ( GC ( I , J ) > J ■ 1, NHM ) 

470 CONTINUE 

WRITE! 6, 480) 
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48 0 FORMAT! //6X,"(GS< I# J) , J*1,NHN). I»1,NSP )"//> 

00 <>90 1*1, NSP 

WR1TE(6,450MGS(I*J), J*i»NHM) 

49C CONTINUE 

*RITE(5*iU0)THCD» AOD, AID, 910 

500 FOR MAT(//tiX," PITCH ANGLE AT SLADE R00T*",F10. 5," DEGREES" //6X, 

1"C0N I NG ANGLE*"* F10. 5, IX, "DEGREES" //6X, 

^"FLAPPING COEFFICIENT, Ai »" ,F 10 L 5, IX, "OEGR EES "/ /6X, 

3"FLAPPIN3 COEFFICIENT, 31*",F10l5, IX, "DEGREES") 

PROGRAM STEP IE. 

START LOOP FOR AZIMUTH STATIONS AT WHICH OUTPUT QUANTITIES ARE 
CALCULATED. 

CTCAL-O. 

CMXCAL»0. 

CMYCAL-O. 

DO 510 1*1,24 
P900U)*15.M 1-1. ) 

P 3 3* ( I-1.)*PI/12. 

C P L *CGS ( P SO ) 

SPi*SIN(PBO) 

CP2*CUS(2.*P8Q) 

SP2*SIN(2.*P30) 

Fl(I)*-THC-TW*RQ/(i.-R0)-2.*31*CPl+(2.*Al+4.*MU*TW/(l.-R0))*SPl 
F2(I)*CPl*(2.*MU*Tw+RO/(l.-ROI>2.+MU*THC)*2.*ttU*AO*SPl 
l + 2.*MJ*81*CP2-(2.*MU*Al + iU*'1U*TW/(l.-R0) >*SP2 
F 3 { I ) » aC-C PI* ( 4,*MU*Ti*/(l.-RQ) + 2«*41)-2.*91*SPl 

FG1*G0(1) 

FG2=G0 { 2 ) 

FG3»G0( 3) 

F N 1 * 0 . 

FN2*0. 

GO 520 11*2,5 
FNl«FNl+FIl(Il)+GO(Il) 

Fn2»FN2+FI2(I1)*G0(I1) 

DO 520 12*1, NHM 
CI2-CGS(I2*P30) 

S I2*S IN ( I2*P30) 

I F ( I l .GT. 2 ) GO TO 530 
FGi«FGi*SC(l, I2)*CI2+GS( l, 1 2) *SI2 
FG2*FG2+GC( 2, 12 ) *C I2 + GS ( 2 , 1 2) *S 1.2 
FG3-FG3+GCI3, 12 ) *C 1 2 + GS ( 3, I 2) *S 1 2 
530 FNi»FNl»FIl(Il)*(GC(Il,I2)*CI2*GS<Il,I2)*SI2) 
Fn2*FN2+FI2(I1)*(GC(I1,I2)+CI2+GS(I1,I2)*SI2) 

520 CONTINUE 

FLT1( I ) *-PI*C l.-RO)*F II ( 1 ) *FG1/ AR+PI* ( 1 ,-R0 ) *(1.-R0 )*FG1*( 1 «-R0) * 
l( l.-RO)*FG2/3.+(l.-RO)*FNl/(2.*AR)+PI*<l.-RO)*<l.-RO)*(l.-RO) 
2*F2(I)/(4.*AR*AR)+PI*(l.-RO)*(l.-RO)*(l.-RO+RO)*F3(I)/(3.* 
3AR*AR)+PI*(l.-RO)*(l.-RO)*<l.-RO)*<l.-RO)*<Fl(I)*2.*TW*(l.+RD) 
4/(l.-R0))/(16.*AR*AR*AR) 

FLT21 I)*FLT1{ I)/(CT*PI/N9) 

FMTI I) ■-?!*( l.-R3)*F12(l)*FGl/AR4PI*(l.-R0)*C.5*(l.-R0*R0)»(l.-R0) 
i*(l.-RO)/6.)*FGl+(l.-RO)*U.-RO*RO)*FG2/6.+(l.-RO)*(l.-RO) 
2*(l.-RO)+FG3/lG.+(l.-R3)*FN2/(2.+AR)*PI*(l.-SO)*(l.-RO)* 

3 ( 1.-RO+RO)*F2(I)/(0.*AR*AR)+PI*(1.-RO)*(1.-Ru)*( 1.-RO*RO+RO) 
4*F3(I)/<i2.*AR*AR)+PI*(l.-RO)*(l.-RO)*(l.-RC)*{l.-RO*RO)*Fl(I) 
5/(32.*AR*AR*ARI*?I*TW*C l.-RO)*(l l .-RO)*Cl.-RO*RO*ROI/C 12.*AR 
6+AR4AR) 

C0L(I)*FMT(I)/FlT1(I) 
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CTCAL*CTCAL*FLTi(I)/24. 

CMXC4L-CMXCAL+FMT(I)*SPl/24. 

CMYCAL*CMYCAL-FMT(D*CPl/24. 

START LOGP F Oh? RADIAL 5 T AT IONS AT WHICH OUTPUT QUANTITIES APE 
CALCULATED. 

DO 510 1 1* 1» N SP 
F23«F2<I)+k30(I1)*F3< I) 

Z3-2.*(RBG<I1)-.5M1.«-RJ) J/U.-RO) 

UT-R30I ID +MU*SP1 
S3Z*S0RT(l.-ZB*Zd) 

FN0»FUN3(R0,AR,RB0(I1)» 1) 

FN1-0. 

FN2-0. 

FNJ-Q. 

00 540 12-2,5; 

NL- 12-1 

PN1-PNM! 12-1, 1, Z3) 

PN2-FUN3UO, AR,R30( II ) , 3) 

FN1-FND3CI I2)*PN1 
FN2-FN2+G0(I2)*PN2 
FN3»FN3+3C(I2)*NL*(NL*1.)*PN1 
DU 540 I 3*1, NHM 
C 1 3-C OS ( I 3*P 8 0 ) 

SI3»SIN(I3*P30) 

FN1-FN1+PN1*(GC(I2, I3)*CI3+GS(I2, 1 3 ) *S I 3 ) 

FN2-FN2+PN2*! GCI 12, I3)*CI3+GS(I2,I3)*3I3) 
FN3-FN3+NL*(NL*1.)*PN1*(3C(I2,I3)*CI3+GS(I2»I3)*SI3) 

540 CONTINUE 

GST(I,Il)=Afi*((l.+Z3)*FGl+3QZ*FNl/(2.*PI)) 

PDIPt I, II )--FNC*F5DFN2/ ( 2,*PI ) 

DGS2I I, ID— AR*fn3/(2.*PI*30Z) 

FL1( I, I1)-<1.-RO)*PI+(POIP( I»IIDGST(I,I1D<1.-R0)*F2 3/ 
l(4. + ARD(l.-R0)*<l.-R0)*(Fl(ID4.*Tw*RB0(Il)/(l.-R0)) / 

2 ( 16 . * AR* AR ) ) / AR 
FL2(I,I1)-FL1II,ID/(CT*PI/NB> 

FM( I, 1 1 ) -GST ( I, Il) + ( 2.*0GS2( I, ID + (1 ,-R0)*( l.-RO )*<F1II)+4.*TW 
1*R30(I1I/(1.-R0)))/(32.*AR*AS) 

FM( I, ID-FMt I,I1)*PI*(1.-RO)*(1.-RO)/(4.*AR*AR) 

FM( I * II ) *FM{ I , 1 1) ♦ ( 1 .-R0 ) *FL1 ( I, 1 1 ) / C 4. *AR ) 
XCP(I>Il)-.25+FM{I,Il)*AR/(FLl( I,I1)*<1.-R0) ) 

510 CONTINUE 

CTCAL-CTCAL*NB/PI 

CMXCAL«CMXCAl*NB/PI 

CMYCAL«CMYCAL*NB/PI 

WRITE(6,541)CTCAL»CMXCAL»CMYCAL 

541 FORMAT! /6X, "COMPUTED THRUST COE FF IC I EN T-", E 10 . 4/ / 6X , 

1 "CO MPUT E 0 MOMENT COEFFICIENT ABOUT ROTOR X-AX I S - " , £ 10 . 4/ /6X, 
2"C0MPUT£D MOMENT COEFFICIENT A30UT ROTOR Y-AXIS- ", E 10 . 4 ) 

C 

C PRINT OUTPUT IN TABULAR FORM. 

C 

WRITE (6, 10) 

-RITE! 6, 550) 

550 F0RMAT(//6X,"TA3LE 1 - SECTIONAL L I FT/ ( R HG* ( 0ME3 A ♦* 2 ) * ( R 1**3 ) ) »/ 6X 
U49UH-I//) 

WRITE(6,560)(R30(I),I»i,NSP) 

56C FORMAT! /12X, "R/Rlt ", 5 ( 1 10 . 4, IX ) 1 
WR I TE ! 6* 5 70 ) 
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570 F0RMAT(/7X*"PSI") 

00 5d0 1*1,24 

WRITE!6»590)PBOD! I),! FL11 1,11), IW.NSP) 

59C FORMA T( /6X,F5.1»8X, 5 (=10. 4# IX) ) 

580 CONTINUE 

WRITE!6,1C) 

WRITE !6,6C0) 

600 FQRMATI//6X, "TABLE 2 - SECTIONAL LIFT*R1/THRUST PER BLA0E"/6X, 

144! H-) // ) 

WRITE (6, 560) !RB0!I>,I*1,NSP) 

WRITS!6,570) 

00 610 1*1, 2<, 

wRITE!b,590)PB0D!I),<FL21I, II), I 1*1 » NS P) 

6X0 CONTINUE 

WRITS(6,10) 

WR ITS ! 6, 620) 

620 FORMAT! //6X, "TABLE 3 - SECTIONAL PITCHING MOMENT 7 ! RHQ* ! GMEGA + + 2 )♦ ! 
1R1**4 ))"/6X, 60! 1H-I/16X,"! ABOUT QUARTER-CHORD)"//) 
WRIT£!6,560)(RB0(I),I"1,NSP) 

WRITE(6»570) 

00 630 1*1,24 

WRITEI6* 590)P800(I)» !FM! I» ID# Il*l»NSP) 

630 CONTINUE 

WRIT£!5»10) 

WRITE <6,640 

640 FORMAT! //6X»"TA3LS 4 - CENTER OF PRESSURE LOCATION FROM LEADING ED 
XGE(FRACTI0N OF CH0RD)"/6X»74!1H-) //) 

WRITE !6,560)(RB0!I)»I-l»NSP) 

WRITE!6»570> 

00 650 1*1*24 

WRITE <6, 590) PBOD! I ) » < XCP! I, ID, U*1,NSP) 

650 CONTINUE 

WRITEI6, 10) 
w R I T£ ( 6* 6e0) 

660 FORMAT! //6X,"TA8LE 5 - TOTAL BLADE LIFT* MOMENT ABOUT HUB AND RADI 
1AL CENTER OF LIFT'VoX* 70! 1H-) //6X* 

2"T0 TAL BLADE L I F T/ ! RHO* ! QME 3A+*2 ) * ( R1 **4 ) ) "/6X, 

3" TOTAL BLAOE LIFT/THRUST PER 3LADE"/6X, 

4"M0MENT ABOUT HUB/ 1 RHO* ! OME GA**2) *( R 1**5 )) "/6X* 

5"RAD I AL CENTER OF L I F T/Rl" 7 /7X* "P SI ", 8 X, "TOT AL BLADE LIFT"* 
65X,"M0MEiNT"*5X,"CENrER"/33X»"A3QijT HUB",4X»"0F LIFT") 

DO 670 1*1,24 

WRIT£<6,680)P800!I),FLT1(I),FLT2(I),FMTII)»C0LCI) 

680 FORMAT(/6X,F5.1»5X,4!E10.4,1X)) 

670 CONTINUE 

WRITE ! 6, 10) 

WR I TE ! 6, 690) 

690 FORMAT! //6X*"TA3LE 6 -SURFACE PRESSURE D IF FE RENT I AL/ I RHO* l O.MEG A»*2 
l)*!Rl*+2) )"/6X,64!lH-)//) 

DO 700 1*1,24 
wRIT£!6,710)PB00!I) 

710 F0RMATI/6X, "AZIMUTH ANGLE •"* F 5 . 1,! IX , "D EGRE ES"/6X* 27 1 1 H-) ) 
WRIT€!6,560)!RB0<I1), Il»l,NSP) 

WR I TE ( 6* 720) 

720 FORMAT! 11X,"X/C") 

00 730 Ii»l*NSP 
F23*F2( I)+RB0(I1)*F3!I) 

DO 740 12*1,10 
CCH«2.*X0UT!I2)-1. 

SCh*SQRT!1.-CCH*CCH) 
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PCJT( II. 12)»-GST< I»I1)*SCH/ (l.+CCH)+2.*SCH*<P0IP(I.Il)*2.*GST 
HMD) +UGS2( X, I1)*SCH*SCH/<4.*AR+AR)+F23*(1.-RJ) *SCH/C2.*AR) 
2+ ( F 1 ( I)+4.*TW*R30l Il)/(1.-R0) ) * ( 1 .-RO ) * ( l.-RO ) *( 1 .-RO ) * 
3SCH*(1.+CCH)/(8.*AR*AR) 

P0UT(I1>I2)»2.*P0JT(I1.I2) 

740 CONTINUE 
73C Ca.NlINUc 

DO 750 I2-1.10 

WRITE(o»755)XQUT( 12). ( P0UT( II. 1 2 ) . 1 1-1, NSP ) 

755 F 3Ri1a T ( fcX» F8 • 5* 5X» 5 ( E 10 • 4» 1 X ) ) 

750 CONTINUE 
IPG*3 

IREH»I-(I/IPG)*IPG 
I F ( IREN .EO, 0 ) WR I TE ( 6. 10) 

70C CONTINUE 
STOP 

PRINT ERRQP MESSAGE IF COEFFICIENT MATRIX IS SINGULAR. 

420 WRITEls.760) 

7fcO F0RMATC/6X, "COEFFICIENT MATRIX IS SINGULAR") 

STOP 

END 

SU3R0UT INE OMIN 1 M d, L AM. 03 . R30 » P30» PL IM, DM A X , I. P ) 

' CALCULATION OF AZIMUTH POSITIONS AT WHICH TRAJECTORY IS DIRECTLY 
OVER A 3LA0E » WITHIN A DISTANCE DM AX . 

REAL MU. LAM 
DIMENSION P ( 20 ) 

Y(X»XI.X2.Y1.Y2)»Y1MY2-Yl)*( X-Xl ) / ( X2-X1) 

1-0 

P1*0. 

X 3 1* 0 . 

P2—0.2 

10 X32»R30*SIN(P2«-Da)+MU + SIM(P2+-P30 + 08)*P2 

I F ( PI .NE. 0 . ) GO TO 20 
30 P1«P2 

IF( °1 .LE. PLIMIRETURN 

P2-P2-0.2 

X31-X32 

GO TO 10 

20 TEST »X31*X32 

IF ( TEST .GT. 0. ) GO TO 30 
PC-YI0..X81.XB2.P1.P2) 

X8C*R80*SiN(PC*D8)*MU*SIN(PC+P80t08)*PC 
RaC-RaO*COS(PC+08)+MU*CaS(PC+PBO*OB)*PC 
D»SQRT(X3C*X8C+LAM*LAM*PC*PC+(RBC-RB0)*(RBC-RB0) ) 

I F ( 0 .GT. DM AX ) GO TO 30 

I-I+l 

P C I ) - PC 

GO TO 30 

END 

SUBROUTINE TRAJIRO.AR.MU.LAM.RBO.PBO.OB.Pa.ZS) 

C THIS SUBROUTINE CALCULATES THE PARAMETERS FOR THE LINEARISED 
C TRAJECTORY CORRESPONDING TO A GIVEN COLLOCATION POINT. 

PEAL MU. LAM 

COMMON /TR A Jl/R, SX.CX.SHP.CHP.ST.CT.SHE.CHE.SP.CP 
X8»R30*SIN(Pa+D8-PB0)+MU*( P 3-P30 ) *S IN ( P 3+0 B ) 

Y3»LAM+(P3-P30) 



Z3*-t 5+ ( l.+KO)+R93*COS{P8+Da-P30»«-MU*(P8-Pao)*COS<P3+OB) 
R»S3RT(X8*X8+YB*Y3) 

SX-Y3/R 

cx«xa/R 

XS»2.*XB/(1.-R0) 

YS«2.*YB/(1.-R0) 

ZS*2.*ZB/C1.-R0> 

RS«2. / { 1 . -RO) 

R1*SQRT(RS*RS + (1. + ZS) + (1.+ZS) ) 
fc2*S3RT(RS*kS+(l.-ZS)*(l.-ZS)J 
CHP »( R1 *-R2 )/ 2 • 

SHP»SORT(CHP*CHP-l.J 

CWS/CHP 

ST-RS/SHP 

R3*SJRT( (XS*AR-1.5)*C XS* AR-1 . 5 1 +YS* YS* AR*4R ) 

R4*SQ*T( ( XS*AR+.5)*< X$*4R*.5)+YS*YS*AR*AR> 

CHE* ( R3+R4) ! 2 , 

I F ( CH E .LE. l.)CrlE»A3S( (R4-R3)/2«) 

SHE»SQRT(CHE*CHE-1.) 

CP*( XS*AR-.5) /CHE 

5P*4R*YS/SHE 

RETURN 

END 

FUNCTION FUN 1 ( R 0» AR» I ) 

C THIS SUBPROGRAM RETURNS THE VALUES OF VARIOUS FUNCTIONS REQUIRED 
C IN FUNCTION FUN2. 

CDMMON/MAINl/N* PI 

COMMON/ TRAJ1/R, SX,CX» SHP> CHP» ST# CT» SHE » CHE , SP» C P 
GO TO <10, 20, 30, 40,50, 60), I 
10 F1-SHP*SHP*(CHP-CT) 

F2*(CHP-CTJMCHP-CTJ 

F3«(CHP-CT)*(CHP+CT) 

F<,«(SHP*SHP*CT*CT-CHP*CHP*ST*ST)f F1-ST*ST*(CHP+CT)/F2 

FUN1»-2.*(SX*SX*F4/F34-CX*CX/F1)/(1.-R0) 

kETURN 

2C PN1*PNM(N,1,CT) 

PN2“PNM { N,2 »C T) 

QN1*QNM(N,1,CHP ) 

QN2*QNM ( N, 2, CHP ) 

F1«1./(ST*SHP> 

F2«SX*SX/(SHP*SHP+ST*ST) 

FUN1»(F1*PN1*QN1+F2*( ST*CHP*PN1 *QN2-C T*SHP*PN2*QN1 ) )/ (PI*(1,-R0) ) 
RETURN 

30 F1»$HP*SHP+ST*ST 
F2»F1*SHP*SHP 

F3»CHP/F1-CHP*CHP*CHP*ST*ST/F2 

F3«F3-2. *CHP*ST*ST*(C HP *C HP ♦CT*CT)/(F1*F1) 

F3*F3*SX*SX/F1 

FUN1-2.*(F3+CX*CX*CHP/F2) /(l.-ROI 
RETURN 

40 FUN1- (l.-RO>*(CX*CX-SX*SX)/(R*R) 

RETURN 

50 F-1«SHE*SHE*CP*CP+CHE*CHE*SP*SP 
F2*(CHE+CP)+(CHE+CP» 

FUNl»2.*AR*(SHE*CPkSHE*CHE+(CP*CP“SP*SP))/((l.-RO)*Fi*F2) 

RETURN 

60 F1»SHE*SHE*CP*CP+CH£*CHE*SP*SP 
F2»SHE*CP*CP-CHE*SP*SP 
FUNI»2.*AR*CCHE-SHE)*F2/( (1.-R0)*F1J 
RETJRN 
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FUNCTION FUN2 l RO. AR, f W, MU, L AM, R 00, P BO, OP, 0 B, X, I ) 

C THIS SUBPROGRAM SETS UP THE INTEGRAND FOR THE INTEGRATION RE3UIRED 
C FOR THE INDJCED VELOCITY IN THE MAIN PROGRAM. 

REAL MU> LAM 
CGMMQN/r.4INl/N,PI 

COMMJN/FUN21/Fl,F2»F3,F4»F5»F6» ZBS 
P-X+D3 

IF ( I .S3. 2) GO TO 20 
IF ( I .GT. 2) GO TO 40 

CALL TRAJCRO, AR,MU,LAM,R30,Pd0,D3,X,ZS) 

ZBS-ZS 

F1-FUN1CR0,AR,1) 

F3-FUN1(R0,AR,3) 

F4«FUN1U0»AR»4) 

F?*FUNl(R0,AR»5> 

F6*FUN1(R0»AR,6) 

FUN2»-F1 

IF(ABS(Z33) .GT, I • ) GO TO 10 
FUN2»FJN2-(l.+Z9S)*F4/2. 

10 IFIUB .E3. C. .AND. X .GT. ( P30-0P ) ) RE TURN 
I F ( A3S ( ZBS ) .GT. 1. (RETURN 
FUN2-FUN2+AR*(1 . + Z3S)*F5 
RETURN 

2C F2*FJNiUC,AR,2> 

FUN2»-F 2 

IF ( A3S ( ZBS ) .GT. 1.IG0 TO 30 
PN1*PNM IN. 1, ZBS) 

FJN2»FUN2-PNl*5Q*ril.-Zt»S*ZBS)*F4/(4.*PI) 

30 IFOB .63. 0. .ANO. X .GT. ( PBO-OP )) RE TURN 
IFUBS(ZBS) .GT. 1. (RETURN 

FUN2»FUN2+AR*S0kT( l.-Z3S*Z3$)*PNl*F5/(2.*PI J 
RETURN 

40 I TEST *1-2 

za*< i.-po)*zbs/ 2 . 

KB»Z3+.5*(1.+R0) 

GO TO (50»60»70»80»90)» ITEST 
50 FN2«2.*MJ*COS(P) 

FN3*0. 

GO TO ICO 

60 FN2«2.*MU*SIN(P) 

FN3 » 1 . 

GO TO 100 

70 FN2— 2.*NU*SIN<2.*P) 

FN3 — 2.+C0SIP) 

' GO TO 100 
80 FN2»2.*MU*C0S(2.*P) 

FN3»-2 . *S IN t P ) 

GO TO 100 

90 FN2— TV*(2.*MU*C0S<P)*MU*MU*SIN(2.*P() /tl.-RO) 
FN3»4.*MU*TH*C0SIP)M1.-R0) 

ICO FUN2*— ( l.-P0)*F3*{FN2>RD*FN3) /( 4.*AR*AR) 

FUN2»FUN2*(1.“R0)*(1.-R0)*FN3*F1/C8»*AR*AR) 

IFCA3S(ZSS) .GT. i.lGO TO 110 
FUN2»FUN24-(l.-R0)*(FM2 + K3*FN3)*F4/( S.*AR*AR) 

110 I F ( 03 .63. 0. • ANO . X .GT. ( P 30-0P )) RETURN 
I F € A 3 S ( Z3S ) .GT. 1. (RETURN 

FUN2»FUN2-(1.-R0)*(FN2»R3*FN3)*F6/(2.*AR( 

RETURN 

ENO 



FUNCTION FUN3(R0,AR,X,I) 

C THIS SUBPROGRAM SETS UP THE INTEGRAND .FOR SOME SPANWISE INTEGRALS 
C REQUIRED IN THE MAIN PROGRAM. 

COMNQN/MAINl/.N»PI 
ZBS*2.*(X— ,5*(1.*RQI 1/(1. -RO) 

Rl«S3RT(l./( 16.*AR+AR)«-.2 5*(l.*ZflS)*(l.*Z8Sn 
R2«SQRT(l./( 16.*ASA4R)«-.25*{l.-ZaS)*(l.-ZBS)> 

CHPI-R14-R2 

SHP1-SGRT(CHP1+CHP1-1. ) 

CTW3S/CHP1 

ST1«1./(2.*AR*SHP1) 

I F C I .Lt. 2 ) GO TO 10 
PN1«PNM(N,1»CT1) 

0Nl»QNi1(.N,l,CHPl) 

GO TO 20 

10 F1»ST1/(SHP1*(CHP1-CT1I) 

IF ( I .EQ. 1IFUN3-F1 
I F ( I .50. 2IFUN3-X+F1 
RETURN 

20 F2«PNi*QNl 

I F ( I .EG. 3IFUN3-F2 
IF( I . EQ. 4IFUN3«X*F2 
RETURN 
END 

FUNCTION PNM(N»M»XI 

C CALCULATION OF ASSOCIATED LEGENDRE FUNCTIONS PNM . 

C RANGES 0 .Lc. N .LE. 4, 0 .LE. N .LE. 3, A3S(X) .lT. 1. 

IF(n .LT. 1 .OR. N .GT. 4 ) GO TO 10 
I F ( M .LT. 0 .OR. N .GT. 3IG0 TO 20 
IFUSSIX) .GE. l.IGO TO 30 
a X *.> QRT <1.-X*X) 

IF(M .EG. 1 > GO TO 40 
IFCM .£0. 2 ) GO TO 50 
I F ( M .EG. 3 ) GO TO 130 
GO TO (60,70, 80, 90), N 
60 PNN-X 
RETURN 

70 PNM».5*(3.*X*X-1.) 

RETURN 

80 PNM».5*(5.*X+X*X-3.*X) 

RETURN 

90 PNM«( 35.*X*X*X*X-30.*X*X*-3. 1/8. 

RETURN 

40 GO TO (IOC, 110, 120, 1301, N 
ICC PNM-SX 
RETURN 

110 PNN«3 . + X*SX 
RETURN 

12C PNM».5*SX*(15.*X*X-3. I 

RETURN 

130 PNM».5*SX+(35.*X*X*X-15.*X) 

RETURN 

50 GO TO (140,150, 150, 170), N 
140 PNM*0. 

RETURN 

15G PNM»3.*SX*SX 
RETURN 

160 PNM«15.*X*SX*SX 
RETURN 

170 PNN-.5*SX*SX*(105.*X*X-15.) 
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RETURN 

160 GO TO ( 190, <!00» 210,220) »N 
190 PNM«0. 

RETURN 
200 PNi1*0. 

RETURN 

210 PNM»15.*SX*SX*SX 
RETURN 

220 PNM«105.*X*SX*SX*SX 
RETURN 

10 RRITEC6,il)N 

11 FQRMAT(6X,*N-*,I5,1X,*INVALID N IN PNM+) 

STOP 

20 WRITc(6,21)M 

21 FDRMATC6X»*M«*» I5»1X»+IN VALID M IN PNM* ) 

STOP 

30 URITEC6.3DX 

31 F3R'UTC6X,RX»*,E10.4,1X#*INVALID X IN PNM* ) 

STOP v. 

END 

FUNCTION 3NM(N»N,x) 

CALCULATION OF ASSCCIATEO LEGENDRE FUNCTIONS 3NM •" 

RANGE* 1 .LE. N .LE. A, 1 .LE. M .LE. 2» A9S(X) .GT. 1. 
ASYMPTOTIC EXPANSIONS USED FOR X .GT. 3. 

I F ( N .LT. 1 .0 ft. N . GT • A)GO TO 10 
I F { M .LT. 1 .OR. M .GT. 2)G0 TO 20 
IF<A3S(X) .LE. 1. ) GO TO 30 
SX«S3RT(X*X-i. ) 

ALX*ALOG( (X*1.)/(X-1. ) ) 

X 2* X* X 

X3-X2+X 

X A*X 3*X 

X5-XA+X 

X6* X5*X 

X7« X6*X 

X3»X7*X 

X9»X8*X 

X10»X9*X 

X11*X10*X 

X12-X11+X 

I r ( M . E3. 2) GO TO AO 
GO TO ( 50>60»70»80).N 

50 IF t X .GT. 3. ) GO TO 51 
QNM»SX*t,5*AL)HX/(SX + SX) ) 

RETURN 

51 0NM»-SX*{2./< 3.*X3)+A./(5.*X5)+6./(7.*X7>+6./(9.+X9>) 

RETURN 

60 I F ( X .GT. 3. ) GO TO 61 
QNM»SX*(1.5*X*ALX-(3.*X2-1.»/(2.*SX*SX)-1.5> 

RETURN 

61 <3NM»-SX*(2./(5.*XA)+4./(7.*X6)+2./(3.*X8»+8./(ll.*X10l) 
RETURN 

70 I F ( X .GT. 3. ) GO TO 71 

3NM«SX*( 1 15. + X2-3. >*ALX/A.-< 5 . *X3-3« *X ) / < 2 . *S.X* S X ) -5 . *X I 
RETURN 

71 GNM»-SX*(6./(35.*X5)«-8./{21.*X7)+16./(33.*X9)+10./(13.*Xll)> 
RETURN 

8C I F ( X .GT. 3 . ) GO TO 91 

QNM»SX*{(35.*X3-15.+X}*ALX/A.-(35.*XA-30.*X2*3.)/(8.*SX*SX) 

l-105.*X2/3.*55./2A.) 
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RETURN 

81 GNM*-SX*(8./< t>3.*X6)*8./n3.*X3)+48./( 143.*X10) ) 

RETURN 

40 GO TO (93,100,110, 120), N 
<50 I F ( X .GT. 3 • ) GO TO -51 
GNH-2./ (SX*SX) 

RETURN 

91 QNH«SX*SX*(2./X4*4./X6*6./X8+8./X10) 

RETURN 

100 I F ( X .GT. 3. ) GO TO 131 

QNM»SX*SX*{1.5*ALX-&«*X/(SX*SX)*X*(3.*X2— l.)/(SX*SX*SX*SX) ) 
RETURN 

101 QNM*SX*SX+(8./( 5.+X5) *24 . / { 7. *X7 ) *16. / ( 3. *X9 ) *80 . /( 11 .*X11 ) ) 
RETURN 

110 I F ( X .GT. 3. ) GO TO 111 

QNM*3X*SX*(15.*X*ALX/2.-(15.*X2-3.)/(SX*SX)+X*(5.*X3-3.*X) 

l/(3X*SX*SX*$X)-5.) 

RETURN 

111 QNrt*SX*SX*(d./(7.*X6)*3./(3.*X8)*48./(ll.*XlQ)*llG./(13.*Xl2)) 
RETURN 

120 I F ( X .GT. 3. ) GO TO 121 

<jNM»SX*SX*( ( 1Q5.*X2— 15. )*ALX/4.-( 35 . * X3-15 . ♦ X ) 7(SX*SX ) 

1 + X*( 35.*X4-30.*X2*3. J / ( 4 . *S X*SX*S X*SX ) -105 . *X/4 . ) 

RETURN 

121 QNM*SX*$X*(16./(21.*X7)*64./(33.*X9)*480./(143.*X11)) 

RETURN 

10 WRITE(5,11)N 

11 FGRMAT(4X**N»*» I5,1X,*INVALID N IN QNM* ) 

STOP 

20 WRIT£(o, 2l)M 

21 FGRNAT(&X,*M.*,I5,1X,*INV4LID M IN QNM*) 

STOP 

3C WRITE(o,31)X 

31 FORMAT! 6X» + X« + * E10.4, 1X,*INVALID X IN QNM*) 

3 TOP 
ENO 

SU8R0UT INE TABSCHIX, N, XT, 11,12* INT) 

GIVEN AN ARRAY X, TO LOCATE THE POSITION OF A VALUE XT. 

IF INT«0, XT LIES BETWEEN XCU) AND X(I2). 

IF INT«1, XT IS GREATER THAN X(N). 

IF INT*-1, XT IS LESS THAN XU). 

DIMENSION X ( N ) 

11*0 
12*0 
NM*N— 1 

00 10 1*1, NM 

I F ( XT . G£ . X(I) .AND. XT .LE. X(I+l))GO TO 20 
10 CONTINUE 

I F ( XT .LT. X( 1) ) GO TO 30 
IF ( XT .GT. X ( N ) ) GO TO 40 

20 INT*0 

I F ( XT • EQ. X ( I ) ) GO TO 21 

I F ( XT , E3 . X{ 1*1) )G0 TO 22 

11*1 

12*1*1 

RETURN 

21 11*12*1 
RETURN 


67 



22 il*I2«I+l 
RETURN 
30 InT *-l 
RETURN 
•*0 I NT * l 
RETURN 
END 
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LISTING OF PROGRAM ASYMP2 
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PROGRAM MAIN (INPUT, OUTPUT* TAPE 5* INPUT, TAPE6*QUTPUT* 
1AF0ATA,TAPE1*AF0ATA) 

C 

C COMPUTATION OF THE UNSTEADY AIRLOADS ON A HELICOPTER ROTOR 3LADE IN 
C FORWARD FLIGHT* USING A SIMPLIFIED VERSION OF THE ASYMPTOTIC APPROACH 
C PROPOSED BY VAN HQLTEN. THE DIPOLE STRENGTH DISTRIBUTION ALONG THE 
C PLACE IS APPROXIMATED BY A PIECEWISE CONSTANT OR PIECEWISE QUADRATIC 
C REPRESENTATION. THE HELICAL TRAJECTORY OF THE FREESTREAM FLUID 
C PARTICLE KtLATIVE TO THE BLADE IS APPROXIMATED BY SUCCESSIVE STRAIGHT 
C LINE SEGMENTS. 

C FOR IDENTIFICATION OF PROGRAM STEPS* REFER USER'S MANUAL. 

C 

REAL MU*LAM* MCL*MINF, MLOC 

DIMENSION R (6), RBO( 5), OR( 5) *X ( 6 ) , FX1 ( 5 ) , FX 2 ( 5 ) , PI < 3) , P2 ( 5).IP( 5>* 
lFG{5),WK<59),aTI4>»GF<4*5),PMIN(20),G0(5),GC(5,5>,GS(5,5),G<24,:>), 
2A<59»59)»a(59»l),IPIV0T(59),p.j00{24),FLTl<2<*),FLT2(24), 
3FMTC24),COC<24)*GL<5>»:iR(5)»FH(24,5)*FL2<24,5),FMl24,5)* 
4XCP(24*5)*X3U.T(10) ,F23( 24, 5 )* POUT ( 5* 10 > 

DIMENSION CL(50»20),MCL(20I»ACLC50> 

CQMMON/MAINl/XBI* <3S,Y3S*R3I*R3S 

YL(XL»XL1*XL2,YL1. YL2 ) « YL 1* < XL- XL 1 > *C YL2-YLl)/( XL2-XL1) 

DATA X/-l.*-.9*-.&,3.,.5,l. / 

DATA NSP*NHM*NAZ*NC0F,0MAX/5» 5, 11*59*0.5/ 

DATA X OUT / . 0 5 * . 1* . 2* . 3* . 4 » • 5* . It • 9 * . 9 5* • 99/ 

PI-4.*ATAN(1. ) 

PROGRAM STEP 1. 

READ AND WRITE INPUT DATA AND ASSOCIATED QUANTITIES. 

READ(5»*)R0,AR»N3»TW»MU*ALR>CT*MINF 
REA0(5»*)N1*N2,.N3,N4 
1F(N1 .6-3. 1 )REA0(5,t)TH: 

IF(N2 .EQ. 0 ) Rfc AQ ( 5* * ) GAMA 
IFIN2 ,CQ. 1)REA0(5»*)AD 
IF ( N3 .EQ. 1)READ(5»*JAI 
IF(N^ .EQ. 1)REA0<5,*)81 
Rfc AQ( 5»*J (R( I ) » I *2 * 6) 

ISEL . EQ • 0 — PIECEWISE CONSTANT REPRESENTATION. 

ISEL .EQ. X — PIECEWISE QUADRATIC REPRESENTATION. 

REAO( 5,*)0P1D,0P20,UTMIN, ISEL»N AFD 

NAF0«0 — AIRFOIL TABLES NOT USED. 

NAFO-1 — AIRFOIL TABLE S USED. 

IFtNAFO .EQ. 0) GO TO 3 
REA0(1*1)NXL*NZL 

1 FORMAT! 30X, 212) 

READ(1*2)(MCL(I)*I»1*NXL) 

2 F0RMAT(7X»9F7.0) 

NLl*NXL/9 
NL2-NL1+1 
00 3 1*1, NZL 
DO 3 J* 1* NL2 
J 1* ( 4-1 ) *9*1 
42* J+9 

I F ( 4 1 .GT. NXL ) GO TO 3 
IF! 42 .GT. NXL) J2«NXL 
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I F { J .EQ. l)READ(l,A)ACL!I)»!CL(I»J3),J3*Jl»J2) 

A FORNAT!F7.0,<)F7.0> 

I F ( J .GT. l)READ(i,5)iCL(I,J3),J3*Jl,J2) 

5 FJRMAT!7X,9F7.0) 

3 CONTINUE 

6 CONTINUE 
DP1*DP10*PI/180. 

QP2-0P20+PI/180. 

WRITE (6, 9) 

9 FORMAT! 1H 1 ) 

WkITE!6,1O)R0,AR,N9,Tw»MU,ALR»MINF 

10 FORMAT! //6X» "ROOT RAOIUS/TIP RADIUS* R0/R1 •"» F 10 . 5/ / 6X , 

1" ASPECT RATIO* ",F10.5//sX»"NUM3£ii OF 8L AOS S»", I2//6X, 

2"LINEAR Twist !RaJT TO TIP) «",F10.5,1X,"DEGREES"//&X, 

3"FQRWAR0 SPEED/TIP SPEED*"* F10. 5//6X, 

A"RO T UR INCIDENCE !FORWARO TILT POSITIVE) *"»F10.5*1X» 
5"0EGREE$"//6X»"FREESTREAM MACH NUMBER* ", F10 . 5 ) 

R!l) «R0 
TW»TW ♦PI/180. 

ALR*ALR»PI/180, 

I F ( N1 . EQ. 0)THC*0. 

I F ( N 2 .EJ. 0 ) A0«0. 

IFIN3 .EQ. 0)A1*0. 

I F { N A .EO. 0 ) B 1*0. 

WRIT£!6,li)CT 

11 FORMAT! /5X, "THRUST COEFF IC I ENT*", F10. 5 ) 

IF !N1 .EO. 1 ) WR ITE ( 6, 12) THC 

12 FORMAT! /8X, "PITCH ANGLE AT BLADE R DOT* ", F10 . 5 , In , "OEGREES " ) 

IFIN2 .EO. 0)WRITE!S»13)GAMA 

13 F0RMATI/6X, "FLAPPING INERTIA COEFF IC I £N T*" , F 10. 5 ) 

IF ! N2 .EJ. 1 ) WR I TE ( t>, 1 A ) AO 

1A F0.RMAT(/6X, "CONING ANGL E » ", F 10. 5, IX , "OEGP E ES" ) 

I F ( N3 .EO. 1 ) WR ITE ( 6, IS ) A1 

15 FORMATI/SX, "FLAPPING COEFFICIENT, A 1-", F 10 . 5, IX, "0E3REE S " ) 

IFINA .EJ. DWRITEIS, 16)81 

16 FORMAT! /6X, "FLAPPING COEFFICIENT, B1 ■", F 10 . 5, IX , "OEGR EES " ) 
THC*THC*PI/180. 

AO* AO*P 1/180. 

Al ■ A1 ♦ P 1/ 1 80 . 

B1*31*P 1/180. 

LAM»MU*ALR+SQRT!.5W!-MUAMU+SQRT!MIJaMU*MU*MU+CT*CT J ) ) 
WRITE!6,2O)LAM,UTMIN,0PlD,0P20 
20 FORMATI/SX, "TOTAL INFLOW RATIO*", F10.5//6X, 

l"M IN IMUM UT»",F10.5,"!ZER0 LIFT CONDITION APPLIED BELOW THIS VALUE 
2 ) "/ /6X, "NORMAL AZIMUTH SPACING*", ^10. 5, IX, "DEGRE E S" / / 6X, 

3"REDUCED AZIMUTH S P AC ING-", F 10. 5, IX, "DEGRE ES ") 

IF ! I SEL .EQ. 0) WRITE!6,30) 

30 FORMAT! /6X, "PIECEWISE CONSTANT APPROXIMATION OF SPANWISE DIPOLE ST 
1RENGTH VARIATI0N"/6X,70!1H*)/) 

1FIISEL .NE. 0) WRITE! 6, AO) 

AO FORMAT! /6X, "PIECEWISE QUADRATIC APPROXIMATION OF SPANWISE DIPOLE S 
1TRENGTH VARI ATI0NV6X, 71! 1H*» /) 

SLCR-1. 

CLO-O. 

IF! N AFD .EQ. 0 ) WR I TE ! 6, A 1 ) 

Al FORMATI/SX, "AIRFOIL DATA TABLES NOT USED") 

IF1NAFD .EQ. 1 ) WR I TE ( 6, A2 ) 

A 2 FORMATI/SX, "AIRFOIL DATA TABLES USED") 

C 

C CALCULATE QUANTITIES NEEOED FOR TRAJECTORY SEGMENT ADJACENT TO 
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C THE BLADE AT THE COLLOCATION POINT. 

C 

ET41P«ALQG(<3.*SORT(5.))/2. ) 

C F1»C JSH( .5*ETA1P) / S I H t «5*£TALP) 

CF2*. 5-fcXP(-£T41P| 

CF3-ETA1P-.5+CF1 

C 

C CALCULATE QUANTITIES NEEDED FOR NEAR FIELD REPRESENTATION. 

C 

DO 50 I«l»5 
n-ACOsuin » 

T2«AC0SIX<I*1>) 

FX1(I)«<T2-T1-(SIN(T2)-SIN(TH) » / ( X ( I *1 J-X ( I ) I 

FX2( I ) »-( (T2-Tl)/2.-(SIN( 2. *T 2 ) -SIN ( 2 . *T1 ) ) / A . I / ( X C 1+ 1 ) -X ( I ) ) 

FX2(I)»FX2(I)*(1.-R0)/(2.*AR) 

50 CONTINUE 

DO 60 I»i#NSP 
S30(I>«U(IM-R(I + l)»/2. 

OR(I)-RII*l)-R( I) 

60 CONTINUE 

IFII5EL .EQ. 0) GO TO 70 
C 

C FOR PIECEWISE QUADRATIC RE PRE SENT A TI 0N» DETERMINE VALUES AT THE ENDS 
C OF THE SEGMENTS IN TERMS OF THE CENTRAL VALUES. 

C 

NSP M-NSP-1 
DO 30 I»1»NSPM 

aTm-3.«(i./DRm»i./0K(i+in 

I F C I .EQ. 1 ) GO TO 30 
3T<I)«3T(I)-1./(BT(I-1)*DRCI)+DR(D) 

BO CONTINUE 

DO 60 I*1»NSPM 
00 90 J«1,NSP 
GF(I»J)-0. 

IF ( J . c Q . I .UR. J .EO. (It’D )GFC I»J)*GF(I>J)+A, / B T ( I ) 

IF l I .EQ. I) GO TO 90 

GF< I» J)»GF< I» J)-GF(I-1, J ) /( 3T{I)*DR(I) ) 

90 CONTINUE 

DO 100 I-2.NSPM 

Il-NSP-I 

DO 100 J-1,NSP 

GF(Il,JI«GF(Il#J)-GF(Il+l,J)/(dT(Ili*DR(Il+lJ) 

ICO CONTINUE 
C 

C PROGRAM STEP 2. 

C BEGIN SETTING UP THE SYSTEM OF SIMULTANEOUS EQUATIONS, 

C 

70 L-l 

DO IIO J « i» N A Z 
C 

C SET THE AZIMUTH STATION FOR THE CURRENT COLLOCATION POINT. 

C 

P30-2.*J*PI/NAZ 
PBODI J)*360.*J/NAZ 
CP1»C0S(?30) 

SPl-SIN(PBO) 

CP2-C0S (2.+PB0) 

SP2»SIN(2.*PB0) 

DO 110 I»1»NSP 
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C SET THE RADIAL STATION FOR THE CURRENT COLLOCATION POINT. 

C 

UT*R3CII)+MU*SP1 
DO IIS M-l»NCOF 
A(L»M)*0. 

115 CONTINUE 

PROGRAM STEP 3. 

TEST THE TANGENTIAL VELOCITY AT THE COLLOCATION POINT# TO DECIDE 
WHETHER NORMAL VELOCITY BOUNDARY CONDITION OR ZERO LIFT CONDITION 
5 HOI) LO 36 APPLIED. 

I F ( U T .GT. UTMINJ30 TO 12D 
wRITE(6#130)RB0(II#Pd0D(J)#UT 
13 0 FORMAT! /6X#"R»"#F9.3#1X,"PSI»"» F3 . 3# "OEGRE E S" » IX * 
1 M UT*"»F9.3#1X#"ZE8J LIFT CONDITION APPLIED") 

3(L»1)*0. 

GO TO 1 AO 

120 DP«1.5*(1.-R0)/12.*AR*UT) 

PLIM--12.+RB01 I )*CP1) /SQRTIMU*MU*LAM*LAM) 

IFINaFO .EQ. 0)G0 TO 133 

PROGRAM STEP A. 

CALCULATE LIFT CURVE SLOPE FROM DATA TABLES FOR THE CURRENT 
COLLOCATION POINT# USING THE LOCAL INCIDENCE AND MACH NUMBER. 

MLOC»UT*M INF/MU 

AL0C*THC-TR*(R80(I )-R0) / C 1 .-R 0 ) +3 1*CP 1- Ai* S PI 
1— IMJ + A0*CP1 + LAM)/'JT 
AL0C-ALQC+1B0./PI 

CALL TA8SCh(MCL#NXL,ML0C# IMCL1# IMCL2# INT) 

IF ( InT .EQ. -1)IMClI*IMCL2»1 

IF ( INT .EQ. 11 IMCL1»IMCL2»NXL 

CALL TA3SCHI ACL >N ZL# ALOC » I ACL1# IACL2# INT) 

I F ( INT .EQ. 0 ) GO TO 131 
IF( INT .EQ. -1 ) GO TO 132 
IACL1-NXL-1 
I AC L2 *N XL 
GO TO 131 

132 I ACL1«1 
IACL2-2 

131 SLC1«ICLI IACL2# IMCL1J-CLII AC LI# IMCL1) ) /( ACLIIACL2) 

1-ACL! IACL1) ) 

CL01»CL(IACL1»IMCLU-SLC1*ACL(IACL1) 

SLC1«SLC1*180./PI 

SLC2« (CL< IACL2# IMCL2 ) -CL C I ACL 1# IMCL2) ) / ( AC L ( I ACL 2) 
l-ACLI IACL1) ) 

CL02«CLUACL1#INCL2)-SLC2*ACL(IACU> 

SLC2*SLC2*160./PI 

IFIIMCLl .EQ. IMCL2)SLCR«SLCl/< 2.*PI) 

I F ( IHCL1 .NE. IMCL2)3LCR"YHML0C»MCLIIMCL1)»MCL{IMCL2)# 
1SLC1#SLC2)/(2.*PI) 

IFIIMCLl .cQ. IMCL2)CL0»CL01 

IFIIMCLl .NE. IMCL2)CL0«YLIML0C»MCLIIMCL1)»MCLIIMCL2)» 

1C LC1, C L 02 ) 

133 CONTINUE 

IFISLCR .EQ. 0. )G0 TO 13A 
Tl — (CF1 + I1.— SLCRJ/SLCR) /'JT 
T2*MU*CP1*CF3*(1.-RQ) / 1 2 . * AR*UT*UT*DR I 1) ) 
T3»-(1.-R0)*CF3/I2.*AR*UT*UT) 


73 



o r> o o o 


S(L> U«-iU + ALR-UT*TW+(RdO( I >-R0) / ( 1 . -RO ) ♦UT + C LO/ (2.*t>I^SLCR) 

GO TO 140 
134 Tl-1. 

T2-T3-0. 

8<L,i)«-UT*DT*CL0/(2.*PI} 

140 CONTINUE 

SET UP THE CONTRIBUTION TO THE COEFFICIENT MATRIX OF THE TRAJECTORY 
SEGMENT ADJACENT TO THE SLADE. ALTERNATIVELY, SET UP THE ZERO LIFT 
CONDITION. 

DO 150 M*1,NC0F 
Ml* (M— 1 ) /NAZ+1 
M2*M-(M1-1)*NAZ 
IFCM1 • EO . (NSP+1 ) )GG TO 160 
M3-M2/2 
M4*M2-M3*2 
WT1*0. 

WT2*0. 

WT3-0. 

IF( JT .LS. UTMINJGQ TO 170 
IF ( Ml .EQ. I)wTl*Tl 
I F ( Ml . EQ. I)WT3«T3 
IFIISEL .EO. 0 ) GO TO 180 
I F ( I .LT. NSP)WT2*WT2 + T2*GF, (I»N1)/DR{M1) 

I F ( I .GT. 1) WT2*WT2-T2*GFI I-l, Ml ) /DR Ml ) 

GO TO 180 
170 I F ( Ml .EQ. I ) WTl*l. 

180 I F ( M 2 .£0. 1 ) A ( L, M ) * WT1+ WT2 

I F l M2 .GT. 1 .ANO. M4 .EQ. O)AIL,M)»(WT1*WT2)*C0S(M3*PBO) 
l-M3*WT3*SIN(M3*PBO) 

I F ( M2 .GT. 1 .ANO. M4 .GT. 0 ) A { L > M ) * ( WT1+WT 2 ) *S I N ( M 3* P80 ) 
1*M3*JT3*COSIM3*P80) 

GO TO 150 

160 IF ( SLCR .EQ. 0. )G0 TO 150 

I F ( UT .GT. UTMIN)WT»-(1.-R0)*CF2/(2.*AR*UT) 

IFCJT .LE. UTMINJwT — { 1 . -R D) / ( 4 . * AR ) 

I F ( M2 .EQ. 1 )ML,M)-WT*<2.*MU*CP1) 

I F ( M2 .EQ. 2)A(L»M)*WT*( 2.*MU*SP1+R8Q( I) ) 

I F ( M2 .EQ. 3 ) M L * M ) * W T+ { - 2 • + M U * S P 2— 2 • * R B 0 ( I ) + C P 1 ) 

I F < M 2 .EQ. 4)A(L,NJ*WT*(2.*MU*CP2— 2.*RB0(I)*SP1) 

150 CONTINUE 

I F ( SLCR .EQ. 0 • ) GO TO 191 

8(L»l)"B(L»l)-WT*TW*(2.*MU*R0*CPl-NU<‘NU*SP2-4.+HU*CPl+R30(I) )/ 
l( l.-RO) 

IF ( UT .GT. UTM IN ) GO TO 190 
191 L«L*1 

GO TO 110 
C 

C PROGRAM STEP 5. 

C START LOOP FOR NUMBER OF BLADES. 

C 

190 DO 200 I3L*1,NB 
I8LD* I3L 

D3*2.+PI*(I8L-1)/N3 

C 

C CALL SUBROUTINE TO DETERMINE AZIMUTH POSITIONS ALONG THE TRAJECTORY 
C AT WHICH THE TRAJECTORY IS CLOSE TO A BLADE • 

C 

I F ( DP 1 .NE. D P2 ) C ALL OMIN{MU>LAM»OB,RBO(I),PBO,PLIM,DMAX, 
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1IMIN,P1IN) 


DEFINE THE AZIMUTH INTERVAL FOR THE FIRST TRAJECTORY SEGMENT. 

J 1- 1 
K1«0 
P32-0. 

I F ( 1 3 L .60. 1 ) P31*-DP 
I F{ 1 9L .5T. 1)P31 — DPI 
210 CONTINUE 

DO 220 1 1-1* 5 
FG( II >*0. 

220 CONTINUE 
FCNF-0. 

FTHNF-0. 

F AONF *0 . 

F 4INF*0. 

F31NF-0. 

PROGRAM STEP 6. 

CALCULATE SLOPE AND INTERCEPT COMPONENTS FOR CURRENT TRAJECTORY 
SEGMENT. 

XB1»R30(I)*SIN( Pal+DBH-MU*P31*SIN(Pai*P30*DB) 

X32*R30(I )*SIN( ?92 + D3) + !iU*PB2*SlN(P62*P80*DB J 
R31»R30( I) *COS( P81 + D3 ) ♦ TJ*P ol *C OS < P31 + P30 + 03 ) 
R32»R30(I)*COS(P32 + D3H-MU*P32*CQS(P32 + PB0<-03l 
X 3 S - ( X32-XB1) /(PB2-P31) 

X3I»X31-X3S*PB1 

IFIP32 .EQ. 0. .AND. (D3 .EO. 0. .OR. D3 .£0. PI))X9I*0. 
RBS*(S32-RS1)/(PB2-P31) 

R B I* R 81-R3S + P31 

Y3S*LAM 

NSPP-NSP+1 

START CALCULATION OF FAR FIELD CONTR IBUT I ON. 

DO 230 11*1, NSPP 
12*11-1 

CALL FFINT<PB1,PB2,R( ID, ISEL, FG1 , FG2, F G3I 
I F ( I 1 .Ei. 1 ) GO TO 240 
I F ( ISEL .£0. 0 ) GO TO 250 

FF1 — ( 1.-R0)*(FG1-FG1M)/(4.*AR*0R< I2)*0R< 12) ) 

FF2— (1.-R0)*(R<I1)*FG1-R{I2)*FG1M) /( 4 .*AR*DR< 12 ) *DR C 12 ) ) 
FFJ— I1.-R0)*(FG2-F52M)/(4.*AR*0R(I2)*0R(I2)) 
FF4*-(1.-RO)*(FG3-F33M)/(4.#AR*DR{I2)*OR(I2)) 
T1M«(-2.*(R30(I2)+R(I1) )+2.*RBl)*FFl*2.*FF2*2.*RBS*FF3-2.*FF4 
T1*(8.*R30(I2)-4,*R3I ) *FF 1-4. *F F2-4 »*RBS*FF 3 + 4. *F F4 
T1P*(-2.*(RB0(I2)«-R( 12) ) *2 • *R 8 1 ) * FF1»2 . *FF 2 + 2 • *RB S*FF 3-2 » *FF 4 
DO 200 13*1, NSP 

IF ( 12 . N£ . 1)FG(I3)»FG(I3)+T1M*GF(I2-1»I3)/0R(I3) 

I F I 13 .EO. I2)FG(I3)»FG(I3)+T1 | 

IF ( 12 .NE. NSP)FG(I3)«FG(I3)+T1P*GF(I2,I3)/DR(I3) 

260 CONTINUE 
GO TO 240 

250 FG(I2)»FG(I2)-(1.-R0)*( FG1-FG1M ) / ( 4. * AR ) 

240 FG1M-FG1 
FG2M-FG2 
FG3M-FG3 
230 CONTINUE 


75 



n n n o o o o o o o 


END CALCULATION OF FAR FIELD CONTRIBUTION. 

CALL SUBROUTINE TO OIVIOE CURRENT TRAJECTORT SEGMENT INTO SUB- 
SEGMENTS ALIGNEO *1TH 3PAUISE SEGMENTS ALONG THE BLADE. 

CALL SUBIVL(R81.R32»PB1»P32>R»IP»P1»P2) 

CBl-COS (P82*P30+03) 

C82-C0S (2.*{Pd2+P30*08) ) 

S81«SIN(P324-PB0+03) 

S32«SlN(2.*(Pa2*PB0+03) ) 

FC-TWM2. *MU*R0*C3 1— MU*MU*S 82-'V. *MU*R B2*CB l)/CI.-RO> 

FTH»2.*MU*Ct)l 

FAO-2.*MU*S31+R32 

FAl--2.*.NU*S82-2.*R82+C81 

Fdl»2.*MU*C82-2.*R32*S8l 

START CALCULATION OF COMMON PART AND NEAP’ FIELD CONTRIBUTIONS. IN 
THE FOLLOWING LOOP* T HE FIRST PASS CALCULATES THE COMMON PART. AND 
THE SECOND PASS THE NEAR FIELD. 

DO 270 11-1,2 

I F ( 1 1 .EJ. 2 .AND. I3L .ED. 1 .AND. P B 2 .cQ. 0.)G0 TO 290 

DO 270 I2-1.NSP 

FG1-0. 

FG2-0. 

FG3*D. 

FG4-0. 

I F ( IP ( 1 2 ) • EQ • 0 ) GO TO 270 
I F ( I I .ED. 2 ) GO TO 290 

CALL CPINT(Pld2)»P2(I2> » ISSL»FG1»FG2.FG3) 

IFdSEL • E 0 • 0) GO TO 300 
F31--(l.-R0)*FGl/(2.*AR>DRd2)*DRd2)) 

FG2 — (l.-R0)*FG2/(2.*AR*0R(I2)*DR(I2) ) 
FG3»-(l.-Rd)*FG3/(2.*AR*0R(I2)*0R{I2) ) 

GO TO 310 

300 FGd2)-FG(I2)-(l.-R0)*FGl/(2.*AR) 

GO TO 270 
290 00 320 13-1.6 

CALL NFINT(P1(I2)»P2( 12). (l.-RO)*X( 13) / (2.*AR>» IS EL » F Nl, FN2, F N 3 ) 
I F ( 1 3 .ED. 1 ) GO TO 330 
FG4-FG4+FX2( 13-1 ) ♦ ( FNI-FN1M ) /PI 
IFdSEL .EQ. 0 ) GO TO 340 

FGl-FGH-FXHI 3-1 1 * ( FNI-FNIN ) / ( PI*DR(I2)*QR (12) ) 

FG2-FG2 + FXKI 3—1 ) * ( FN2-FN2M )V( PI*DR ( 12 ) *0R (12) ) 

FG3-FG3 + FXK 1 3-1 ) * ( FN3-FN3N )VfPI*DR (12 )'*0R ( 12) ) 

GO TO 333 

340 FG(I2)-FG(I2)+FX1(I3-1)*(FN1-FN1N)/PI 

330 FN1M-FN1 
FN2M-FN2 
FN3M-FN3 
320 CONTINUE 

IFdSEL .EQ. 0 ) GO TO 270 

310 T 1(1-2. *( (RBI— R( I2tl ) ) *( R8 I*RBO( 12 ) ) ♦FGl + ( 2 . *RBI-RBO ( 1 2 ) 

1-R( 12*1) )*RBS*FG2+R3S*R3S*FG3) 

Tl.-*.* < (RBI-R(I2) )*(RBI-R( 1 2*1 T) *F Gl* ( 2 . * R6 1-R ( I2)-R (12*1) ) 
1+RBS*FG2*R3S*R9S*FG3) 

TIP-2 .+ ( ( RBI-R( I2))*(R3I-R30(I2))*FG1+(2.*R3I-R30(I2)-R(I2)) 
1*R3S*FG2+«BS*RBS*FG3) 

00 350 I3-I.NSP 

I F ( 12 .GT. l)FG(I3)»FG(I3M-Tlrt*GF(l2-l.I3)/0R(I3) 
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I F ( I 3 .£0. I2)FG(I3)»F3(I3)+Ti 

I F ( I 2 .NE. NSP)FG(I3J«FG< 13) »T1 P*GF ( 12* 1 3 > /Ok ( 1 3 ) 

350 CONTINUE 

I F ( I 1 .EO. 1 ) GO TJ 27G 
FC.NF»FCNF + FG4*FC 
FTHNF«I-THNF*FG4*FTH 
FAJNF«FA 0NF4FG4+F4O 
FAlNF«FAl,NFtFG4*F41 
F31NF«FB1NFhFG4*FB1 
27C COST IN J E 
280 CONTINUE 

A00 THE CONTRIBUTIONS* CALCULATED FOR THE CURRENT t RAJcC TOR Y SEGMENT* 
TO THE COEFFICIENT MATRIX ELEMENTS 

DO 350 M* 1* NC OF 
MlMM-ll/NAZ + l 
M2«M-(M1-1)*NAZ 
I F ( M 1 .El. (NSP+1) 133 TO 370 
iM 3 * M 2 / 2 
M 5* M2-M 3*2 

IF ( M2 .EO. 1)A(L*M)«A(L»M)+FG(M1) 

I F ( .1 2 .GT. 1 .ANO. M4 , £ 0. 0 ) A ( L* M I « A ( L * M ) +FG ( Ml ) 

1*C0S(.M3*(P32*PB0 + D8) ) 

I F ( M 2 .GT. 1 .AND. ,N4 ,3T. 0 1 A < L * M ) « A ( L» M ) +FG ( Ml ) 
1*SIN(M3*(P82*PR0+D3)) 

GO TO 360 

370 I F ( M 2 .EO. 1 > A { L* M ) => A l L» M ) «-F THNF 
I F ( M2 • E 0 • 2)A(L*M)»A(L*M)'fF40NF 
I F ( M 2 .El. 3 )A(L*M)*A(l*M)+FA1NF 
I F ( M 2 .£3. 4)A(L*M)»A(L»M)+FB1NF 
360 CONTINUE 

B(L»1)*<J(L>1)-FCNF 

PROGRAM STEP 7. 

REDEFINE THE AZIMUTH INTERVAL. FOR THE NEXT TRAJECTORY SEGMENT. 

TEST TO SEE IF THE FINAL SEGMENT HAS SEEN CALCULATED. ALSO TEST 
TO SEE IF THE NEXT SEGMENT IS CLOSE TO A BLADE* IN WHICH CASE 
REDUCEO SPACING IS TO 3E USEO. 

P 32 « P B 1 

I F < P3 2 .LE. PLIM130 TO 230 
I F ( 0? 1 .EO. DP2IG0 TO 390 
I F ( J 1 .GT. IMINTGO TO 390 

IFU1 .EO. 0 .AND. PB2 .LE. ( PM IN ( J 1 H-DP 1 ) ) GO TO 390 

I F ( K1 .EO. 1 .AND. PB2 .GT. (PMI^(Jl)-QPl) ) GO TO 390 

IF (K1 .EO. 1 .ANO. PB2 .LE. ( PMIN { J 1 1-DP 1 ) J GO TO 400 

IFIP32 .GT. (-3.*0P1) .AND. I8L .EO. 1 ) PB1«PB1-DP2 
I F( °32 .LE. C— 3.*0P1 > .OR. I3L .GT. 1 ) P 81- P Bl-OP 1 
IFIK1 .EO. 0 .AND. P81 .LT. ( PM IN ( J 1 ) ^OPl )) PB1-PMIN ( J 1 ) +DP1 
GO TO 210 
380 P81-P9I-DP1 

GO TO 210 
390 P9WB1-DP2 
K l* 1 

GO TO 210 
400 PB1-PB1-0P1 

K 1* 0 
Jl-Jl+1 
GO TO 213 
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c PROGRAM STEP 3. END IF LOOP FOR NUMBER OF BLADES. 

W 200 CONTINUE < 

c SEt u, THE ELUENTS CORRESPONDING 10 THE * .OXIU.RT UNKNOWNS 

c 

M»n($P*NAZ + 1 
A{L»M)*AlL»M)-UT 

M ■ M ♦ 1 

A(L.N)»A(L»M)+MU*CP1 

A(L* M) ■ A(L»M )+RBO{ I ) * SP 1+ . 5*MU* ( 1.-CP2 ) 

au!m)«A(L»M)-RB 0 (I)*C?I-. 5 *M'J*SP 2 
L-L + l 

PROGRAM STEP 9. END OF COLLOCATION LOOP. 

110 CONTINUE 


SI? tt Srr5I E |x“; A ED’JAT IJN$ NEEDED TO CEOSE THE STSTEN. 

DO A 10 I*l> A 
do A20 m*i,ncof; 

A(L»M)*0. 

420 CONTINUE 
3 l L > 1 ) * C . 

I F ( I .cO. 

• EQ. 

.EO. 


IF( I 
IFt I 


.AND. Nl 
.AND. N2 
• ANO. N3 


.E 3. 
• £ 3 1 

. £ 0 < 


1) SO TO 430 
1) GO TO A AO 
1 ) GO TO 450 
1 ) GO TO A60 


500 

A90 

460 


520 

510 

A70 

530 


1 
2 

3 - 

I F < I .EO. A .AND. NA . = 3. 

I F ( I .LE. 2)11-1 
IF ( I .GT. 2)11*1-1 
DO 470 Ml* 1» N SP 
I F C I .GT. USQ TO A80 
DO 490 M2*1>NSP 
M-t M2-1)*NAZ+I1 

jp(J|I^^Ht^' 3 l)A(L*Ml*A(L»M)-DRlM 1 )*GF(Ml-l»M 2 )/(O.* 0 R<M 2 )) 

\f\tl \lu 3ip>itL^iTtI*^JlT-5Ji5i , »*6fcAx,rt2»/C6.*ORlM2» » 

IF(M2 !e3. M1)A(L»M)-A(L»M)-DR(M1) 

CONTINUE 
GO TO 470 
DO 510 M2-1.NSP 
M*IM2-1)*NAZ+I1 

IE!?; ; i s ;Ua"“«i^NT-;rN“i , ;iJ’«” R Nu 1 .i«;i>NA''( 6 ..D»(»2. 1 

GO TO 510 
I F ( M 2 . E3 ■ 

CONTINUE 

CONTINUE , 

GO TO ( 53G> 5A0»550»560) » I 
M*NSP*NAZ>2 

A IL» M) * ( l.-RO)* t l.-RO + RO ) 7 ( 3. *AR ) 


Ml)AlL,M)*A(L,M)-RBa(Ml)»DR(Ml) 
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a(L*l)«AR*CT/<N3*<l.-R01) 

THE ABOVE EQUATION EQUATES HE TOTAL LIFT DUE TO ALL THE BLADES, 
AVER AGE 0 uVEk THE AZIMUTH* TO THE THRUST COEFFICIENT. 

L*L*1 
GO TO 410 
430 M*N3P*NaZ+1 
A C L* M ) * 1. 

6 (L» 1 1* THC 

THE ABOVE EQUATION SETS THE COLLECTIVE PITCH TO THE GIVEN VALUE. 

L*L* 1 
GO TO 410 
540 M»NSP*NAZ+2 

A(L*M1 «( l.-R0)*(i.-R0*R0*ROl/{12.*ARl-2./GAMA 
G ( L * 1 ) *0. 

THE Act JVE EQUATION REPRESENTS THE ZEROTH HARMONIC COMPONENT OF 
MOMENT EQUILIBRIUM ABOUT THE HUB. 

L-L + l 
GO TO 410 
440 M«NSP*NAZ*2 
A ( L* M ) «1 , 

9 ( L* 1 ) » AO 

THE ABOVE EQUATION SETS THE CONING ANGLE TO THE GIVEN VALUE. 

L*L* 1 
GO TO 410 
550 M-NSP+NAZ+1 

A(L*M)* MU*(1,-R 01*11, -R0*R0)/(4,*AR) 

M»M*2 

A(L»M)«-a.-RO)*(l.-RO*RO*SO)/(6.*AR) 

B(L»1)«-Tr*MU*<RO*<1.-RO*RO>-4.M1.-RO*RO*RO>/3. 1 /(4.*AR) 

THE ABOVE EQUATION REPRESENTS THE FIRST HARMONIC COSINE COMPONENT 
OF MOMENT EQUILIBRIUM ABOUT THE HUB. 

L-L + l 
GO TO 410 
0 M*NS P*N AZ + 3 
A ( L * M ) * 1 . 

B ( L» 1 1 * Al 

THE ABOVE EQUATION SETS THE CYCLIC PITCH COEFFICIENT, Al, TO THE 
GIVEN VALUE. 

L-L + l 
GO TO 410 
560 M»NSP*NAZ*2 

A(L,M)«MU*(1.-R0)*(1.-R0*R01/(4.*AR) 

M»M + 2 

A(L,M)«-(1.-P0)*(1.-R0*R0*R0)/(6.*AR) 

8 < L » 1 1 » 0 . 

C 

C THE ABOVE EQUATION REPRESENTS THE FIRST HARMONIC SINE COMPONENT OF 
C MOMENT EQUILIBRIUM A30UT THE HUB. 
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L*L*1 
GO TO 410 
460 M*NSP*NAZ*4 
A(L,M)»1. 

S(L,l)*9l 

THE ABOVE EQUATION SETS THE CTCLIC PITCH COEFFICIENT, 31, TO THE 
GIVEN VALUE. 

L-L + l 

410 CONTINUE 

PROGRAM STEP 11. 

SOLVE THE SYSTEM CF SIMULTANEOUS EQUATIONS AND PRINT THE SOLUTION. 

CALL GELIM(NCOF,NCOF, A, 1,3, IP I VOT, 0, WK, IERR > 

I F ( I ERR .EQ. 1 ) GO TO 570 
L * 1 

DO 580 I * 1 , NS P 
GC( I)*3(t,l) 

L*L*l 

DO 530 J* 1, NHM 
GC(I,J)*3CL,1) 

L * L* 1 

GS(I,J)*3(L,1) 

L*L + 1 

580 CONTINUE 

THC*3(NSP*NAZ*1,1) 

THC0*THC*i3C./PI 

AC*9(NSP*NAZ+2»1) 

40D*A0*18C./PI 

A1*3!NSP*NAZ+3»1) 

A 10*4 1*130. /PI 
ai-B(NSP*NAZ+4,l) 

810-81*130. /PI 
WRITE (8,9) 

WRITc(6,590> 

590 F0RMAT(6X, "SOLUTION FOR COE FF IC I ENTS"/ 6X, 2 5 ( 1H- ) / /6X , 
l"(GO(I),I-i,NSP)"» 

WRITE (6, 600) (GO (I), 1*1, NS?) 

600 FORMAT!/ 6X, 5(E10«4, IX)) 

WRITE ( 6, 610) 

610 FORMAT! / /6X, " ( ( GC ( I , J ) , J - 1 , NHM ) , I -1, NS P ) "/ / ) 

DO 620 1*1, NSP 

WRITE (6,600) (GC(I,J),<I S 1,NHN) 

620 CONTINUE 

WRITE! 6,630) 

630 FORMAT! / /6X, " ( ( GS ( I , J ) , J - l, NHM ) , I * 1, NSP ) "/ / ) 

00 640 1*1, NSP 

WRITE(6,600)(GS(I,J),J*1,NHM) 

640 CONTINUE 

WRITE (6, 650)THCD, AOD, AID, 310 

650 F0RMAT(//6X, "PITCH ANGLE AT SLADE R DOT *", F 10. 5, 1 X , "DE GR E ES "/ /6X , 
1"C0NING ANGLE*",F10.3,1X,"0EGREES"//6X, 

2"FL AP P ING COEFFICIENT, A 1 F 10 . 5, 1 X, "DEGR EE S"/ / 6X, 

3 "FLAPPING COEFFICIENT, 3 1 *", F 10 . 5, 1 X , "DE GRE ES " ) 

C 

C PROGRAM STEP 12. 

C START LOOP FOR AZIMUTH STATIONS AT WHICH OUTPUT QUANTITIES ARE 
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C CALCULATED. 

C 

CTC AL*0. 

C NXC A L * 0. 

CMYCAL-C. 

00 670 1*1.26 
?BG0(I)*lb.*CI-l» 

P30»Pa0Q( I)*PI/180. 

cpx-casxpao) 

SPi*S IN ( PBO) 

cp2»cgs<2.*pbo) 

SP2»SIN(2.*PB0) 

F2*CPl*(2.+NU*TW*RQ/( 1 • -R 31 *-2 . *NU*THC ) +2 , + MU* AO* S Pl + 2 . *MU*31*C P2 
l-SP2*(MU*MU+TW/(l.-R0)+2.*iU*Al) 

F3»A0-CP1*<6.*MU*TW/<1.-R0) +2.*A1>-2.*B1*SP1 
GI1-0. 

G 1 2 * 0 . 

DO 690 1 1* 1» N SP 
G ( I. 1 1 ) *G0 ( 1 1 ) 

00 690 12* 1. NHM 

G(I»I1)*S(I.I1)+GC(I1»I2>*C0S(I2*PB0»+GS<I1.I2)*SIN(I2*P90) 

690 CONTINUE 

DU 700 I1*X»NSP 

IFTISEL .NE. 0) GO TO 7X0 

GI1*GI1+0R<I1)*G<I»I1> 

GI2*GI2+0R(I1)*RB3<I1)*G<I. II) 

GO TO 700 

710 GI1*GI1+6.*DR< U)*G< X. ID 76. 

GI2*GI2+6.*0RTI1)*R30(I1)*G(I,I1) 76. 

I F ( I X ,EQ. X ) GL ( 1 1 ) * 0 . 

I F ( X 1 .GT. 1 ) GL ( 1 1 ) *GR ( 1 l-l ) 

GRt 111*0. 

I F ( I X • E 0 • NSP ) GO TO 720 
DO 730 12*1/ NSP 

GR(I1)*GR{I1)+GF(I1»I2)*G(I»I2)/DR(I2) 

730 CONTINUE 

720 GIX*GI1+0R(II)*(GL(I1)+GR(I1))76. 

GI2"GI2+0RCIlJ*(R{Il)*GL(ll)+R(Il+l)*GRCIl))/6. 

700 CONTINUE 

FLTim*(l.-R0*RQ)*A072. + CPl*(2.*l1U*THC*<l.-R0J-2.*MU*T* 

1-A1*( l.-RO+RO) ) +SP1*(1.-R0»*(2.*NU*A0-81*( l.+RO) ) 
2+2.*MU*Bl*(l.-RO>*C?2-$P2*(2.*MU*Al*Cl.-RO>+;iU*NU*TW> 

FLT1 ( I ) *-FLTU I )*( l.-RO ) / (4.* AR ) 

FLT1( 1 1*- <FLT1( I)+GI1J*PI*C l.-RO) /AR 
FLT2(II»FLT1CI)/(CT*PI/N8) 

FMT( I)»Cl.-R0*R0*R0)*A3/3.+CPl*(NU*THC*<l.-R0*R0)-2.*Al* 
1<1.-R0+R0*R0) 73. +MU+TW* ROM l.+RO )74.*MU*TW+{1.+R0+R0*R0) 73. ) 

2 + SPlMNU*AO*(l.-RO*RO)-2.*31Ml.-RO*RO*RO)/3.)+CP2*HU*Bl* 

3( l.-RO* RO ) -S P2* ( NU* Al* ( l.-RO+RO) + HU*MU*TWM l.+RO ) /2. ) 
FHTm—FHTm*Cl.-RO)/(*.*AR) 

FMT( IJ*-(FNT( I)+GI2)*PIM l.-RO) /AR 
C0L<I)*FMT(I)/FLT1(I) 

CTC AL*C TC AL+FLT1 (1)726. 

CNXCAL*C.NXCAL+Ft1T{I)*$P17 26« 

CMYC4L*CNYCAL-FMT< D + CP1726. 

START LOOP FOR RADIAL STATIONS AT WNICH OUTPUT QUANTITIES ARE 
CALCULATED. 

JO 760 II-l.NSP 
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1 ? i > •/*;»** J* * 

“ G 5;;U.“U»lo(n!.ic«dO(ni.«»o(ui 
GO TO 790 
780 GGUT*G( I# ID 

7,0 ■EE!cililU-5t5!Si“o»*?GO.J T -«».-*OI*FS3«I.Il»/C^.**R»W*R 

FL ?t I ;^ 7 'pTia'Uo)MD”o ! !*ll.-RO)*f23(DlD/(16.*4R*AR*4«) 

740 CONTINUE 
670 CONTINUE 

CTCAL-CTCAL+N8/PI 

CMXCAL»CNXCAt*N0/PI 

CMYCAL*CMYCAL*N3/PI 

( *ftITcI6»791)CTCAL»CNXCAU» CrtYCAL enT , h , E10 .<,/ / 6X, 

791 FQRMAT( /6X»"C0NPUT E _ i^OUT ROTOR X-AX 1 5 »" * E 10 • *» 1 1 6X ' 

SSSSB KSS ISfEKKK BX! *o,o» 

PRINT ALU OUTPUT QUANTITIES IN TA8ULAR FORM. 

WRITE(6,9) 

800 illl !1(/^;:«TA3LE 1 - SECTIONAL L I FT , < RHO* ( OME 3 A **2 ) * ( R 1~ 3 , > «, 6X 

1 MRITE(6l310»(RaC(I)fI»i»i«SPI 
310 F0RNAT(/12X,"R/SD ".»( tlO.** M> > 

WRITE(6*8il) 

811 F0RMAT(/7X,"PSI") 

»RlTi?6»330)P800(I»»(FLlII*ID»Il m l»^5*** 

330 F 0 RMAT(/ 6 X»F 5 « 1 * 8 X, 3 (E 10 . 4 » 1 X)) 

820 CONTINUE 

WRITE (6*9) 

;55;!;; , ,5"!-T..LE * - SSCTIOIUI UfT.RHTH.vm PER EE.DE-/6P. 

1 ypiT£tE»aloltR0o(l)»I-l»RSP> 

W RITE ( 6 , 311 ) 

WRIT£?6* 8 3CH P800(I)*(FL2(I*IlDll m l*NSP) 

850 CONTINUE 

WRITE ( 6* 9 ) 

WRITE(6,36G) SECTIONAL PITCHING MOMENT /( RHO* ( QrtEGA**2 > * ( 

WRITE(6,810)IRB0(I)»I*1*NSP) 

WRITE(6,8ll) 

WRITE ( 6 ^ 330 IP 300 ( 1 )* (FM( 1 , 11 > , 1 1* 1» NSP ) 

870 CONTINUE 

WRITE (6,9) 

8S „ ;si;svt*/?«:-T..LE * - ce«« * «5»»« fRon 60 

10E(FRACTI0N OF CHORO)"/6X,M(1H-) //> 

WRITE(6,810) (R30(I)*I-1»NSP) 

WRITc(6,811) 


82 



o o o o non 


00 S90 1-1*24 

-RITE(6,930)PaG0(I),<XCP(I»Il)*Il-l*NSP> 

890 CONTINUE 

w»ITE(6,9) 

W° I TE ( d, 90C ) 

9C0 F3RMAT( //6X. "TABLE 5 - TOTAL BLADE LIFT* MOMENT ABOUT HUB AND RADI 
1AL CENTER OF L I FT"/ 6X , 70 ( IH-) //6X * 

2" TOTAL BLADE LIFT/ ( RHQ*{ 0MEGA**2) *( Rl**4) ) "/6X, 

3"T0TAL BLADE LIFT/THRUST PER BL ADE"/6X* 

4"M0MENT ABOUT HUB/ ( RHO* ( 0MEGA**2 ) * { R 1**6 ) ) " /6X, 

5"RADIAL CENTER OF LIFT/R1"//7X, "P SI "* 8 X* "TO T AL BLADE LIFT". 
55X*"NO.MENT",5X* H CENTER"/39X,"AdOUT HU3",4X,"DF LIFT"! 

DO 910 1*1,24 

WRITE ( 6* 920) P 800 ( I)*FlT1<I)*FLT2(I)*FMT{I)*C0L(II 
92U FOkMAT (/6X*F5.1*5X»4( E10.4*1X) ) 

910 CONTINUE 

WRITE (6, 9) 

-RITE (6,940) 

940 FORMAT (//6X, "TABLE 6 - SURFACE PRESSURE DIFFERENTIAL/ (RHO* (OMEGA** 
12)*(RI**2) )"/6X,64(lH-) //) 

00 930 1-1,24 
WPITE(6,960)P80D(IJ 

960 FORMAT! /6X, "AZIMUTH ANGLE •"» F 5 . I* IX* "DEGRE E S"/6X* 27 ( 1H-) ) 
WRITE(6,310) (RBO(Il) , 11*1, NSP) 
wRI TE ( 6* 970) 

970 FORMAT! 11X*"X/C") 

DO 980 11*1, NSP 
DO 980 12-1,10 
CCH-2.*X0UT( 121-1. 

SCH-SQRT( l.-CCH*CCH) 

P0UT( II, I2)*-G( I, Il)*SCH/( l.+CCH)»F23( I* II ) *{ l.-RO) *SCH/ ( 2. * AR ) 
POUT (I 1,12) -2.* POUT (II, 12) 

980 CONTINUE 

00 990 12-1,10 

WRITE(6,99UX0UT(I2),(P3UT(Il,I2),Il»i,NSP) 

991 F0RMAT(6X*F8.5,5X,5(E1J.4»1X)) 

990 CONTINUE 

IREN* I— ( I /3) *3 
I F ( IR EM .£9. 0) WRI TE ( 6* 9 ) 

950 CONTINUE 
STOP 

PRINT ERROR MESSAGE IF COEFFICIENT MATRIX IS SINGULAR. 

570 WRITE ( 6,930) 

930 F0RNAT(6X, "COEFFICIENT MATRIX SINGULAR") 

STOP 

END 

SUBROUTINE DM IN (MU, L AM, 03*R80, P 80, PL IM, DM A X* I, P ) 

CALCULATION OF AZIMUTH POSITIONS AT WHICH TRAJECTORY IS DIRECTLY 
OVER A 3LA0E, WITHIN A DISTANCE OMAXj 

REAL MU, LAM 
DIMENSION p ( 20) 

Y(X,X1,X2,Y1,Y2)*Y1*( Y2-Y 1 ) * ( X-Xl ) / ( X2-X1 ) 

I-D 

Pl-O. 

XBI-O. 

P2--0.2 
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10 X02«R3O*SIN<P2*Oa>*NU*SIN{P2*P3O*DB)*P2 
IF<P1 • NE • O.JGO TO 20 
30 P1-P2 

I F ( P 1 .Lc. PLIMJRSTJRN 
P2-P2-0.2 
X 3l*X 82 
GO TO 10 

20 TEST*X61*X32 

IFtTSST .GT. O.JGO TO 30 
PC«r<0.,X81»X82,Pl»P2) 

X3C-R30 + SI.NIPC + 03J MU*SIN( PC*P30*D8 J*PC 
R6C»R30*CuS(PC«-03J *MU*COS ( PC+P30+DB J -PC 
Q-S3RT(X3C*XBC + LAN*LA.N*PC-PC+(P. 8C-R80J*(I<3C-RB0) ) 

I F ( 0 .GT. DNAXJGO TO 30 
1 - 1*1 
P ( 1 1 - PC 
GO TO 30 
ENO 

subroutine subivl(rbi,r32,psi>pb2»r,i*pi»p2> 

DIVISION OF TRAJECTORY SEGMENT INTO SUB-SEGMENTS ALIGNED WITH BLADE 
SPANWISE SEGMENTS. 

DIMENSION k(6J, I (5) . P 1C 5 > , P 2 ( 5) 
Y(X»X1»X2*Y1>Y2J-Y1*(Y2-YIJ*(X-X1)/CX2-X1) 

DO 10 J -1» 5 
I(J J-0 
P1UJ-0. 

P2( J ) -0. 

10 CONTINUE 

DO 23 J * l» 5 
Rl-RCJJ 
R 2- R ( J*l) 


I F ( R3 1 

.LE. 

R 1 

.AND. 

R32 

• LE. 

R1 J GO 

TO 

20 

IFCRfll 

.GE. 

R 2 

.AND. 

RS 2 

.GE. 

R2JG0 

TO 

20 

I F ( R 3 1 

.GE. 

kl 

• AND. 

R31 

.LE. 

R2 .AND. 

R32 

1 • AMD. 

R82 . 

LE. 

R2JG0 

TO 

30 




I F ( R3 1 

.LE. 

R 2 

.AND. 

R32 

.GT. 

R1JG0 

TO 

40 

IFCR31 

.LT. 

R2 

• AND. 

R32 

• Gc • 

R2J GO 

TO 

50 

1 F ( R3 1 

.GT. 

R1 

.AND. 

RB2 

.LE. 

R1JG0 

TO 

60 

IF( R31 

.GE. 

R2 

.AND. 

R3 2 

.LT. 

R2JG0 

TO 

70 


30 KJJ-l 

PKJJ-PBl 
P2(J)-PB2 
RETURN 
40 I(JJ«1 

P1(J)-Y(R1»R81»RB2»P31>P321 
I F ( R32 .GT. R2)P2(JJ-Y(R2»R31»RB2>P31»PB2J 
I F ( R32 .LE. R2)P2( JJ-P32 
GO TO 20 
50 KJJ-l 

P2(JJ-Y(R2»R81»R32#PB1#P82) 

I F ( R 3 1 .LT. R1JPK JJ«Y(R1»R31#RB2»PB1#PB2J 
I F ( R3 1 .GE. R1JPKJ J-P31 
GO TO 20 
60 KJJ-l 

P2(JJ-YIR1»RB1»R32»P81»P32J 
I F ( R 3 1 .GT. R2)P1(J!-Y(R2»K81»R32»PB1>PB2) 
I F ( R31 .LE. R2JPKJJ-P31 
GO TO 20 
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70 I ( J ) « 1 

Pl< J J»Y(R2,R81f RB2,PB1»PB2) 

IF(R3 2 .LT. R1)P2(J)*Y(RI»RB1»RB2»PB1,PB2) 

IF ( R82 .GE. R1)P2(JJ«PB2 
20 CONTINUE 
RETURN 
ENO 

SUBROUTINE NF IN T< P i» P2» < . ISEUT 1, T2» T3 ) 

INTEGRATION OF NEAR FIELD PRESSURE GRADIENT. 

DIMENSION F0(2)»Fi(2)»P2(2)»F3(2) 

C0M10N/HAIN1/XI»XS»YS»RI,RS 

oo»tx-xn*(x-xn 

01»-2.*(X-XI)*XS 

02*xs*xs+ys+ts 

Q«A.*U0*i)2-Dl*Di 
IF<3 . NE. 0 . ) SQ*SaR T ( a ) 

00 10 I* l# 2 
I F ( I . EO. DP.PI 
I F ( I .EO. 2 » P-P2 
00*00 + 0 l*P + 02 + P*P 

I F ( J .NE. 0.)FC(I)-2.*ATAN{(01+2.*D2*P)/SQ)/SQ 
I F ( 0 .EO. 0.)F0(I)«-1./(02*P» 
F1(I)»(ALCG(D0)-01*F0(I) )/(2.*02> 

IF ( ISEL .EQ. 0) GO TO 20 
F 2 ( I ) - C P-01+FH I)-00*F0< I) >/D2 
F3(I J*( ?*P/2.-01*F2(I)-30*Fl(II )/02 
GO TO 10 
20 F 2< I ) *0 • 

F3( I) *0. 

10 CONTINUE 

T1*(X— XI) + {F0(2l-F0(l))-XS*(Fl(2>— F1CU) 
T2*(X-XI) + (F1(21— F1(1D-XS*(F2(2)— F2C1D 
T3- (X-XI»*(F2 (2 »-F2(l))-XS*CF3( 2) -F3( ID- 

RETURN 

END 

SUBROUTINE CPINT(P1»P2*1S£L»T1,T2»T3) 

INTEGRATION OF COMMON PART PRESSURE GRADIENT. 

OIMENSION F0(2),Fl(2)»F2(2),F3(2),F4<2) 

COMMON /MAIN1/XI»XS»YS,RI»RS 

RO* XI*X I 

R1-XI*XS 

R2«XS*XS+Y3+YS 

0«4.*XI*XI*YS*YS 

IF ( 3 .NE. 0.)S0*SQRTCQ» 

00 10 1*1,2 
I F ( I .EQ. 1JP-P1 
I F ( I .EQ. 2) P*P2 
RR»R0+2.*R1*P+R2*P+P 

I F ( Q . EQ. 0. .AND. P .EQ. O.JGJ TO 20 
I F C Q .NE. 0,)T*2. + ATaM( 2.*R1+2.*R2*P)/SQ)/SQ 
I F ( a .EO. 0. J T* 0. 

IF ( Q .NE. 0. »F0(I)*2,*U1+R2*P> / I 84RR ) +2. * R2+T/3 
I F ( Q .EO. 0.)FG(IJ»-l./{3.+R2+R2*P*P*P) 
F1(H«<-1./(2.*RRI-R1*F0(I))/R2 
F2(I)*(-P/RR + R0*F0(I) I/R2 
I F ( ISEL .EO, 0) GO TO 30 
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F3( I )«{-R0*Fl(I)-2.*Rl*F2( I) ) /R2* ( ALOG ( PR ) 12 . -R i * T ) / < R 2 * R2 ) 
IFIQ .EQ. 0.)F3(I)«AL0G(P*P»/<2.*R2*R2) 
FA(I)«(P*P*P/RR-A.*R1*F3I I )-3 . *RO*F 2 ( I ) I/R2 
G3 TO 10 
20 Fom«o. 

F 1 { I ) »0 • 

F2 ( I ) *0. 

30 F3(I)»0. 

FMIJO. 

10 CONTINUE 

DF0«F0(2)-F0(1) 

0F1-F 1 C 2 » -FI C 1 » 

DF2»F2(2)-F2(1) 

DF3-F3C 2I-F3C X» 

CFA«FA<2)-FAI1> 

T 1»R0+3F0+2.*R1*0F1M XS* X 3-YS* YS ) *DF2 
T2«R0*0F1 + 2.*RI*0F2MXS*X$-YS*YS)*DF3 
T3»RO*0F2 + 2.*Rl*DF3f( XS*X$-YS*YS )*0FA 
RETURN 
END 

SUBROUTINE FFINT(P1»P2»R* ISEL»T1» T2*T3) 

INTEGRATION OF FAR FIELD PRESSURE GRAOIENT. 

DIMENSION Fl(2)»F2(2)»Fi(2) 

C0*H0N/HAIN1/XI,XS»YS>R I»RS 

R0«XI*XX 

R1«XI*XS 

R2»XS*XS+YS*YS 

D0»R0+(R-kI)+(ft— RI) 

Dl« Rl- ( P-R I ) *R S 

D2»R2+RS*RS 

S02*S0RT(D2) 

FN1*A9S (XI*YS > 

FN2*(R-RI)+R2+RS*R1 
FN3»Dl*Di-00*D2 
Q“4.*02*FN1*FNI 
I F ( 3 .NE. 0. )SQ«S3RT<0) 

DO 10 1 - 1,2 
I F C I .EQ. 1>P«P1 
I F ( I .EQ. 2) P*P2 
F 1 ( I ) *0 . 

F2( I ) *0. 

F3(I)aO. 

xa«xi+xs*p 

YB«YS*P 

R8«RI+RS*P 

RR»XB*XB*YB*Y8 

U*R-RB 

D«SQRT(RR+U*U) 

T«Dl + D2*P-*-S02*0 
G1«AL0G(T*T)/(2.*S02) 

G2«(D-D1*G1> /D2 

I F C X I .EQ. 0. .AND. P .EQ. 0.1G0 TO 20 
I F ( X I . LQ . O.JGO TO 33 
IF ( IS EL .EQ. 0) GO T3 AO 

G3-ATAN(2.*(T*(S02+RS)-FN21/SQl-ATAN(2.*(T*<iD2-RS>»FN2) /SQ) 
G3"A.*SD2*G3/SQ 

G A“ ( ALOGI ( D-U ) / ( 0 U I ) -2 . *RS*GI-R1*G 3 ) /R 2 
G5»(-R0*G3-2.*R1*5A+-2.*CR-RI)*G1-2.*RS*G2)/R2 



G6»< 8.*S02*FN3MSQ2*D-RS*U) / ( T*RR ) ♦8.*S02*FN2 + U/ RR 
H-A.»02*FN2*G3)/Q 
G7«(-2.*0/RR-Rl*G6-RS*G3> /R2 

Ga-(-R0*36-2.*Rl*67f2.*3H-2.*(R-RI)*G3-2.*RS*G4WR2 
G9«(-R0*G7-2.*R1*G3+2.*S2*2.*<R-RI»*G4-2.*RS*G5) /R2 
Fl< I)«<2.*GU(R-RI)+03-RS*G4-YS*YS*G6)/2. 

F2( I ) »(2.*G2 +IR-RI )*GA-RS*G5-YS*Y S*G9) /2. 

G13»2.*ATAN( (2.*CSD2-RS)*T*2.*FN2)/Sa)/Sa 
Gll*AL0G(ABS(<SD2-RS)*T*T+2.*FN2*T-<S02+RS>*FN3n 
G12»(-2.*RH'SD2*G10+Gll-< S02-RS)*G1)/R2 
G13»(P+2.*RO*( XS*XS-YS* YS ) *S02*G10/R2-2. *R 1+G11/R 2 
1 ♦ ( 2. *K l*D2-( SD2+RS)*FN2)*31/(D2*(S02 + RSn+RS*D/02)/R2 
F3( I)»P*ALGG(0 +U)-k1*G12-<5*XS*G13 + RS*G2 
GO TO 10 

AC Fim»RS*YS*(YS*0 + SD2 + Y3)/(R2*SD2*T*D)-X3*XB* 

1( kS*0-SD2*UW(R2*T*RR)-Y3*Y8*U*(l.-Rl/T)/I J*K2*RR» 

F 2 C X » -0 . 

F 3 ( I ) *0 . 

GO TO 10 

30 IF( ISSL • E Q. 0)G0 TO 50 

G1A«(T-01-SD2*ABS<R-RI> ) /( T-01*S02*ABS l R-R 1 1 ) 
G1A«ALuG(G1A*G1A)/2. 

F 1 ( I)» (XS*XS-YS*YS)*(-0/P+D2+Gl+01*GlA/ABS (R-R I) I / ( R2*R 2 ) 
1+ YS *YS *G1 /R2 

F2(I)»(X3*XS-YS*Y3)*(0+AB3(R-RI)*G1A*01*G1)/(R2*R2> 

i+YS*YS*G2/R2 

G15»(P+RS*0/02) /R2-(H-RI)*Gl/02 
F3(I)«P*ALGG(D+U)-XS*X3*G15+RS*G2 
GO TO 10 

50 Fl( I ) * ( XS*XS-YS+YS)*0/(R2*R2*(R“RI)*P) 
1+YS*YS*P/(R2+(R-RII*0) 

F2(IJ»0. 

F3(I)-0. 

GO TO 10 

2C IFIISEL .EQ. OJGO TO 60 
G16-AL-3GI T/(2.*U*J) ) 

Fit I) «< XS*XS-YS*YS)*(02*G1-01*(1.-G16) /ABS<U))/(R2*R2J 
l+YS*YS*Gl/R2 

F2<n-(XS*XS-YS*YS I*(01*G1 + ABS(U»*(1.AG16) )/(R2*R2) 
1+YS*Y3*G2/R2 
G17*RS*D/(R2*D2)-U*Gl/02 
F3( I)»-X$+XS*G17*RS*G2 
GO TO 10 

60 FI d ) •-( XS*XS-YS*YS) *RS/(R2*R2*ABS(UH 
F2( 1 1 *0. 

F3( I ) »0. 

10 continue 

T1«F1C2)-F1(1) 

T2*F2(2)-F2(1) 

T3«F3(2>-F3(1) 

RETURN 

END 

SUBROUTINE TABSCH(X,N#XT»I1»I2# INT) 

C 

C GIVEN AN ARRAY X, TO LOCATE THE POSITION OF A VALUE XT. 

C IF I NT »Q » xr LIES BETWEEN XIII) A.NO X(I2». 

C IF INT*1» XT IS GREATER THAN X(N). 

C IF INT — 1» XT IS LESS THAN X(l>. 

C 


DIMENSION X f N 1 



I1«0 

12-0 

NM « N-i 

00 10 I*l#NN 

IFUT .Gt. XII) . AND • XT .Ls. X(I*1>)G0 TO 20 
10 CONTINUE 

I F ( XT .LT. XIIHGO TO 30 
IF ( XT .GT. X ( N H GO TO AO 

20 INT-0 

IFIXT .EO. XIIHGO TO 21 
IFIXT .60. XII + 1HG0 TO 22 
11*1 
12 * 1*1 
RETURN 

21 11 * 12*1 
RETURN 

t 22 11 * 12 * 1*1 

RETURN 
30 INT — 1 
RETURN 
40 INT*1 
RETURN 
PNO 
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